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Preface 

This user and programmer reference manual together with the Progress 
Report for the Period Ending 31 October 1977 constitute the Final Report 
for Contract NAS5-20974. This manual describes the computer programs 
developed to implement the ATS-6 navigation model on both the AOIPS at 
G3FC and the McIDAS at SSEC. Tlie above mentioned progress report 
describes the theory of the navigation model and the results obtained 
with it. Some of the Appendices of that report have been updated and 
included here. 
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I. The ATS-6 Navigation Problem 
A. Introduction-Objective 

If a satellite image or sequence of images of the earth is 
to be useful for quantitative measurement, it is necessary to be 
able to determine the earth location (latitude and longitude) 
that corresponds to any given picture element (line and element) 
in that image or sequence of images. Development of a navigation 
system for a given satellite involves two procedures: 1) Defining 

an algorithm for converting a satellite picture element location 
to earth location and vice versa; 2) Defining a procedure for 
measuring the set of constants needed by the algorithm in 1) above. 
Previous progress reports for this project describe how the ATS-6 
navigation model was developed. This user manual briefly describes 
the current version of the navigation model (this section) and 
how to use the computer programs developed for it (following sections) 
In the process of developing a navigation model for satellite 
images we must establish criteria of validity for the resulting 
model* Since there is no way to precisely relate the image 
location to other satellite sensors, the criteria of validity 
must be referenced to some measurement derived from the image 
itself* For this work, two tests were used* First, the line and 
element position of identifiable earth points (landmarks) was 
measured and compared to the values predicted by the navigation 
model* Ideally these should always differ by less than one pixel. 

c 

In general, this navigation model is capable of predicting positions 
within an error of the order of one pixel. Second, a set of cloud 
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tracer winds was derived from a three-image time sequence* These 
winds were compared to measurements made on the same clouds in 
SMS“1 images from approximately the same time* Winds, averaged 
over approximately one hour, from each of the two satellites 

t 

showed agreement on the order of one meter per second. In the 
< following sections we describe the model used in making these 
measurements. 

B. Nature of the Problem 

A three-axis stabilized satellite contains an attitude 
measuring and correcting system which attempts to keep the satellite 
and its camera pointed at some selected point. For ATS-6 the 
camera is generally aimed at the point on the earth equator 
which is closest to the sub satellite point. In order to keep 
the satellite camera pointing in its desired direction within 
its requir^d~range of accuracy, the satellite's attitude may be 
changed by significant and mp re die table amounts about all three 
axes several times during the scanning of an image. If wind speeds 
are to be measured accurately, the attitude changes must be measured 
and accounted for. 

In addition to attitude changes of the satellite, we find that 
we must also account for image distortions caused by the scanning 
mechanism. We refer to this problem as mirror-scan nonlinearity. 
Ideally the images would be generated by sampling at equal-angle 
intervals (angular position a linear function of picture element 
number) as the camera^ s mirror scans across the earth. However, the 
ATS-6 camera scans in both left to right and right to left directions 
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and therefore the camera^ s scanning mirror- must change direction 
and accelerate for each scan. Two resulting effects are evident 
in the images: scan lines of opposite direction are offset by 

7 to 11 pixels and each scan line varies from equal- angle sampling 
by 2 or 3 pixels maximum over a 200 pixel range. 

C. The Solutions to the Problems 

In this section we describe the methods currently used for 
attitude determination and for correcting mirror-scan nonlinearity 
in ATS-6. 

The satellite has attitude sensors and data from these sensors 
are recorded on magnetic tape along with the image data. Unfortun- 
ately, we found that these attitude telemetry data did not accurately 
reflect changes in satellite attitude as seen in the images. As a 
result we developed a method to correct for attitude changes using 
the images themselves. This is the earth edge displacement technique 
described in Section III of reference 3. These earth edge displace- 
ment measurements enable us to compute the changes in attitude with 
time. We must also compute a reference attitude, usually the mean 
attitude for the first image in a sequence of three to be studied. 
This reference attitude ^is computed from landmark measurements; the 
technique is described in references 1 and 2. 

The problem of mirror-scan nonlinearity cannot be solved 
completely, but a correction scheme adequate for our purpose of 
cloud tracking has been developed. Ideally, we would like to 
define some reference coordinate system to use in converting the 
images to equal-angle sampling. Unfortunately, we have no way 
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to define such a coordinate system which would be independent 
of the images. Thus, we use one scan direction (odd numbered 
scans) as a reference and shift scans of the opposite direction 
to match* This method corrects the alternate scan offset but 
leaves an uncorrected nonlinearity of up to two to three pixels 
at some points in the image. This error is reasonably constant 
from one image to the next and is small enough and varies slowly 
enough that the images may still be used for cloud tracking. 
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II. Using the ATS-6 Navigation System 

A. Introduction 

This section gives a general discussion of the use of the 
ATS-6 image correction and navigation software* For specific 
commands and data inputs to the various programs see Section III. 

The general organization of the main programs is indicated in 
Fig. II. 1. 

B. Image Correction 

The alternate scan element shift function AE(E) which is used 
to correct for mirror-scan offsets (Appendix D) is computed from 
infrared image data. The program OFSTG computes values of AE and 
a weighting function at points at equal intervals across the image. 

This procedure is described in detail in Appendix D. Since the 
image is not viewed before running this program it is not knom if some 
of the computed offset points will be computed for locations off 
the earth. Some offset points will be averages of computations 
from data entirely on the earth, some entirely off the earth and 
some mixed. Because the brightness off the earth is almost uniform, 
the weighting factor, which is the brightness range in the line 
segment used in the correlation computation, will cause the off- 
earth points to have negligible effect on the averaged offset. Data 
points for correlations done entirely from off-earth data must be 

eliminated from the curve-fitting process. To do this the user 

— y ' * 

must scan through the listing of offset and weight values from 
OFSTG. There is a fairly sharp discontinuity at both earth edges. 

The element location of these edges are then input to program OFSTF. 
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Figure II. 1. General organization of ATS-6 iiaage correction and 

navigation software. Names of corresponding computer 
programs in brackets. 
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The offset data points generated are not only discontinuous 
but also very noisy. Therefore we fit a smoothing and interpolating 
polynomial to the data points. This is done by program OFSTF. As 
stated above OFSTF must have the limits of valid data. The program 
OFSTF sets up arrays for subroutines APCH and APFS. These least- 
squares polynomial fitting subroutines are standard procedures 
from the IBM Scientific Subroutine Package (SSP) . 

Once the alternate scan correction polynomial has been deter- 
mined, program LATSP may be used to read image segments from the 
tapes. This program requires input data specifying the date and 
time of the image for image identifier entries. This program reads 
a standard size image area which is 512 lines by 512 elements on 
AOIPS or 500 lines by 672 elements on McIDAS. The location of this 
segment with respect to the whole image is specified as input data 
in terms x>f line and element of the upper left comer of the image 
segment. 

It should be noted here ' that ‘ this alternate scan correction 
method only shifts data from one scan direction to match that of 
the other. It does not correct the data to true equal^angle sampling. 

C* Attitude Determination 
/ 

The initial satellite attitude is determined from landmarks 
and the satellite’ s orbit. The attitude computation is done using 
several widely spaced landmarks on the first image (t^) of the 
sequence to be used. A minimum of three landmark areas should be 
used. One or, preferably, more landmark points may be measured in 
each area. Landmarks should also be measured for other images in 
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the sequence to serve as a check on the computations and on the 
stability of the attitude* 

Attitude changes from one image to the next are determined 
from earth edge displacement measurements. These displacements 
are measured using the regular AOIPS or McIDAS cloud tracking 
program (WINDCO) . For this procedure a pair of images at the 
earth edge are loaded (say t^ and ^ 2 ^^ Then the wind tracking 
procedure using a correlation tracking metric is used to track 
the displacements of the edges* The line direction lag size is 
set to zero so that the correlation peak search is done by moving 
the target grid only horizontally (i.e* in element direction) 
within the search grid- This measurement is made approximately 
every 5 to 10 lines along the edge. It need only be done for the 
same range of lines as is covered by the area to be used for 
cloud tracking; but must be done for both edges* As in the case 
of the alternate scan data^ these measurements form a discontinuous 
and slightly noisy function. The program EDGFT is used to fit a 
set of Chebyshev polynomials to the data points for each edge* 

The coefficients, number of coefficients, scaling factors and 
valid range of use in terms of time of day at a ^given line are 

^ stored in common block NAVCOM. Reference 3, especially Section III, 
gives more of a discussion of the earth edge correction techniques. 

D. Measuring Cloud Tracer Winds 

Once the attitude and earth edge displacement polynomial coefficients 
have been computed and placed in common block NAVCOM along with 
the orbit and frame constants, the coordinate transform functions 
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are ready for use in cloud tracking. The ftmctions available are 
SE for satellite image coordinate (line, element) to earth 
coordinate (latitude, longitude) transformations and ES for the 
inverse. For compatibility with SMS systems a subroutine SATEAR, 
which links to ES and SE is supplied. 





III* Program Command Format 

This section describes the command formats for using the ATS-6 
navigation programs* There are two versions, one on the AOIPS and 
one on the McIDAS. Part A covers the use of the AOIPS version. 

AOIPS uses a prompting system; the appropriate responses to each 
of the prompting requests are explained. McIDAS uses a command 
input system; the command sequence and input parameters are explained* 

A. AOIPS Version of November 1977 

There are currently six main programs plus the main menu 
driver for ATS- 6 processing on AOIPS. These are set up to 
work with the METPAK driver MET2 ox Terminal 2. The system 
is started up by mounting the METPAK disk on BD0:, then installing 
the tasks* The main menu driver expects them to have a name 
ending in 2 (e.g* LATSF2)* The installation command file 
RD0 : [1,162 ] ATS INS* CMD is available* Operation must be initiated 
by first running METPAK to restore the global common* Once the 
common has been restored, the EXIT option is taken and the ATS--6 
menu driver ATSF2 is initiated (i.e* RUN DB0 : [350,62 ]ATSF2) * The 
system Is now ready to run- 
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A, AOIPS Version, of November 1977 

1. A6INT - Initialize COMMON/NAVCOM/ 

Initiate by requesting option: 1 - INITIALIZE NAVCOM 

on the ATS-6 PROCESSING menu. 

Input request : 

NAVCOM TO BE PRINTED? (1 = YES, DEFAULT = NO) 

Response ; If a 1 is entered, the updated contents of 
COMMON/NAVCOM/ will be output to the line printer just 
before A6INT exits. Otherwise, not. 

Input request : 

ENTER DAY NUMBER 

Response : Enter the Julian day number (1-363) of the data 

to be worked with. Year is assumed to be 1974. 

Input request : 

ENTER 3 PICTURE START TIMES - HHMMSS 

Response ; Enter the picture start time, to nearest second, 
for the three images to be worked with. As currently set 
up COMMON/NAVCOM/ can only hold enough coefficients for 
earth edge correction for three image times. 

Input request ; 

ENTER FIRST ORBIT POSITION: T(HHMMSS), X, Y, Z (KM) 

Response : Enter a satellite orbit position vector as read 

from orbit data on the experimenter history tape. T is the 
time of the position in packed integer format (hours, minutes, 
seconds), X, Y, Z are the location in the earth inertial 
co-ordinate system in floating point kilometers. 

Input request : 

ENTER SECOND ORBIT POSITION: T(HHMiBS), X, Y, Z (Ipl)- 

Response : Same as above for a second satellite orbit position. 
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Default ; Except as noted, the default response causes a 
return to the ATS-6 processing menu. 

Program output ; A6INT will enter values into COMMON/EAVCOM/ 
which is part of the saved global common. 

2. OFSTG - Generate alternate scan offset data points. 

Initiate OFSTG by requesting option: 

2 - GENERATE ELEMENT OFFSET DATA FROM E.H.T. 
from the ATS-6 PROCESSING MEND. 

Set up : Mount an ATS-6 Experimenter History Tape on tape 

drive MM0 or MML . 

Input request ; 

MOUNT TAPE. ENTER DRIVE NUMBER. 

Response ; After mounting the tape, enter a 0 or 1 for MM0 
or MML respectively. 

Default ; The default is tape tmit MM0. 

Program output ; OFSTG will print a table of element numbers, 
offsets and weights on the line printer. These same values 
are stored in global common in COMMON/BUFFER/ and COMMON/BUFFI/ 
for use by program OFSTF. 

Note ; This program is a heavy user of CPU cycles. 

-v/ 

s 

3. OFSTF - Fit a set of Chebyshev polynomials to the data points 
generated by program OFSTG. 

This program is automatically initiated at the end of program 
OFSTG. It may also be initiated by requesting option 3 - 
CURVE FIT TO OFFSET DATA. 

Input request ; 

ENTER LEFT, RIGHT ELEMENTS FROM PRINTER LISTING. 

Response ; The printer output of offset values should show a 
fairly continuous variation in the mid-portion of the element 
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range with a noticeable discontinuity near each end. Enter 
the elemeitt number of the points at each end of the continuous 
midsection -of the data. For example, for day 74195 164223Z 
data these values were 180 and 2270. 

Default ; On default entry the system returns to the ATS-6 
PROCESSING menu. 

Program output ; The program fits a set of Chebyshev polynomials 
to the data points and store the coefficients in the global 
COMMON/NAVCOM/. 

4. LATSF - Read image segments from ATS-6 experimenter history 
tapes . 

Initiate by requesting option 4 - READ IMAGE SEGMENT FROM 
E.H.T. The programs A6INT, OFSTG, and OFSTF should have 
been run to set up COMMON/NAVCOM/ entries. Mount the E.H.T. 
on tape drive M^^^ or MMl. 

Input request ; 

ENTER PICTURE TIME (HHMMSS) 

Response t Enter the picture start time in packed hours, 
minutes, seconds format. This must be exactly the same 
as one of the three entries made while Initializing NAVCOM. 

(See A6INT above.) 

Input request ; 

ENTER START LINE AND ELEMENT 

Response ; Enter the image coordinates for the upper left 
hand comer of the desired image. 

Input request ; 

ENTER ZOOM AREA (1-7) 

Response ; Select a number from 1 to 7 to designate this region 
of the earth. Use the same niuriber for a given area for all 
three picture times. 
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Input re.quest ; 

ENTER TAPE UNIT (0,1) 

Response ; Enter 0 or 1 according to which drive the tape is 
mounted on (NM0 or MMl) . 

Default ; If a default entry is made, the system returns to 
the ATS-6 PROCESSING menu. 

* 

Program output : This program loads 512 x 512 element image 

segments onto digital disk files. Both visible and infrared 
data are loaded for' each request. 

5. EDGFT “ Fit a set of Chebyshev polynomials to earth edge 
measurements stored in the wind file. 

Initiate this program by requesting option 5 - CURVE FIT 
TO EARTH EDGE DATA. Earth edge measurements should be 
stored in the wind file. There are no requests for input. 

Program output ; This program fits a set of Chebyshev 
polynomials to left and right earth edge data. The 
coefficients are stored in COMMON /NAVCOM/. 

6. ATSNV - Compute a nominal attitude from landmarks. 

Initiate this program by requesting option 6 — RUN NAVIGATION. 
Landmark measurements from ATS-6 should be stored on the 
landmark file. (Currently this program reads test landmarks 
from the file DB0: [350,62]A6LMKS.DAT. , however, the program 
should be modified before general use.) There is no request 
for keyed in data. 

Program output ; This program computes yaw, pitch and roll 
values and stores them in global COMMON/NAVCOM/. In addition, 
these values are displayed on the operator terminal. The 
landmarks and computed residuals are output to the line printer. 
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B. McIDAS Version of November 1977 

1. Set up and initializing COMMON/NAVCOM/ . 

The current ATS-6 navigation system uses the files OFFSTD 
and ATSGOM. These two twenty sector files should be created 
before attempting to use the McIDAS ATS-6 navigation. 

The navigation program should be run' to initialize COMMON/NAVCOM/ 
and to set the three picture start, times used for earth edge 
corrections. To zero NAVCOM enter; 

AN 0 0 0 0 NEW 

To set picture start times enter; 

AN 0 0 0 0 HHMMSS^ HHMMSS^ HHMMSS^ 

Orbit, frame geometry and scan perigd are to be stored in the 
, McIDAS navigation file using standard McIDAS commands DQ, DS 
and ON. Entries are: 

DS SSYYDDD scan period (Psec) 

Ce.g. DS 1474195 1200000) 

DQ FIRST SSYYDDD HHMMSS Y^ (decameters) 

DQ SECOND SSYYDDD HHMMSS Y^ Z^ (decameters) 

(e.g. DQ FIRST 1474195 175531 -1198560 4041970 -43760) 

DQ SECOND 1474195 164223 133380 4214050 -5950) 

ON SSYYDDD Line-angle Total— lines Element-angle Total-elements 

Angle are in +DDDM1:KS format 

(e.g. ON 1474195 195512 2400 200412 2400) 

2. Generate alternate scan offset data points. 

Set-up ; Mount an ATS-6 experimenter history tape on a tape 
drive . Enter ; MT 14 /. 0 . 

Running program ; Enter the two letter keyin (currently GC) . 

There are no parameters. 



III-7 


Program output ; This program will print a table of element 
numbers, offsets and weights on the line printer. These same 
values are stored in the file OFFSTD. 

3. Fit a set of Chebyshev polynomials to the data points stored 
in file OFFSTD. 

Set-up ; The program to generate the^ data points must have 
been run first (see 2 above) . 

Rtmning program ; Enter the two letter keyin (currently GC) . 

The parameters are the left and right ends of the valid element 
range. These values are determined from the printer listing 
of offset values generated in 2 above. 

Program output ; This program does a least squares fit of a 
set of Chebyshev polynomials to the alternate scan data. The 
j coefficients are stored in COMtlON/NAVGOM/ . 

4. Read image segments from ATS-6 experimenter history tapes 
into digital areas. 

Set-up ; Mbunt tape. COMMOK/NAVCOM/ should have been set 
up by previous programs. 

Running program ; Enter the two letter keyin and parameters : 

BK SSYYDDD HHMMSS Area Line Element 

Program output ; This program loads a visible, infrared or 
combined area onto digital disk areas. Only standard size 
areas (500 x 672) may be used. 

5. Fit a set of Chebyshev poljmoraials to earth edge measurements 
stored in the wind file. 

Set-up ; Use WINDCO with image coordinates (WC I) and line lag 
size zero (LS 0 X) to measure earth edge displacement of both 
left and right earth edges between times t^ - t^ and t^ - t^. 

Running program ; Enter the two letter keyin (currently GC) 
to initiate the program. There are no parameters. 
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Program output ; This program computes coefficients for a 
set of Chebyshey polynomials for left and right earth edges. - 
These coefficients are stored in COMMON/NAVCOM/ . 

6. Compute- attitude from landmarks (navigate). 

Set-up ; Several landmark measurements should have been made 
from the ATS-6 images and stored in the regular IfclDAS 
navigation file. The t^ landmarks, from at least three 
different locations, should have computation code 0. Landmarks 
for later times should use code 300 and be used as a check on 
the navigation. 

Running program ; Enter the keyin and parameters : 

AN SSYYDDD 0 0 (P) 

‘ Program output ; This program computes a satellite attitude 
and stores the yaw, pitch and roll values in COMMON/NAVCOM/ 
and, via an SQ call to DX, in the McIDAS navigation file. 
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IV. Software Internal Description 

This section contains descriptions of the computer programs and 
subroutines developed for the ATS-6 navigation model. These programs 
are available on NASA’s AOIPS, on SSEC’s McIDAS and most are also 
available on the University of Wisconsin ',s Univac 1110. There are 
some variations in the main programs to allow for peculiarities of 
each system. The code for the tape read subroutines and for subroutines 
CEKTHR and CRKATS are unique to each system but yield identical results. 
Subroutines APCH, APFS and CNPS are from IBM’s Scientific Subroutine 
Package (SSP) and are not documented here. Appendix E^of this report 
contains source code listings for the McIDAS and AOIPS versions of the 
mailt programs and for most subroutines. Fig. IV. 1 illustrates the 
coordinate system used in thes^ programs. This section contains three 
parts : 

A. Description of procedure used by main programs 

B. Entries in common block NAVCOM 


C. Subroutine function descriptions 
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1 d-2400 



Figure IV. 1. ATS-6 Image Coordinate Systems. The ATS-6 satellite 
scans from south to north with 1200 scans per image. 

In a visible image each scan consists of 2 lines for 
a total of 2400 lines per image. The satellite coordinate 
lines used are numbered 1 to 2400, north to south to be 
consistent with the convention for SMS images. Satellite 
coordinate elements are numbered 1 to 2400 left to right. 
Infrared image data only have one line per scan (but 
2400 elements per line) . Infrared lines are repeated 
on image sector loads by LATSF to keep the aspect ratio 
1:1 and to keep the coordinates the same for visible 
and infrared. 
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A. This section gives a description of the procedure executed by 
each of the main programs. 

A6INT - Initialize NAVCOM 

1. Make selected entries in GOMMON/MAVCOM/. 

2. Print all entries in C OMMON/ NAVCOM/ . 

3. If two satellite orbit positions entered, call GASORB 
to convert to position and velocity.'" 

4. Note: This program only used on AOIPS. 

) 

OFSTG Generate offset data 

1. Advance to first even numbered scan to be used (scan 800) by: 

a) Read tape record (lOTPIN) 

b) Crack out scan number (CRKTHR) 

t c) If desired. scan missed print message and modify, 
start scan number 
d) Loop back to a. 

2. 'Check to see that following record contains next lower numbered 
odd scan. If not, modify desired even scan number and return 
to 1. * 

3. Desired even and odd scans' found. Back up and read whole 
records into arrays and crack out IR brightness values (CRKATS) . 

4. Compute offsets for this scan pair. See report of 31 October 
1977 Appendix G for details. 

5. Add compute'd offsets and weights to accumulated values for 
each element position. 

6. Print table of element numbers, offsets, and weights. 

7. Store number of points and element numbers in common block 
BUFFI. Store weights and offsets in common block BUFFER. 

On MclDAS these values are stored in the file OFFSTD. 

QFSTF - Fit polynomial for offset data 

1. Pick up valid element range as input. 

2. Transfer selected range of element numbers, with corresponding 
offsets and weights to array" DATI. 

3. Call subroutine APCH to set up matrix for least squares fit. 
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4* Call subroutine API'S to invert matrix and. compute coefficients 
for least squares set of Chebyshev polynomials, 

5, Store coefficients, scaling factors, number of coefficients 
and valid element range in common block NAVCOM, 

LATSF - Load ATS- 6 image segment 

1, Input information on data request (coordinates, etc,) 

2, Set up image label and write to disk, 

3, Call subroutine GENOFF to set up a table of offsets, for 
every element position to be read in, by evaluating the set 
of Chebyshev polynomials stored in common block NAVCOM, 

4, Advance to first data record, 

a) Skip two header records 

b) Read a partial record from tape 

: c) Call CRKTHR to crack out scan number 
d) If desired start scan not reached, go back to b. 

^5, Back up one record so entire record can be read. 

6. Read image segment 

a) Read a record* 

b) Check scan number. If less than last scan to be read in, 
terminate image load. 

c) Pass Visible-2 data to subroutine LINGRB to select and 
repack desired line segment. 

d) Write image line segment to digital disk area. 

e) Pass Visible-1 data to subroutine LINGRB to select and 
repack desired line segment. 

f) Write image line segment to digital disk. 

g) Pass Infrared data to subroutine LINGRB to select and 
repack desired line segment. 

h) Write infrared image line to digital disk. 

i) Write infrared image line to digital disk a second time 
so that visible and infrared coordinates .match. 

Done now; rewind tape. 


7 . 



IV-4 


EDGFT - Fit polynomials to earth edges 

1. Do curve fit for t^ — t^ image pair then for t^ — t^ image pair. 

a) Bead a wind from disk file. 

b) Check for valid year, day, times, error code. If invalid, 
return to step a. 

c) If element position of vector less than picture center 
element, store scan number, shift in array for left edge. 

d) If element position of vector greater than picture center' 
element, store scan number, shift in array for right edge. 

e) Loop back to a till end of wind file encountered. 

f) Fill array DATI with scan numbers, shifts for left edge. 

Pass array to subroutine APCH to set up matrix. 

g) Pass matrix from APCH to sub routine. APFS to invert matrix 
and compute coefficients for set of Chebyshev pol 3 niomials . 

h) Store coefficients, scaling factors in common block NAVCOM. 

i) Repeat steps f, g, h for data from right edge. 

j) Determine valid argument range for left and right edges. 
Determine the overlapping portion of valid argument range 
for left and right edges and store this overlapping portion 
in common block NAVCOM. 

ATSNV - Navigate ATS-6 image from landmarks 

1. Read landmarks from landmark file. 

2. Convert integer values to floating point. 

3.. Pass landmark data to subroutine ATTTUD to compute satellite 

attitude. 

4. ' The lIcIDAS version stores the attitude in the navigation file 

by SQing DX. 

5. Compute and list residuals. 

a) Pass picture time, latitude, longitude of landmark to 
subroutine ES to compute line and element. 

b) Compute residual equals measured value minus computed value 
for lines and for elements. 

c) List values and loop back to a through all landmarks. 
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B. Entries 

NATO 

INAV 

lYR 

IDAY 

TOTLIN 

DEGLIN 

TOTIEL 

DEGELE 

PICLIN 

PICELE- 

TMRSCL 

lOYR 

lODAY 

TM 

RiX 

RIY 

KLZ 

RIDX 

KLDY 

RIDZ 

PITCH 

ROLL, 

YAW 


in common block NAVCOM. 

Navigation sequence number. 

Flag to indicate type of navigation. 

Year of date for which navigation is valid. 

Day of year for which navigation is valid. 

Total number of lines in image (= 2400 . ) 

Total sweep angle in line direction (= 19.92 degrees). 
Total number of elements across image (= 2400.). 

Total scan angle in element direction (= 20.07 degrees). 
Picture center line. 

Picture center element. 

Scan period (nom. .02 minutes) . 

j 

Year of date for orbit values (lOYR = lYR) . 

Day of year for orbit values (lODAY - IDAY) • 

Time of orbit location (minutes , GMT). 

X“Component of location of time TM (earth radii) . 

Y component of satellite location of time TM (earth radii) 

Z component of satellite location of time TM (earth radii) 

X component of satellite location of time TM (earth 

radii /minute ) . 

Y component of satellite velocity at time TM (earth 
radii/minute) ♦ 

Z component of satellite velocity at time TM (earth 
radii /minute) . 

Pitch angle pf rotation from BC to PF coordinates. 
Satellite attitude (radians) . 

Roll angle of rotation from BC to PF coordinates. 

Satellite attitude (radians) . 

Yaw angle of rotation from BC to PF coordinates. 

Satellite attitude (radians) . 
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PTIM(3) Picture start titles for the three images for which 

earthedge correction applies. PTIM(l) is the time 
of the reference image. (Time in minutes, GMT) 

Minimum time of time over which earthedge correction 
valid. Subscript corresponds to image number as in 
PTIM. 

TMX(3) Maximum time of time over which earthedge correction 

valid. (e.g. earth edge correction may be used for 
scan -time t within image starting at 'PTIM(2) only if 
TMN(2) _< t £TMX(2). Units are minutes, GMT. PTIM(i) < 
TMN(i) _< TI«(1) . “ 


No^: For earthedge correction arrays, dimensions of value 2 refer 

to image pair. 1 => PTIM(l) - PTIM(2) shifts; 2 => PTIM(l) 

- PTIM(3) shifts. 


NLC0EF(2) 
NRC0EF(2) 
SCLL0 (2) 

SCLL1(2) 

ELC0EF(11,2) 

SCLR0 (2) 

SCLEl (2) 

ERC0EF(11,2) 

NASCEF 

SCLAS0 

SCLASl 

lELEMK 


Nutriber of coefficients in polynomial for left edges. 

Number of coefficients in polynomial for right edge. 

Offset for scaling argument value (scan number) for left 
edges . 

Multiplier for scaling argument value (scan number) 
for left edges. 

Coefficients of set of Chebyshev polynomial for left 
edges. 

Offset for scaling argument value (scan number) for 
right edges. 

Multiplier for scaling argument values (scan number) 
for right edges. 

Coefficients of set of Chebyshev pol 3 momial for right 
edge. 

Number of coefficients in polynomial for alternate 
scan correction. 

Offset for scaling argument (scan number) for alternate 
scan offset. 

Multiplier for scaling argument value (scan number) for 
alternate scan correction. 

Minimum element number for which alternate scan correction 
applies . 


ORIGINAL PAGSsISc 
0? P(X)R QUALITY / 
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lELEMX 

ASC0EF(16) 


Maximum element number for which alternate scan 
correction applies. 

Coefficients fo set of Chebyshev polynomials for 
alternate scan correction. 
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C. Subroutine Function Description" 
Name ; BCTOPF 

Call ; CALL BCTOPF (X, Y, Z, IDIR) 


Input Parameters : 

X, Y, Z Satellite orientation in body centered (IDIR = 1) or 
picture frame (IDIR = 2) coordinates 

IDIR Direction to rotation 

IDIR = 1 Body centered to picture frame 
IDIR = 2 Picture frame to body centered 


Re turned Values : 

X, Y, Z Satellite orientation in picture frame (IDIR = 1) or 
body centered (IDIR = 2) coordinates . 


Algorithm : 

BCTOPF uses subroutine ROTATE to create a rotation matrix from the yaw, 
roll and pitch angles computed by the navigation program. It then 
multiplies the vector (X, Y, Z) by the matrix (IDIR = 1) or its transpose 
(IDIR = 2) . 


Reference : 

For discussion of the coordinate systems see "Design and Testing of the 
Navigation Model for Three Axis Stabilized Earth Oriented Satellites Applied 
to the ATS-6 Satellite Image Data Base" progress report for 17 Nov. 1975, 
appendix pp* 2-4 or progress report for 31 June 1976, appendix pp. 2-4. 
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Name: 

CRKA-TS 

Call: 

CALL CRKATS (N 


Input Values ; 

N - Number of data words to crack out, 

S - Source array of 12-bit words stored in sequence as a continuous 
bit string after being read from tape. 

Returned Values : 

P - Destination array of full words* 


Function ; 

ATS-6 experimenter history tapes are written with the image data stored in 
the 9 least significant bits of successive 12-bit words. CRKATS extracts 
the 8 most significant data bits (bits 8-1 of a 12-bit word numbered 11-0) 
and stores the resulting data value in a full computer word. 
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Name ; CEKTHR 

Call ; CALL CRKTHR (N, S, D) 

Input Values ; 

N‘ - Number of 12-bit words to crack out. 

S — So.urce array of 12— bit words stored in sequence as a continuous 
bit string after being read from tape. 

Returned Values ; 

D — Destination array of whole words. 


Function: 


ATS-6 experimenter history tapes are written with the data stored in 
12.-bit words (actually thirds of 36-bit words) . CRKTHR places these 
12-bit words in whole words (16-bits for PDP-11) . 
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Name : EDGCOR 

Call : CALL EDGCOR (PTIME, ALIN, DELLIN, DELELE) 

Input Parameters : 

PTIME - Picture start time (minutes, GMT)* 

ALIN - Line coordinate value at which correction is to apply. 

Returned Values : 

DELLIN - Line correction value. 

DELELE - Element correction value. 


Function : 

EDGCOR requires that left and right earth edge shift polynomials be stored 
in COMMON/NAVCOM/. The current version only allows for three picture start 
times to be in use at once. 

EDGCOR evaluates the earth edge shift polynomials for left and right edges, 
then converts these values to line and element shifts. The polynomial 
only applies to the line range for which earth edges were measured. Outside 
this range, the value zero will be returned. 
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Name ; ERTOER 

Call ; CALL ERTOER (XLAT, XLON, SE, YE, ZE, IDIR) 


Parameters : 


XLAT 

XLON 

XE, YE, ZE 
IDIR 


Latitude (degrees) of a point on the earth. 

Longitude (degrees) of same point on the earth. 

Cartesian coordinates in, coordinate system rotating with 
earth (kilometers) . 

Direction of transformation. 


Function ; 

Converts coordinates of a point on the earth from latitude, longitude 
to rotating Cartesian coordinates (IDIR = 1) or vice versa (IDIR = 2) . 


References ; 

See progress report of 17 Nov. 1975 or 31 June 1976, appendix pp. 2-4 
for coordinates. 


See also subroutine ERTOST. 



IV-13 


Name ; ERTOST 

Call : CALL ERTOST (XE, YE, ZE, X, Y, Z, IDIR, TIME) 

Parameters ; 

XE, YE, ZE 

X, Y, Z 

IDIR 

TIME 

Function. ; 

ERTOST a unit vector pointing from the satellite to a given point on the 
earth (IDIR = 1) or given a pointing vector it computes the location, if 
any, on the earth's surface that the vector is pointing at (IDIR = 2). 

Reference : 

See progress report of 17 Nov. 1975 or 31 June 1976, appendix pp. 2-4 
. for coordinate. 

Related subroutine ERTOER, STTOLV. 


Cartesian coordinates (kilometers) of a point on earth's 
surface in rotating coordinate system. 

Unit vector in inertial coordinates pointing from satellite 
to point (XE, YE, ZE) on Earth's surface. 

Direction of transformation 

IDIR = 1 => (XE, YE, ZE) ->• (X,‘ Y, Z) 

IDIR = 2 => (X, Y, Z) (XE, YE, ZE) 

Time of day (minutes, GMT) 
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Name : ES 

Call ; CALL ES (PTIME, ^AT, XLON, XLIN, XELE) 


Input Parameters : 

PTIME - Picture start time (minutes, GMT). 

XLAT - Latitude (degrees) of a point on the earth^s surface. 

XLON - Longitude (degrees) of that point. 

Returned Values ; 

XLIN - Line coordinate in the image picture frame coordinate system. 
XELE - Element coordinate in image picture frame coordinate system. 


Function : 

Subroutine ES does an earth (latitude, longitude) to satellite (line, 
element) coordinate transform based on the satellite’s attitude, orbit 
position, and, if available, attitude correction based on measurement 
of earthedge shifts. 


References : 


See subroutine SE. 
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Name ; FLALO 

Call ; XLAT = FLALO (ILAT) 

Input Value ; 

ILAT = Latitude {or longitude) integer value in the format ^WDDtlMSS. 
(For PDP-11 an INTEGER*4 value.) 

Return Value ; 

XLAT = Latitude (or longitude) in degrees (floating point). 

Function ; 

Converts an angle stored as an integer in degrees, minutes, seconds 
format to floating point degrees. 



Name : FLIP 


Call : CALL FLIP (A, B, I, N. ALTRET) 

Parameters : 

A - an NxN matrix 
B - an NxN matrix 

I - row on which to perform operation 
N - dimension of A and B 

ALTRET - flag indicating an error (LOGICAL) 

Return : 

Aj B are returned in modified form* 

Function : 

Ail rows greater-than I are added to row I. The same operation is 
performed on matrices A and B. 

Reference : 

This subroutine is used only by subroutine INVERT. 
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Name ; FTIME 

Call ; TIME = FTIME (ITIME) 

Input Parameters ; 

ITIME = Integer time of day in the form HHMMSS. (INTEGER*4 on the PDP-11,) 
Returned Value ; 

TIME = Time of day in minutes (floating point) • 
j^onction : 

Converts a time of day in the packed integer format hours, minutes, 
seconds to time of day in minutes (floating point). 
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Name ; GASORB 

Call ; CALL GASORB (Rl, Tl, R2, T2) 

Input Parameters' : 

Rl Position vector of satellite at time Tl in earth inertial 
reference frame. 

!? 

Tl Time (minutes, GMT) at which satellite is a position Rl. 

R2 Position vector of satellite at time T2. 

T2 Time (minutes, GMT) of position R2. 

Returned Values ; Note - results are stored in COMMON/NAVCOM/ . 

TM - Time of position. TM = Tl 

RIX, RIY, RIZ- Position of satellite at TM. 

(RIX, RIY, RIZ) = Rl/RE 
Where RE = radius, of earth. 

* 

RlDX, RIDY, RlDZ - Velocity of satellite at time TM. 

Ftmctlon : 

Given two position vectors and their corresponding times, GASORB computes 
the position- and velocity of the satellite at the time of the first given 
position vector. The method used in an f,g series from the method of Gauss. 

Reference ; 

Escobal, P. R. Methods of Orbit Determination , John Wiley & Sons, 1965, 
pp. 196, 197. 

See also subroutine ORBIT. 
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Name: GENOFF 


> 

Call: CALL GENOFF (lUELE, NELES) 


Input Parameters : 

lUELE First element of offset array. 

NELES Niimber of elements in offset array. 


Fmetion ; 

The subroutine GENOFF evaluates the alternate scan correction polynomial 
for all element values across the image to be read in. The values are 
stored in an array in COMMON/OFFSET/. This array is then used by the 
program which loads ATS-6 images (LATSF or LDATSF) . 


References ; 

See information on programs LATSF and OFSTF and on COMMON/NAVCOM/ . 
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Hame ; INVERT 

Call ; CALL INVERT (AA, B, N, ALTRET) 

Input Parameters ; 

AA - an NxN matrix 
N - dimension of AA and B 

Return Values : 

B — the inverse of AA 

ALTRET - a flag to indicate AA is singular (LOGICAL) 

Function ; 

INVERT returns in B the NxN inverse of the matrix AA. If AA is singular, 
ALTRET is set to .TRUE. 

Reference ; 


See also subroutines FLIP, MINMIZ. 
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Name : LINGKB 

Call : CALL LINGRB (INDATA, ISDIR, lELE, NELES, IBDF, lOUTD) 


Input Parameters : 

INDATA An array containing the bit string for either the Visible-1, 

Visible-2 or Infrared sensor as read from tape. Array actually 
starts 72 bits before first data word. 

ISDIR Indicates scan direction. 0 => Even numbered scan. 

1 => Odd numbered scan. 


lELE First element of desired line segment. 

NELES Number of elements in line segment. 

IBDF Sampling factor. (IBDF = 1 only) 


Returned Values : 

lOUTD Output array with data .packed one pixel per 8-bit byte. 


Function ; 

This subroutine unpacks pixel data from the bit string read from magnetic 
tape. It then selects out the desired line segment (512 elements on AOIPS) , 
shifts the even scans based on the evaluation of .the alternate scan correction 
polynomial, and stores the pixels in an array to be written in the image file. 


References : ' • 


See main program LATSF. 
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Name: LS 

Call : CALL LS(X, Y, VAL, DD, DIR) 

Input Parameters t 

X Starting point for line search (a vector) . 

DD Unnormalized directional derivative VS(Y) • DIR. 

DIR Direction to do search (a vector) . 

Returned Values ; 

Y The selected point (a vector) . 

VAL Value of the objective function S evaluated at Y. 


Function ; 

This routine performs an Armijo line search from the point "X" in the 
direction "DIR" and returns the selected point in "Y" and the objective 
fimction v alye S (Y) in "VAL". On call, "DD" is the imnormalized directional 
derivative VS • DIR. This line search routine returns in Y the point 
X + 2-N*DIR where N is the least nonnegative integer such that -S(2 ^*DIR) 
represents at least 40% of the functional drop in the linearization of S 
at X in inoving from X to X + 

References : 

See also subrountine MINMIZ and function S. See the report for 17 Nov. 1975 
Appendix section IV or report for 31 June 1976 Appendix section II. C. 
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Name ; MINMIZ 

Call ; CALL MINMIZ (PTIN, PTODT, GNOEM, VAL, ITN) 

Input Parameter ; 

PTIN Starting point for search for minimum value of objective function 
(PTIN is a vector) 

Returned Values ; 

PTOUT .Optimal point. The objective function has a minimum at PTODT. 
(PTOUT is a vector) 

GNOEM The norm of the gradient of the objective function at PTOUT. 

VAL The value of the objective function at PTOUT. 

ITN Number of iterations done. 


Function : 

MINMIZ finds the point PTOUT at which an objective fxmction is a minimum. 
In this case PTOUT is the satellite's attitude (PITCH, ROLL, YAW). The 
objective function is defined in the reports (report of 17 Nov. 1975 
Appendix eqn. 15, report of 31 June 1976 Appendix eqn. 17). Basically 
MINMIZ serves as a driver for PRTIAL, INVERT and LS. 


References ; See also subroutine ATTTUD. 
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Name ; NKMLIZ 

Call ; CALL NRMLIZ (VX, VY, VZ, WORM) 

Input Parameters ; 

VX, VY, VZ Cartesian components of any vector. 

Returned Values : 

VX, VY, VZ - Normalized components of the input vector ♦ 

VNORM “ Length of the input vector - 

Function ; 

NRMLIZ computes the length of the vector with components (VX, VY, VZ) 
then divides each component by that length to return a unit vector. 
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Name ; ORBIT 

Call ; CALL ORBIT (X, Y, Z, T) 

Input Parameters ; 

T Time of day (minutes, GMT) 


Returned Values ; 

X, Y, Z - Position of satellite at time T (kilometers) 


Function ; 

Given the position and velocity of the satellite at some reference time 
as computed by GASORB, the subroutine ORBIT computes the satellite's 
position at the time T. 


Reference ; 

Escobal, P.R. Methods of Orbit Determination , John Wiley & Sons, 1965, 
pp. 427, 428. 
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Name : PFTOTC 

Call : CALL PFTOTC (XLIN, XELE, X, Y, Z, IDIR, INIT) 

Parameters ; 

XLIN Line number of a point on the ATS-6 image. 

XELE The element number of that point. 

X, Y, Z Unit vector, in the picture frame coordinate system, pointing 
at location defined by (XLIN, XELE) . 

IDIR Direction of coordinate transformation. 

INIT Initialization flag. Set INIT = 1 before first call. 

Function : 

PFTOTC converts from picture frame cartesian (X, Y, Z) coordinates to 
picture frame image (LINE, ELEMENT) coordinates (IDIR = 1) or vice versa 
(IDIR = 2) . 

References : 

See progress report of 17 Nov. 1975 or 31 June 1976, appendix, pp. 2-4, 
for coordinates. 

See also subroutine BCTOPF. 
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Name : PRTIAL 

Call: CALL PRTIAL (PT, GRAD, HESS) 


Input Parameters : 

PT Point (vector) at which to evaluate GRAD and HESS. 

Returned Values : 

GRAD The gradient of the objective function, evaluated at PT. 
HESS The hessian of the gradient function evaluated at PT. 

Function : 

PRTIAL .computes values of gradient and hessian for a given function. 


References ; 


HINMIZ. See the reports of 17 November 1975 Appendix 
page 10 ff or report of 31 June 1976 Appendix page 13 ff. 
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Mame ; ROTATE 

Call ; CALL ROTATE (A, R, IR, IDERIV) 

Parameters ; 

A A matrix. 

R An angle of rotation (radians) 

IR Axis number (1, 2, 3) • 

IDERIV Derivative of rotation matrix. 


Function : 


This routine returns in "A” the product of the input matrix "A" and a 
matrix EM, where, if "IDERIV=1", EM represents a rotation through an 
angle "R" (in radians) about the axis "IR". IF IDERIV=2, the first 
derivative of EM is operated on A, and if IDEEIV=3, the second derivative 
of EM is used. 


References ; See also subroutine PRTIAL. 



IV-29 


Name; S 

Call ; SVAL = S (PT) 

Input Parameter ; 

PT The point (a vector) at which S is evaluated. 

Returned Value ; 

S The objective function which is minimized in computing the 

satellite attitude. 


Function ; 

S is the objective function which is mini m ized by MINMIZ in computing 
the satellites attitude. See the report Of 31 June 1976 Appendix 
equation 17 or report of 17 Noveiriber 1975 Appendix equation 15. 


References : 


See also the subroutines MINMIZ and LS. 



Name : SATEAR 


Call ; CALL 

Parameters ; 

PICTIM- 

XLIH 

XELE 

XLAT 

XLON 

ITYPE 

INAV 

BETAIN 

BETDOT 

ATPRAC 

Function: 


SATEAR (PICTIM, XLIN, XELE, XLAT, XLON, ITYPE, INAV, RETAIN, 
BETDOT, ATFRAC) 


Picture start time (minutes, GMT) 

Satellite image line (master coordinate) 

Satellite image element (master coordinate) 

Latitude (degrees, 4-North, -South) 

Longitude (degrees, 4East, -West) 

Type of conversion 

1 => Satellite to earth 

2 => Earth to satellite 

Dunffiiy variable used so call will match SMS version. 


. Subroutine SATEAR calls subroutines SE or ES. 
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Name : SE 

Call: CALL SE (PTIME, XLIN, XELE, ALOT, ALON) 


Input Parameters : 

PTIME - Picture start time (minutes, GMT). 

XLIN - Line number of a point on the ATS-6 image. 
XELE - Element number of that point. 


Returned Values : 

ALAT - Latitude (degrees) of the point on the earth* s surface at (XLIN, XELE), 
ALON - Longitude (degrees) of that point. 


Function : 

Subroutine SE does a satellite (line, element) to earth (latitude, longitude) 
coordinate transform based on the satellite’s attitude, orbit position, and, 
if available, attitude corrections based on measurement of earthedge shifts. 


References : 


See also subroutine ES. 
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Name: STTOLV 

Call : CALL STTOLV ( X, Y, Z, IDIR, Tip) 


Parameters : 


X, Y, Z 
IDIR 


TIME 


Unit pointing vector in the satellite inertial coordinate 
system or in the satellite local vertical coordinate 
system. • 

Direction of transformation. 

Satellite inertial to local vertical (IDIR =1) 

Local vertical to satellite vertical (IDIR = 2) 

Time at which transform applies. 

(minutes , GMT) 


Function : 

STTOLV uses the satellite's position at time TIME, and converts a unit 
vector in the satellite Inertial coordinate system to the satellite's 
local vertical system (IDIR = 1) or vice versa (IDIR = 2) . 


Reference : 

See progress report for 17 Nov. 1975 or 31 June 1976, appendix, pp. 2-4. 
See also subroutines ERTOST, LVTOBC, and BCTOPF. 
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Name ; UNIT 

Call ; CALL UNIT (A) 

Returned Values ; 

A - a 3x3 identity matrix 

Function ; 

Subroutine UNIT returns a 3x3 identity matrix in array A. 



References 


1, Kiihlow, W. W. , Design and Testing of the Navigation Model for Three 
Axis Stabilized Earth Oriented Satellites Applied to the ATS-6 
Satellite Image Data Base , Progress Report for Period Ending 

C» 

17 November 1975, Contract Number NAS5-20974, University of Wisconsin. 

2, Kuhlow, W. W. , G. C. Chatters^ Design and Testing of the Navigation 
Model for Three Axis Stabilized Earth Oriented Satellites Applied 
to the ATS- 6 Satellite Image Data Base , Progress Report for Period 
Ending 31 June 1976, Contract Number NAS5^20974, University of 
Wisconsin, 

3, Kuhlow, W, W. , G, C. Chatters, Design and Testing of the Navigation 
Model for Three Axis Stabilized Earth Oriented Satellites Applied 
to the ATS-6 Satellite Image Data Base , Progress Report for Period 
Ending 31 October 1977, Contract Number NAS5-20974, University of 


Wisconsin, 
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NOMENCLATURE 


Coordinate Systems (C.S.) 

El = Earth- centered Inertial 

ER = Earth- centered Rotating 

LV = Local Vertical 

(satellite-centered) 

1 

BC = Body Centered 

(satellite-centered) 

PF = Picture Frame 

(satellite-centered) 

Orthogonal Matrices 

= rotation into LV C.S. 

= rotation into BC C.S. 

Rpj, = rotation in PE C.S. 

R^ = optimized "navigation" 
matrix 

R(e,k) = rotation about axis 
k (k-1,2,3) in a ccw 
sense by an angle 6 

Vectors 

y\. 

r = unit vector 

satellite radius vector 
= landmark pointing vector 
^2 “ earth-coordinate, El C.S. 

= earth-coordinate, ER C.S. 
’^LV ~ Poiiiting vector in LV C.S. 
^BC ^ pointing vector in BC C.S. 
^PF pointing vector in PF C.S. 

63 - (Q.C,!)"^ 


Other 

®y»6^,6 = attitude yaw, roll, 

pitch angles 

- geodetic latitude, longitude 

“ sidereal time of Greenwich prime 
meridian 

~ nadlans/line, radians/element 

X = mirror step angle 

6 = mirror sweep angle 

~ navigation parameter 

— optimized navigation parameter 

LIN, ELE = satellite image coordinates: 
Line, Element 

= picture-center coordinates 
(line, element) 

e = eccentricity of earth oblate 
spheriod model 

n „ = earth's equatorial radius 
Tp = earth’s polar radius 
® “ distance from satellite to landmark 
t = time 

= epoch time, lies within image 
frame interval 

• IJ = Euclidean norm- 

•| = absolute value operator 
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I. PRELIMINARIES 
A. Coordinate Systems 

All coordinate systems used in the model are 3-D right-handed orthogonal 
coordinate systems. There are five all together. Two have their origins 
placed at the dynamical center of the earth. 

The plane formed by the x- and y-axes of the earth-centered inertial 
coordinate system (El) lies in the earth’s equatorial plane. The x-axis 
points at the vernal equinox (y) which is assumed to be inertially fixed 
and the z-axis points north. Rotating relative to this inertial frame is 
the earth-centered rotating coordinate system (ER) with its x-axis passing 
through the Greenwich meridian and its z-axis coincident with the El z-axis 
(Fig.> A.l). 

In the local vertical (LV) system, the z-axis points to the center of 
the earth, i.e. the unit vector representing this axis at time t is given 
as z = -r^ (t)/[|r^ (t)jj , where r^ is the satellite radius vector (Fig. A.l). 

The x-axis is parallel to the earth’s equatorial plane and nominally points 
east. 

Fixed in the satellite body is the body-centered (BC) coordinate 
system whose axes are nominally coincident with the LV system. Departures 
from this alignment are measured by the yaw, pitch, and roll time dependent 
angles which in part make up the attitude telemetry data. These rotations 
are explicitly defined later on. Also nominally coincident with these 
coordinates systems is the last to be defined, the image or picture frame 
(PF) coordinate system (Fig. A. 2). The z-axis -points to the picture center 
(earth image center) which for ATS-6, is the image point occurring at the 

V 

midpoint of the mirror sweep angle for the mid-mirror' scan number of a full 
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image frame scan* (For ATS— 6, one coinpXete image scan consists of 1200 
mirror scans with 2400 samples or elements per scan line*) The PF x*-axis 
is parallel to the center mirror sweep scan at the picture center and 
nominally points east.^ 


Convention for Orthogonal Matrices 

' Two basic orthogonal transformation representations in transforming 
a vector from one coordinate system to another with common origin are given 
here. Their forms arise naturally in the development of the navigation 
model depending on convenience or available information. 

^1^^^ two coordinate systems with a common origin whose axes are 


unprimed xyz and primed x*y^z* respectively, let R represent the orthogonal 
transformation of the vector r, (whose components are expressed in the 
unprimed system) to the vector whose components are expressed in terms 
of the primed system, i,e. r* = Rr. 

The first form R can take is, in matrix representation, R = [xlylz]"^. 

As 

where x, y, z are the unit column vectors of the primed coordinate axes 
whose components are expressed in the unprimed coordinate system. The 
^ to the transpose of the matrix. To see that this transformation 

is valid, one need only to carry out the operation implied, 


r’ = [xly|z] r = (xt, yr, z*r) . 


Thus, x*r, y*r, z*r represent the projections of r onto the x^, y’, 

axes respectively. This is precisely the representation of r in the primed 
system. 


The second form of R is expressed in terms of rotation angles where in 
its simplest form R = R(e,k) represents a rotation by an angle 0 counter- 



A-6 


clockwise about the k-th axis as viewed from above where k=l,2,3 refers 
respectively to the x,y,z-axes. Thus, using* the conventions defined above, 
if the z and z’ axes were coincident and the x'y’z’ system were rotated by 
an angle 6 counterclockwise relative to the xyz system about the z-axis, 

cos 6 sin 0 0 

-sin 0 cos 9 0 

0 0 1 




In general for R(9 ,k), 

\i ^ \k '*■1 if X = k^ ’ \i ^ “^1i ^ * 


i-J Jt 
i < j and i,j 5 ^k 


\?ith this notation, the sequence of Euler rotations can he represented 

in a compact form. For example, the three rotations about the ATS--6 body 

axes to define the ATS--6 attitude relative to the LV coordinate system are, 

in sequence, a ccw rotation by an angle 0^ about the BC z-axis (yaw) , 

followed by a ccw rotation by an angle 0^ about the BC x-axis (roll) , 

followed by a cw rotation by an angle 0 about the BC y-axis (pitch) . Thus 

P 

a vector expressed in the LV system is transformed to the BC system 


(r ) by the operation: 
BC 


= R(-0p,2) R(e^,l) R(0y,3) , [LVTOBC] 


where the angles 0^, 0^, 0^ are generally dependent on time. 


C. Orbit Methods Used in Predicting Satellite Position, r (t) 

s 

The orbit parameters used in the model are derived from the ephemeris 

data available from the ATS-6 magnetic tapes. These data are in the form 

of position (r ) and velocity (r ) components expressed in the El coordinate 
s s 

system approximately every three seconds of ephemeris time. Por a two body 

*Note: Label in brackets (e.g. [ROTATE]) is corresponding PORTRAN subroutine. 

See Appendix E for listing. 



orbit, a position and velocity vector at any time is theoretically 
sufficient to uniquely determine the orbit • However, the velocity vectors, 
unlike the position vectors, are not given with sufficient accuracy for 
the purposes of this model. Therefore, a Gaussian orbit determination 
method is used in which two position vectors at different times are used 
to calculate an accurate velocity vector corresponding to one of the selected 
position vectors. In symbols, 

♦ 

), t ^ t . [GASORB] 
S 1 S 2 sisi 12 

The vector spread between r (t ) and r (t ) must be less than 70 Once 

. . Si s 2 

r (t ), r (t ) are determined, the satellite’s position for any other time 
is determined by an iterative f, g computational procedure: 

r (t) = f r (t ) + g r (t ) . [ORBIT] 

S S ^ * S 

The details of the computational algorithms for these two procedures are 
given elsewhere. A thorough testing of these two routines which are 
incorporated into the ATS-6 model has indicated that any errors generated 
by them or by uncertainties in the satellite’s position are negligibly 
small compared to errors arising from other sources. 

II. THE NAVIGATION MODEL 

A. Earth Coordinates to Satellite Image Coordinates (Subroutine ES) 

Let t be the instant at which a point on the earth (landmark) is imaged 

by the' satellite scanning system and let r^ (t) be the vector which points 

from the satellite-center to the landmark in question. However, t is not 

\ 

known. Therefore, we make a guess then use the derived line to get a better 
value and iterate till the solution converges to within a line. The relation 

P. R. Escobal, Methods of Orbit Determination , J. Wiley and S'ons, New York, 
1965. Gaussian orbit algorithm, pp. 196-197, f, g method, p. 423. 
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of this vector to the satellite radius vector r and the landmark r is 

s 2 


given by: 


r (t) = r (t) - r (t) , [ERTOST] 

1 2 S 


( 1 ) 


where r^, and therefore r^, are expressed in El coordinates (Figure 1). 

The position of the satellite r (t) is determined from the orbit routine 

s 

discussed above, while r is derived as follows: 

2 

- K(-e^,3) \ , 

o T 

where r^ = p(cos 6^ cos iji^ , cos 6^ sin <|)^, (1 - e^) sin 6^^) , 


[ERTOST] 


( 2 ) 




L :l 

ir /^l - e^sin^ 0^ , 
eq L 

the geodetic latitude and longitude of the landmark. 


e “ eccentricity of the oblate spheriod earth model = 8.1812X10’'^ , 

r = earth’s equatorial radius = 6378.15 km. 
eq 

The landmark r^ expressed in ER coordinates is transformed into the El 
E 

coordinate system via the transformation R(-6_,3) where, 0 , the angle 

E E ' 

between the x-axes of the two systems at time t is given by 

6 = a + a * DDD + a * t [ERTOST] 

E 1 2 3 


where, 6„ is the sidereal time of the Greenwich prime meridian, 
E 


DDD is the day of the year , 
t is universal time , 

a ,a ,a are constants derived for a specific Julian date; for January 0, 1974 
12 3 

= 99.59477026 degrees 
= .985647336 degrees/DDD 
a = .2506844773 degrees /decimal minute. 



It is convenient at this stage to transform the pointing vector r 

into a unit vector, r (t) = r (t)/|r (t)||. Transformation of r into 

1 1 ' 1 ' 1 

the LV frame is accomplished by 


[STTOLV] 

where, from the definition of the LV coordinate system, 

a . -?^(t)/|pjt)|| - 

X = (z X e )/|lz X e II = i~y/d, x /d, 0)^, 

3 3 , s s 

y = z X X = (x z /d, yz/d, -d^)^, 

S S o o 


• (3) 

\ 


d =v/x2 + y|. 


where 1. 

s -^s s 

In matrix form. 


®LV 


1 

d 



X 

0 

s 

s 


X z 

y z 


s s 

s 


-dx 

s 

CO 

•o 

1 

-•dz, 

s 


The transformation from LV to BC coordinates is accomplished by 
using the attitude telemetry data. However, since we have found these 
data to he unusable for our purposes, the LV and BC coordinates are set 
equal to each other. Thus: 


BC 


^LV 


(4) 


The transformation of r__ into the image frame vector, r_„, "is 

A T ^ 

^PF ~ ^^F’^PF’^PF^ ^F ^BC * 
where = R(?p,2) R(?^,l) R(e’y,3), 


(5a) 


(5b) 
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and 0^, 9^ are ^the optimized angles restilting from the least-squares 
navigation technique employing landmark measurements* Since is the 
transformation from the satellite's body-centered frame to its picture 
(image) frame, it is expected that this transformation will be reasonably 
time- in depan den t . The effect of thermal and mechanical stresses oti the 
stability of 0^, 6^, 0^,over a period of time on the order of days is some- 
thing that can only be deduced indirectly by updating the navigation* 

With r^^ determined from equations (5a) and (5b) and adding line and 
element corrections deduced from earth edges, the satellite image coordinates 
can be calculated 

LIN = L + (sln-1 y-ni 7 )/pT + [PFTOTC] (6a) 

c ir-b L e 

ELE = + (tan(Xp^/ZpP)/p^ + AE^ and [EDGCOR] (6b) 

where, 

LIN = line number, 

ELE = element number, 

AL = line correction due to attitude shifts, 
e ^ 

AE = element correction due to attitude shifts, 
e ' 

L = picture center line = 1200 for visible, 600 for IR ATS-6 image 
c 

data. 


E = picture center element — 1200 for ATS-6 image data, 
c 


p = number of radians /line 
L 

pg == number of radians /element 


nominal ATS-6 scanned field is 20^ X 20*". 


Summarizing equations (1) - (5), 

^PF “ ^F \c \v * 


[ES] 


where r^*||r^|j * ^^(t) - r^Ct) = R(-0g,3) * 

Equations (7a) and (7b) effectively define the earth coordinates to 


(7a) 

(7b) 


satellite linage coordinates-transformation. 
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B. Satellite Image Coordinates to Earth Coordinates (Subroutine SE) 

For a- given LIN, ELE, it is apparent from equations (6) and Figure 2 

that 

I'pj. = (“COS Xsin 5, -sin X, cos Xcos 6)^ [PFTOTC] (8a) 

where 

^ “ LIN + AL^) = mirror step angle (8L) 

(E^ - ELE + AE^) pg = mirror sweep angle (8c) 

From equation (7a) , 

h - ») 

■(The three successive transformations in equation (9) are orthogonal 

matrices; therefore the inverse of each is equal to its transpose). 

Now r = r + sr , (10) 

® 1 

where s = j|r^ (t) jj equals the distance from the satellite-center to the 
landmark. The solution of s is achieved by using eq. (10) and the equation 
of the earth spheroid, 

where , 

= earth's equatorial radius = 6378.15 km, 

r = earth's polar radius = 6356.77 km, 

P 

^E’^E’^E vector components of .the landmark, r^^, in the ER frame. 


Equations (2), (lO) and (11) represent a system of equations of four 

unknowns (x_,y„,z_,,s) • The solution of s can easily be accomplished as 
E £ ri 

follows ; 


Divide the x and y components by r and the z component by r in (10) ♦ 
^ P 
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This results in the equation. 


r*=r*+sr* 
2 s 1 


( 12 ) 


where. 


r * - (x /r _>y /r ,z /r )^, 

1 1 ®«1 1 eq 1 P 

r * = (x„/r ,y /r ,z /r . 
s s eq s eq s p 

Nowlir *|P = II R(-0x,»3)|| ^ • ||r„*|| - 1, since R(8 ,3) is an orthogonal 

2 -K ii ii 

matrix, and||rJP equals the left side of equation (11). Therefore 
£• 


the equation 


r^*|| = r=||?/ + sr^*\ 


(13) 


contains only s as an unknown. The solution of (13) , expressed in a 
form to minimize computational round-off errors, is; 

s == -(B + V^)/2A, 

where 

RAD = b2 - 4AC, 


(14) 


A = F + (1 - F) z 2, 

1 

B = = 2(x X + y y ) F + 2z z , 
I s 1 s 1 s 

C = (x^ + y^) F + 2^ - r^, 
s s s p 

F = r^/r^ , 

P eq 


[ERTOST] 


J 


With the solution of s, it follows from (1) and (10) that 


tg « R(6g,3)(r^ + sr^) * (x^,y^,Zj,) , [ERTOST] (15) 


and hence 


= tan l[z^/((l - e^)^ x| + y|)] 


[ERTOER] (16a) 
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“ tan-1 [y^/x^], [ERTOER] (16b) 

where 

9j^ = geodetic latitude, 

= longitude 

~ ~ ^^^th's eccentricity squared. 

Equations (9 ) , (15 ) , and (16 ) constitute the basic - satellite image 
coordinates to earth coordinates — transformation* 

C* Navigation Optimization Procedure (Subroutine ATTTUD] 

Navigation of the satellite image' data base consists of finding ' 

a time dependent transformation to predict the earth coordinates from 

the associated satellite image coordinates. This is accomplished by 

using landmark measurements from the data base and the model discussed 

above to determine the optimal transformation in a least-squares sense. 

A landmark measurement consists of the earth coordinates, (8 ), the 

JL tf 

associated image coordinates, (LIN,ELE), and the time, t, at which the 

landmark was imaged. 

Let 

. th 

= 1 parameter to be optimized, 

R. (a.) = transformation associated with the a.. 

Ax X 

^ 1 th 

“ unit pointing vector of the k landmark derived from (9 , 

^ Xj t ic 

r, = unit pointing vector of the k landmark derived from (LIN,ELE), , 

K. 

S(ap = - R^(ap ^ '||2 [S] (17) 

k 

which is a sum over all landmark measurements included in the optimization. 
The navigation is complete when a set of parameters a, are found (call 
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them a^) which minimizes S, The optimized values are then used in 
the transformation to predict (0^, <!),.) from an arbitrary (LIN,ELE). 

Xi L 

As an example, consider equation (5a) , 
where [BCTOPE] 

and let us assume that the attitude telemetry data is known with a 

s 

reasonable degree of accuracy but that the orientation of the PF frame 

relative to the BC frame is not. Thus, we wish to optimize 0,9,0 

y r p 

in the transformation Kpj.* Using the landmark measurements and the atti'tude 

telemetry data, calculate r^^ as given by equations (1) to (4a) and r„„ 

150 irr 

from equations (8) for each landmark; therefore S in this case takes the 
form 


and the values 0,0,0 which minimize S are the ones then used in 
y r’ p 

equation (5b) . 

The method used to minimize S (u^) is an iterative procedure which 
uses a modified Newton’s method. 

Let 

a = {«.} for convenience, 

X 

value of a resulting from the n iteration, 

N = total number of parameters 

H = N XN matrix (Hessian) whose ij component is 
VS = gradient of S, (VS)^ = 3S/8a^* 


The iterative procedure works as follows: 

1) Start with m = 0 and increase m by 1 until a value M is found 


such that 
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S(a^) - S(a^ - 2 H ^VS(a^)) (,4)2 ^ H“^VS(a’^) • VS(a"), [LS] 

2) When the inequality is satisfied, set 

= a“ - 2~^ H"^VS(a”), [MCNMIZ] 

3) Check to see if the following convergence criteria is met 

ls(a^’^^) - S(a”)l ? 10-l®ls(a^)| . 

If not, set -f- and go to l);,if yes, then = a and the 

procedure is finished* 



APPENDIX B. ANALYSIS OF METHOD WdlCU DETERMINES ATS-6 SSP IMAGE COORDINATE 
DISPLACEMENTS BETWEEN SUCCESSIVE IMAGES RESULTING FROM ATTITUDE 
CHANGES 

This appendix presents an analysis of the technique used to calculate 
the ATS-6 Sub satellite Point (SSP) image coordinate changes between successive 
images resulting from attitude changes during the image -scan time. Appendix C 
shows how these measurements are used to account for the attitude changes 
in the process of computing accurate cloud displacements* 

1. Method 

Let T]^, T 2 designate two successive ATS-6 data images where is the 

^’earlier” of the two images and E^) (Line scan and Element numbers) 

are its SSP image coordinates (Fig. B.l) determined from the ATS-6 navigation 

model* Let (L,E) be the Tj[ image coordinates for a point on the right 

earth edge and AE , AE the measured displacements along line L of the 
R L 

Right and Left T 2 - earth edges relative to the Ti earth edges. (These 
measurements arc obtained in practice on McIDAS using the infrared ATS-6 
data images and an image-matching technique which is constrained to measure 
displacements of the earth edge only along a scan line* It is worth 
emphasizing that the image -matching method indeed measures displacements 
of the geometrical earth edge and not features near it - such as clouds* 

This is not surprising since the greatest contrast is between earth and 
space - not within features near the earth’s edge.) 

The object, then, is to compute the displacement coordinates (AL,AE) 

r 

of the T 2 “SSP relative to the Tj- SSP at line L using the measured values 

of AE ,AE, . 
r L 
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FIGURE B.l* Earth-edge Displacement Measurement 
Geometry 


Let points 1 5 be on line L containing the earth edge points "indicated 
in the figure, and points 0 and O' the T^-SSP and T2-SSP, Furthermore, 
point 5 is the bisect point for the T2— chord coinciding with line L, and 
both and X2-earth circles have equal radii a* 

Letting be the distance (always non-negative) between points 

'i and j, AL,AE can be derived from two simple geometrical identities: 

' The first is 



(E + AE ) - (E + AE) - (E -f AE) - ( 2 E - E + AE, ) , 

c c c L 

solving for AE, 

AE = (AEj^ + AEj^)/2 . (1) 

The second identity is 

“ 1*2 ~ ~ ^3 ~ ’ where 


( 2 ). 
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• ^4 - "2 - 2 - ^ 0 ->' - '^5 - 

- 2V(ppr^)2- (pprp^)2 

= 2 V(L - L^)2 + (E - E^)2“- [L - (L^ + L)]2 , and (3) 

P 3 - \ = 2(E - E^) . (4) 

Substitution of (3) and (4) into (2) and rearranging yields 

6 = Vx2 + - (Y - AL)2 - X , where (5) 

6 = (AEj^ - 4Ej^)/2 

X = E - E 

c 

Y = L - L . 

c 

Since 6 and (L,E) are obtained from the AE ,AE measurements and a nominal 

R L 

value used for (L ,E ), AL can be solved in (5) to yield 
c c 

AL == Y + VRAD , where ( 6 ) 

RAD = Y^ - 2 X 6 - 6 ^ . 


The choice of sign in ( 6 ) is determined by substituting the expression for 
6 in Eq. (5) in RAD; the result is 


AL = Y + V(Y - AL)^ = Y + , Y > AL 

« Y - , Y < AL . 


(7) 


The peculiar condition implied by (7) that AL must be knoxm before AL can 
be calculated is really not a problem since generally AL is small (+ 5 lines) . 
and |y| is usually » |al|. For most cases then, the conditions are 

AL « Y + VrAD , Y < 0 , - ( 8 ) 


RAD=Y^-2X6-62 ^ 



B-4 


For the case where |y[ approachs (AL) in value there are other problems; 
these are discussed at the end of this appendix. 

Equations (1) and (8), then, define the displacement in image coordinates 
of the T 2 -SSP relative to the T^-SSP for line L. By repeating this entire 
process for other scan lines the SSP--shifts (AL,AE) as function of T^-line 
position are obtained. In practice the right and left earth--edge displace- 
ments are first measured over the entire range of scan lines of interest. 
Curves are then separately fit to the right and left edges measurements 
resulting in a vs L and a AE^^ vs L curve over the scan-line range of 
interest. The curve values themselves are then used in (1) and (8> to 
compute ALjAE curves. 



C-1 


APPENDIX C. ANALYSIS OF ALGORITHI^ IfflICH ACCOUNTS FOR RELATIVE ATTITUDE 
CHANGES IN SUCCESSIVE ATS-6 DATA IMAGES USING EARTH-EDGE 
SHIFT MEASUREMENTS 

This appendix provides a semi-rigorous analysis o£ the algorithm used 
to account for the relative attitude changes in a sequence of ATS— 6 data 
images using the earth-edge measurements derived by the technique discussed 

•Cl 

in the main portion of the report and APPENDIX B. In effect, this appendix 

shows why the algorithm "works”; the following appendices discuss in more 

detail the limitations of and the errors associated with this method. 

Consider a sequence of n ATS-6 images designated by . . .. T , 

12 n 

from which the displacements in earth coordinates of a feature (cloud) 
between successive images are to be determined* Let 

(L^,E^) — image coordinates (Line, Element) of the cloud feature 
for the image, 

= f„(L.,E,) - unit pointing vector in Picture Erame coordinates 

JrrX irr x X 

\ 

derived" from (L^,E^) using the ATS-6 scan-camera geometry, 

= unit Local Vertical pointing vector associated with 

- 3x3 rotation matrix which transforms without error 

into at the time the feature was scanned, i.e. 

PEx ^ 

’^PFi "" \ "^LVl • 


1* A Two-Image Sequence 

Consider now the first two images T^ , in the sequence whose relation-- 
ship between LV and PF coordinates are given by 


^PF2 - ^LV2 ’ 




( 2 ) 

(3) 
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rhe relationship of the cloud feature in LV coordinates for image 
can be- related to as follows: 


^^LV2 ^LVl 




( 4 ) 


where is the displacement of the cloud in LV coordinates from to 

and results from two effects: 


(i) apparent change in the earth's orientation over the time 
interval caused by a non-zero angle of inclination in 

the satellite's orbit (the eccentricity of the ATS-6 orbit 
was so small (~ lO"**) that no significant - change in angular 
size occurs) . . ,u'- ' - - 


(ii) motion of the cloud relative to the earth's surface. 


Now the orientation of the earth relative to the LV frame depends only on 
the orbit and the earth's position relative to celestial coordinates; thus 
xf r^^^ and r^^^ were accurately known, the displacement of the cloud over 
the time interval associated with these two vectors could- be calculated 
with an accuracy limited only by the equations describing the d 3 mamical 
relationship between the satellite orbit and the earth^s position relative 
to celestial coordinates • It is assumed that for ATS-6 that this transformation 
produces negligible error. 

Returning to equations (2) and (3) , the rotation matrices differ 

from each other slightly because of an attitude change of the PF frame 

relative to the LV frame between T and T : they are related to each other 

by an infinitesimal rotation I + where I is the unit matrix and 

an antisymmetic matrix whose elements are the small-angle differences between 

the PF axes at times T and T . Thus 

1 2 
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\ = a + Rj ; (5) 

Slibstituting (4) and (5) in (3) and expanding 

The first term on the RHS of (6) is r „ , the last term is negligible 

rr 1 

(~10^ times smaller than the first order terms) , the second term is the change 
in PF coordinates due to an attitude change and the third term is the change 
in PF coordinates due to the two changes in the LV frame discussed above. Using 
(4) to combine the first and third terms, equation (6) then can be rewritten as 


^PF2 ^1^LV2 *** ®12^PFl 


(7) 


Comparison of (7) with (2) implies that if and were known, 

and therefore the cloud displacement “ could he computed with 

a high degree of accuracy* 

A good approximation to is obtained from the earth-edge displacement 

measurements which provide the displacement in image coordinates AL^^(L), 

AE, (L) of T relative to T, as a function of line number L. The line shift 
12 2 <1 

AL^^ is proportional to a small change in roll PE-X axis and 

the element shift is proportional to small changes in pitch about the 

PE-Y axis, i.e* 

ALi^ = 4R^j,/p , (8a) 

AEi 2 = AP^^/p , (8b) 

where p is the angular size of a pixel = 1.45 X 10~** radians. 

^ T 

Letting - [X, Y, Z] where T is the transpose, the explicit form of 

and the last term in (7) is 
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1 

o 

o 

o 


"x“ 

E r_„ = 

0 0 AR 


Y 

12 PFl 





-AP -AR 0 _ 


_z. 


= [APZ, ARZ, - (AFX + ARY)]^, (9) 

where it is understood that AR, AP refer to the to changes in attitude 
at the times the cloud was scanned. 

Note Ei 2 contains no non-zero yaw angle element (small rotation about 
the axis passing through the image SSP) since there is no way that the earth- 
edge displacement measurement technique can provide this angle; however, the 
analysis given, in Appendix E of reference 3 shows that the effect on the 
accuracy of wind measurements (displacements) for ya\<r changes is generally 
small. Analysis of landmark measurements and the ATS-6 wind sets sho^m. in 
this report, indicate that for the data images studied thus far, the yaw 
changes are negligible. 

With and hence the last term in (7) determined, equation (7) can 
be rewritten as 


r' - = r. 


E. r„„. = R r. 


PF2 PF2 12 PFl 1 LV2 ’ 


( 10 ) 


where ^pj .2 interpreted as the cloud pointing vector with the attitude 
changes removed. 


2. Computing in Practice 

Cloud displacement measurements and computation of accomplished 

as follows: 


(1) the Tj image coordinates (L^,E^) of the cloud are recorded. 
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(2) the total displacement of the cloud in image coordinates 

T T 

(AL , AE ) from to are measured on McIDAS using an 
image matching technique^ 

(3) from previously determined earth-edge displacement parameters, 

the displacement coordinates relative 

T T 

attitude changes are computed and subtracted from AL , Ae , 
These differences are added to , E^ and the associated 
PF vector computed. That this vector is a good approximation 
of r’ ^ is shown below. 


By definition. 




( 11 ) 


^PFa = ^FF<^ + ^ = ^FF»a'®2> 

Expanding ippO-j + “ *^12 ^^ 2 ^^ 

T 

+ AE ) keeping first order terms, we have 


"ff<S \ ^ "pf'S-^ 2> ^'^la - 4 ^PFtta'^2> “12 ' 


From the geometry of the ATS-6 scan-camera, 

A T 

Xpp = [-cosXsind, -sinX, cosXcosfiJ 
» [X, Y, Z]"^ 

where 

X = (L - L)p - mirror step angle, 

'' c 

fi = (E - E)p = mirror sweep angle, 
c 

(L^,E^) ~ image center line and element value, 
p = angular size of a line or element. 


( 14 ) 


( 13 ) 
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Taking the partials of r^j, in (14) with respect to L,E, and retaining 
first order terms with the partials evaluated at (L 2 ,>E^): 


3L 3E ^PF "^^12- 

= [0, PZ, ~PY]^AL^^ +' IPZ, 0, -PX]^AE^2 , (15) 


substituting (8a) and (8b) into the above expression and adding the two 
vectors in (15) yields 

[ZAP , ZAR , - (XAP + YAR )]^ , 

12 12 12 12 ’ 

which is equivalent to (9); therefore (15) is the same as Since 

the first term in (13) is have 




r “ E r = 

PF2 12 PF PF2 


by 


(16) 


comparison with (10) . Thus the method discussed above correctly yields 

^ f 

""pF2* 


3. Transforming PF to LV coordinates 

With the attitude changes between T^ and T^ remqved, we have 

’^PFi "" ^I’^LVl 

^PF2 “ ^^LV2 


' (17) 
(18) 


' Let R be the transformation from LV to PF coordinates derived from 
T^ landmark measurements in a least squares sense assuming a constant 
attitude over the time interval the T^ landmarks were scanned, i.e. R 
is a constant matrix. R then will differ slightly from R^^ due to small 
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attitude changes over the scan time and, in addition, will differ from 
(assumed to be the error-free transformation) due to mirror scan 
nottlinearities. Similarly, as above, we can relate R to R^^ by an infinitesimal 
transformation 


R = R, (I - e) , (19) 

1 

where e is the infinitesimal antisymmetric error matrix which is a function 
of line and element. It is important to realize that the elements of e 


contaxn absolute e rrors which the AK-6 model can not acc ount f or. . - - 

However, our previous work using landmark measurements has provided us with 
bounds on these errors and the rate at which they change as a function of 
image coordinates. The worst-case estimates correspond to about +2.5 pixels 
(+3.6 X 10~** radians). The errors tend to be oscillatory as a function of 
line or element position with a period of 100 to 200 pixels. Thus for 
typical feature displacement measurements (corresponding to tens of pixels) , 
e can be considered to be a constant locally. (It would be appropriate to 
remind the reader at this point that the absolute mirror scan nonlinearities 
as a function of image coordinates repeat from image to image.) 

Taking the transpose of (19) and operating on (17) and (18) 




A -y 

= r + r 
LVl e 


^ ^PF2 **" ’^LV2 ’^e ’ 


where 



(20b) 

( 21 ) 


Application of R then, results in the correct LV vectors to within the same 

additive vector constant. It is sho^m in Appendix C of reference 3 that this 

constant has a negligible effect in computing cloud velocities or displacements. 
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In addition location errors on the earth due to is given in this 
appendix* 


4* General Case of n- image Sequence 

The case for images ... simply a generalization of the 

2-image case. The error-free LV to PF frame matrix for can be related 

to R by the image pair relative- attitude transformation matrices E, , - 
1 1 ^ l*hl 

where 


R == (I + E - ) •*• (I + J(I 4- E,^)R 

n n-l,n' ^ 23 12 1 


■ 


where only first order terms are retained. Thus the analog to (10), generalized 
to n images} is 


’^PFn ^PF ^i,i+lj’^PFl ^i^LVn ’ 


where in practice is computed by evaluating the expression 


r =5 r 
PFn PF 


( ' n-1 n-1 . \ 

L - 2 AL^ E - Z AE. (23) 

n . i,i+l 1+1 n , 1 , 1+1 1+1 y ' * 


where L. ,(L _), E, ..i(L.,-) are evaluated- from the T, to T,,_ earth-edge 

1,1+1 1+1 i,i+l 1+1 1 i+1 

displacement measurements evaluated at L where L.,- is the T. ,, line number 

1+1 1+1 i+1 

/ 

of the "center of gravity" of the feature being tracked. 

Application of (19) to the results in 


^ t ^ I 

^■TT74 = ^TT7- + ^ 

LVi LVi a 


i = 1, 2, ... n 
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APPENDIX D. METHOD OF OBTAINING MIEPOR-SCAN-OFFSET CORRECTION CURVES 

The nature of the mirror-scan nonlinearity problem has been discussed 
in previous reports. In order to correct the alternate scan offset caused 
by this effect, we need a table of offset values as a function of element mraiber 
(i.e. a AE(E) function). For our initial efforts to produce a AE(E) function, 

a 

we used the McIDAS cloud-tracking program to measure the displacement of 
a feature seen in the odd number scans of an image to its position in the 
even numbered scans of the same image. The method described here is a . • 

somewhat more automated scheme and does not require viewing the ATS-6 images. 

The method has been applied to an IR image (74195 173134Z) and gives good 
agreement with the old method with less scatter of individual points about 
a polynomial least-squares fit curve (Fig. D'.l. ) . 

In the automated method of computing the AE(E) function a correlation 
value is computed for the match between a small segment of an odd scan and a 
shifted segment of an adjacent even scan. The amount of shift which gives 
the best correlation for that small line segment is taken to be the AE 
value for that scan line and element location. Values are computed over 
many scans and elements, then averaged over the scans. The result is a table 
of AE values as a function of E (element). The table values are then smoothed 
by using a least squares fit pol 3 momial. In more detail the method is as 
.follows: 

Let: 

e = element number (1 ^ e £ 2400) 
s = scan number (1 ^ s ^ 1200) 

p^(s,e) = pixel digital value at a given even-scan, element position 
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p^(s,e) “ pixel digital value at a given odd--scati, element position 
6 - an element shift value 

t — a small number defining the length of the scan segment used 
for a correlation 

C(s,e,6) = a measure of the match between line segments from adjacent 
odd and even scans. 

Now define C byt 

e+t 

C(s,e,6) = 2 [p (s,e’) ~ p (s + 1, e* + 5)]^ (s odd) 

e’=e-t ° ® 

The element shift, AE, for a given point on the image (s,e) is the value 
of 6 that gives a minimum value to C; 

C(s,e,AE(s,e)) _< C(s,e,6) for all values of 6 

After an array of these values is generated, a weighted average over scans 
is taken: 

I w(s,e)AE(s,e) 

AE(e) = 

I w(s,e) • 

The weighting factor is the range of brightness values from the odd scan 
jused in computing C. 

w(s,e) = MAX(p^(s,e*)) - MIN(p^(s,e')) 
where e-t^e*_5.e + t. 

The reason for using this weighting is simply that high contrast features 
should, in general, produce a correlation value of more significance than 
features of nearly uniform brightness. 



The computer program to do this used the values: t - 7 and ^ ^ ^ 1^* 

Values of AE(s,e) were computed for every eight scans for the middle third 
of the image (s = 399, 407, 799) and for every 10 elements across most 

of the width of the image (e = 100, 110, 2300), Using .this line and 

element range some AE values will be coii 5 )uted for points off the earth. TJhen 
off-earth and on-earth points are averaged together the weighting factor will 
give only a very small or zero weighting to the off-earth points. For some 
element values no on-earth points are encountered. We found that these 
points could he determined from the table of AE(e) values. There is a 
discontinuity between the on-earth and off-earth values. 

The AE(e) function is then used to create a corrected image by shifting 
even scans. An improved method of correcting the. alternate scan offset has 
also been developed. In the old method a fixed shift value was used across 
the width of a McIDAS image (672 pixels). Although this gave fairly good 
results, there could be a one or two element alignment error near the edges 
of the image. The_new method, as before, uses the odd scans as a fixed 
reference. However, the required amount of shift is computed, or looked up 
in a table-, for each pixel in the even scans. Thus the amount of shift varies 
from one side of the image to the other. Pixels are dropped or doubled, as 
appropriate, between regions of different shift values. 
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FIGURE D.l. Mirror-Scan Nonlinearity Curve AE(e) using newly developed procedure. 
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Appendix E- 


Saurce Code Listings (FORTRAN) 
of Programs and Subroutines 


This appendix contains FORTRAN source listings for the 
main programs and subroutines for the ATS--6 image correction 
and navigation system. For main programs both AOIPS and McIDAS 
versions are given. For subroutines only the McIDAS version is 
given as the AOIPS version would be identical. 



ATSF2/A0IPS 


0XM£NSI0i4 MIM(IO) 


711 

FORMAXI///// 

712 

FORMMC 

74H 

713 

FOkKMC 

?4H 

714 

FURM4XC 

74ri 

‘715 

FORrtftTC 

4 * 

7 4H 

716 

'PURHA'TC 

74ri 

717 

FOrMaIC 

74H 

718 

FOkMAK 
IT = 5 

74H 


ATS 6 PROCESSING 
1 iNllIALIEE NAVCOM 


J 


2 GENERATE BGEMENT OFFSET DATA FROM E.H.X, 

3 CURVE FIX TO OFFSET DATA 

4 READ image SEGMENT FRqM E.H.T. 

5 CURVE Fir TO EARTH EDGE DATA ' 

6 RUN navigation 

7 ExH ATSfe processing 


1 

2 

3 

4 
b 
6 
7 

990 

999 


WRITE(IT.7in 
• rtRlTSCIT,712) 

WRlTEtlT,713) 

rtRirE(lT,714J 

WRHEIH',715) 

wRIXEClT.7163 

wRlXECIT.717] 

*»RlXEtiX,718) 

NI.=8 

MENUsl 

CAUL INCOHC* > ,0.NL,1 ,M1H,HEN0,1*1J 
■ IF t«ENo.LE.O) GO TO 990 
IJsMiN'm 

GO TO tl.2,3.4,5,b.7J .IJ / 
call REQOESIRADsOC A6INT2' J) 

GO TO 999 . 

CALL REQUES(RAD50( ' 0FST62* J J-' 

GO TO 999 '' / 

CALL R£QUES(RAD50C'OFSTF2' )J 

GO It) 999 . 

CALL REQ0EStRAD50C'LATSF2' )) 
oO TO 999 

call R£QUE5(RA050C'EDGFT2')3 
GO rO 999 

CALL REGgES(RAO50('ATSNV2' 7) 

GO 10 999 

call RBOUeS CRAD501 ■ rfET2' ) ) 

GO TO 999 
CONTINUE 

CALL REQObS(RAD50t'rtE12' J.) 

CONTINUE 

END. 


i 

) 

7 

7 

7 

7 

•7 
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ATSW/AOIPS 


C 


2 

10 


730 

40 


708 

50 

51 

70b 


780 

100 




program to NAVlGATe, Al&b IMAGES FROH LAKOHaRKS, 

IHIEGER+4 lLAr»lljUN,liD,lTIME 

OImEmSIOn PTiMEClSO) . AL1MU50J , AE.DE( 150) .ALATt 150) , A0t)«(160) 
DIMENSION IliMEClsO) ,101N(150) .IBOEC150) .lt<ATU50) ,n.0NU50) .ICOOE 
♦C150) 

DXMeNSIUH (UN(10).«0UXC24) 

Cuh«ON/iUVcUM/NA^«, InAV.lYR, lOAI.TOTblW.OEGIilN.TOTlEL.DEGEljE.PlCl,! 
=^l'l ,PlCE,LE,I'MPSCI»,iUIB, lOuAl ,1M,R1X, RlI,Rli,RlDX,KlDI ,’Rl OX , PITCH . ROL 
X Art .P'llH (3) .TMn{ 3) , 1MX(3) . NbCOEP 1 2 ) . nRCOeF ( 2 ) , SCObO ( 2 ) , 5CLL1 C 2 ) . 
4EbCO£F(ll,2) .SCbhOta) , SCbRl (2 ) , ERCOEF ( 1 1 , 2 ) . W AfaCEF.SCLASO , SCLASl . 

♦ lEb£M«,lEbEfJX,ASCUEfC16) 

DATA Pi/3. 1415025535/ 

Data IT/5/.GP/6/. 

RADdEg=P1/180.0 

bU«=10 

OPEMCONir=l.ON,«AHE='DtJO: 1350, 62jAfabMK.S. DAT* .TlPh='ObD' .RBAOOKLI) . 

CONTliNOE 

Nb=0 

CONriNOE 

«b=Nbtl 

REAOCbUN,730.ENO=40)IlD,lTlM6;(N£.).lCoDECND) .ILlN(NL) .lEliElNL) ,XbAI 
+CNL) .IbGNCNb) 

FURMATtZX.SIO) 

GO TO 10 
CONIINOfi 
Nb = Nti-l 

CbOSEtOAlTsliUN) 

KC=2HDb 
00 50 1=1, «b 

WRJTb(bP,708)lTXMEU) .ICOOECl) ,lbl«CI) ,I£LEtI),lbAltl),XbONa)- 
f 0RMAmx,8X9) 

COMPlNOE 

CONTXNOE. 

hRlTE(IT,706)XDAy,ffli 

FOrMATUX, 'DA1=' ,X5, • number of IiAnDharKSs’ ,13) 

NbrtK=0 


DO 100 1=1, NL 

If (MOOCICOOElIj/iOO,10),NE.O) Go TO 100 
NbMK=NbMR+l 

PTI«EtNbMK) = FTIME(irit4ECl) J 
AblMCNbPlK)=FbOATtlblN[X)) 

AEbE(«I,MK) = FuOAmEbECX)l 
AbA'l (NbrtlC)=FLAbO(lbArCI)) 
AbON(HbMK)=FbAbiKlbOhtI)) 

rtRXiblbP, ?80jPl l»ECNbMK) ,ACjXN(NbMK) ,AECiECNl,HlO, 
4AbAl CNbMR) , AbONlNOMX) 

FORHAl Cli,6Fl5.5) 

CONxXi'iOE 

.CbOSE(ONir = t,P) 

. IFCSOMtC.tiE.O) GOTO 990 

CAbb AiXrOoCPriHE, AbXN.AtbE, AbAT.AbON.NtiMK) 

PD=PlrCH/KADOEG 

RDsRObb/RADOEG 




ATSNV/AOIPS 


lOsYAW/HADObti ' 
wRirblll, 70iJPD,RO,lO 

701 tORHAlCiaO.iX, ‘PlrChs' ,fcl6.9. ' R01iL=' ,E16.9. ' XAw= • ,E 1 6.9) 

196 CUNllNOE 

Oa 200 1=1. KL 

PTH = FHrtt.ClllHE(X)) 

XliAsFLALUClliAlCI)) 

XL0 = Fl,MbDCl6O«(l)J 

CALL LS(PrK.XLA.XLO.XblN.XELC) 

RIiIh = Uii«tI)-X£jIN 
Keb£=l£Lfc,tlJ-X6.L6* 

HKIIt:(LP, /lO)iriMECI) .iCODttl) .RtilN.REIiE 

200 continue 

201 CONTINUE 

710 f0RHAr(3X,l6,5X,l5,5X,2F10.2) 

GO TO 999 

990 CONTINUE 

999 continue 

CAOL REQUES[RAD50C'ArSF2')) 

END 



★★★A** ATSNV/MCIOAS ****** 


SELTfL 
ELT007. 
000001 
000002 
000003 
000004 
000005 
000006 
000007 
000008 
000009 
000010 
GD0D11 
000012 
000013 
' GOOOU 
000015 
000016 
000017 
000018 
000019 
000020 
000021 
000022 
000023 
000024 
000025 
000026 
0D0027 
000028 
0D0029 
000030 
000031 
000032 
000033 
000034 
000035 
000036 
000037 
000038 
000039 
000040 
000041 
000042 
000043 
000044 
000045 
000046 
000047 
000048 
000049 
000050 
000051 
000052 
0-00053 
000054 
000055 


AF*ATSNV/HC 
RLIB62 12/2 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 ' 

000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
GOO 
000 
000 
000 
000 
000 


IDAS 

2-1 7:01 :39-(0,) 

$JOB ATSNAV U3200 
^OPTION .8,9,20 
SFORTRAN 

SUBROUTINE MAIN 

C PRObRAM TO NAVIGATE ATS6 IMAGES FROM LANDMARKS* 

LOGICAL OPTION 
INTEGER POS,GL£B 
EQUIVALENCE (I COM ,NAVN1 
DIMENSION IC0MC1) 

DIMENSION PTIMEC150);alIN( 150) jAELEdSO) ,AlAT(150) ,AlON( 150) 
DIMENSION IT lHE(150)tlLlN(150), IE LE (150), ILAT(l50)fIL0N(150),I CODE 
*(150) 

DIMENSION MSQC10) 

DIMENSION MIN( 10) ,M0UTC24) 

DIMENSION 6LES(672) 

. COMMON/NAVCOM/NAVN,INAV,lYR,IDAY,TOTLINtDEGLlN,TOTIEL,DEGELE,PlCLI 
*N, PICELE ,TMPSCL,IO YR, lODAY fTM,R1XtR1Y,R12,RlDX,RlDY,KlD2,PiTCH,R0L 
*L,YAW,P7IM(3) ,TMN(3),TMX (3) ,NLC0EF(2) , NRC OE F ( 2) , SCLLO ( 2) , S C LLl C2> , 
★ELCOEF (11 ,2) ,SCLR0(2) ,SCLRU2) , ER C 0£ F ( 1 1 , 2) , NASC E F , SCL AS 0,S CLAS 1 , 
★IELEKN,IELEMX, ASC0EF(16) 

DATA Pl/3. 1415926535/ 

DATA MSG / 6 HD E f AID, 8*0 A 
DATA ICODE/150*0/ 

DATA M1N/6HATSNAV,8*0/ 

CALL IQCMIN) 

NDAY=MINC1) 

J0UT=1 

IF(OPTION (MIN(4) ,3H P)) 

IF (.NOT. OPTION (H1N(2) ,3h 
IPT1=MIN(5) 

IPT2=:MIN(6) 

IPT3=MIN(7) 

NAVDAY=MOD (NDA Y , 1 00000 ) 

ISS==NDAY/100000 
ISS-(lSS/2)*2 
KDAY=ISS* 100000+NAYDAY 
RADDEG-Pl/180.0 
IFCMINCS) .NE.3HNEW) GO TO 2 
DO 3 1=^1,203 
► 1COM(I)=0 

CALL WCOM 
CONTINUE 

CALL GETNAV(KDAY, lEXXST) 

1FC1PT1.N£.0)PTIMC1)=FTIME(IPT1) 

1F(1PT2.N£.0)PTIMC2)=FTIMECIPT2) 

IF(IPT3.NE.0)PTIMC3) = FTIME CIPT3) 

IFCIEXIST.GE .1G)bO TO 990 
CALL -H£DD£RCKDAY,6LEB , POS) 

1 F(6LEb(P0S*l ) .NE.KOAY .OR.GLEB CPOS+2) .LT.O)GO TO 990 
1ND£X=6*GL£B (P0$+2)+88 
NL-0 

10 CALL SCRA(30, INDEX) 

CALL RbADW<30, 672 ,6LEB) 

' DO 30 <1 = 1,661,6 
KIND = GLEB CD/1 OOGOO 


J0UT=2 

A) .AND . .NOT. (MIN(2> .EQ.O 


))6D TO 195 


OF FOOB QPALIW 



****ft<t ATSNV/fnCIDAS 


****±it 


000056 

000 


IF(KIND.EQ.O) 60 TO 40 

000057 

000 


IFCK1N0.NE.1J 60 TO 30 

000058' 

000 


IF (NL.bE.1 50) GO TO 30 

000059 

000 


NL=NL+1 

000060 

000 


iTIH£(NL> = 6L£8a + 1 ) 

000061 

000 


I CODE (KL)=«OD( GLEB (I) ,100000) 

000062 

000 


ILIN (NL)=GLEB(H-2) 

000065 

000 


lELE (NL)=6LEB{I+3) 

000064 

000 


ILAT (HL)=GLEB< 1+4) 

000065 

000 


IL0N(NL)=GLEB(I+5) 

000066 

000 

30 

CONTINUE 

000067 

000 


J J=GLEB(672) 

000068 

000 


IFO J .LT.O) 60 TO 40 • 

000069 

000 


1NDEX=6*J J+88 

000070 

000 


60 TO 10 

000071 

000 

40 

CONTINUE 

000072 

000 


KC=2HDL 

000073 

000. 


IF(.N0T.0PTI0N(HIN(3),3K L)> GO TO 51 

000074 

000 


DO 50 1=1, NL 

000075 

000 


ENC00E(132,708,K00T)KC,NDAY,ITIHE(I),IC0DEa),lLINII),IEL£(l), 

000076 

000 


*ILAT (I),ILON(I) 

000077 

000 


CALL TP< JOUT,HOUT) 

000078 

000 

708 

FORHATC1X,A2,8I9) 

000079 

000 

50 

CONTINUE 

000080 

000 

51 

CONTINUE 

000081 

000 


ENCODE C132, 706, HOU DID AY, NL 

000082 

000 


CALL TP(JOUT,HOUT> 

000083 

000 

706 

F0RKAT(1X,'DAY=',I5,' NUMBER OF LANDHA8KS = ‘',13) 

000084 

000 


NL«K=0 

000085 

000 


DO 100 1=1, NL 

000086 

000 


IF(KODaCODE(I)/1O0,10).NE.0) 60 TO 100 

000057 

000 


NLHK=NLHX+1 

000088 

000 


PTIME (NL«K)=FTIMECITIME( I) ) 

000089 

000 


“ALIN CNLMK ) = FL0AT(ILIN (I) ) 

000090 

000 


AELE (NLMK)=FLOAT tIELE (I) ) 

000091 

000 


ALAT (NLMK)=FLAL0 (ILAT( I) ) 

000092 

000 


ALON CNLKK) = FLALO<ILONtI) ) 

000093 

000 

100 

CONTINUE 

000094 

000 


I F(NLHK.LE.O) GO TO 990 

000095 

000 


CALL ATTTUO(PTIKE, ALlNjAELE, ALAT, ALON, NLMIC) 

000096 

000 


PD=PITCH/RADDEG 

000097 

000 


RD=ROLL/RADD£G 

000098 

000 


YD=YAW/RADDE6 

000099 

000 


ENCODEC132,701,HOUT)PD,RD,YO 

000100 

000 


CALL TP(JOUT,«OUT) 

000101 

000 

701 

F0RKAT(1H0,1X,'PITCH=',E16.9,' ROLL = ' , El 6. 9 , ' YAK=',E16.9) 

000102 

000 


CALL WCOH 

000103 

000 


HSQ (3)=KDAY 

000104 

obo 


MSQ(4)=ILAL0(PD) 

000105 

000 


MSQ(5)=ILAL0(RD} 

000106 

000 


«SQ(6)=ILAL0(YD) 

000107 

000 


CALL SQ(HSQ) 

000108 

000 

195 

CONTINUE 

000109 

000 


IF ( .NOT.OPTION (HIN (2) ,3H R) . A ND . .NOT , <HI N (2) .EO .0 ))GO TO 201 

000110 

000 


CALL GETNAV(KDAY,IEXXST) 

Q00111 

000 


■ DO 200 1=1 ,NL 

000112 

000 


PTH=FTIHE (ITIHEII) ) 
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000113 

000 


XLA=FLALO(ILAT<I)) 

000114 

000 


XLO=FLALO(ILON Cl>) 

000115 

000 


CALL ES(PTM,XLA,XL0,XL1N,XEL£> 

000116 

000 


RLIN=-ILIN(I)-XLIN 

000117 

000 


RELE=I£LE(1)-XELE 

000118 

000 


ENCOD£n32,710 ,HOUT)ITIHE<I> t ICOt> E< I > , RUN, RELE 

000119 

000 


CALL TP( JOUT.HOUT) 

000120 

000 

200 

CONTINUE 

000121 

000 

201 

CONTINUE 

000122 

000 

710 

FORMAT(3X,I6,5X,iS,5X,2F10.2> 

000123 

000 


RETURN 

000124 

000 

990 

CONTINUE 

000125 

000 


CALL EHESSC3HREQ,NPAY) 

000126 

000 


RETURN 

000127 

000 


END$ 

000128 

000 

$FI LEHA 

000129 

000 

DELETE ATSNAV,GORP 

000150 

000' 

SINCLUDE HEDDER 

000151 

000 

SCATAL06 

000132 

000 

NA«E= 

ATSNAV»5rRtW,D 

000133 

000 

TYP E= 

EG 

000134 

000 

LIB=ATS FLBtLL 

000135 

000 

BEGIN 


000136 

000 

SEO J 

V. 


END ELT. 

3HD6,P ****** ATTTUD ****** 



****** atttud/hcioas ****** 


SELT,L AF.ATTTUD/MCIO AS 


ELT007 

RLIB62 12/22-1 7 

:01 :41-(0,) 

000001 

000 


SUBROUTINE ATTTOD CPU M E t ALIN , A ELE t ALAT t ALON , NLHK > 

000002 

000 

c 

SUBROUTINE TO CONPUTE ATS6 ATTITUDE FROM LANDf^ARKSo 

000003 

000 

c 

YAW, PITCH-, AND ROLL VALUES FOR A SHALL OFFSET OF THE PICTURE 

000004 

000 

c 

FRAME COORDINATES FROM BODY CENTERED COORDINATE SYSTEM. 

000005 

000 

c 

7 JUNE 1977 G. C. CHATTERS 

000006 

000 


DIMENSION PTIHE (1 ) , ALIN Cl) ,AELE(1 ) , ALAI Cl ) , ALONCI) 

000007 

000 


DIMENSION X(3,150) ,YC3,150),TIKEC150),PRY(3) 

000008 

000 


COH«ON/NAVCOH/NAVN,INAV,IYR,IDAY ,T0TLI N , DEGL IN ,T0TIEL, DEGELE , PICLl 

000009 

000 


*N,PICELE,TMPSCL,IOYR,IODAY,TM,RlX,R1Y,R1Z,RlDX,RlDY,RlDZrPlTCH,ROL 

000010 

000 


*L,YAW,PTIMC3) ,TMNC3),TMX C3) ,NLC0EFC2) .NRCOEF C2) ,SCLU0(2) ,SCLL1 C2) , 

000011 

000 


*ELC0EFC11,2),SCLR0C2),SCL'R1C2),ERC0EFC11 ,2) ,NASCEF ,SCLAS0,SCLAS1, 

000012 

000 


*IEL£MN,IELEMX, ASCOEFC1 6) 

000013 

000 


C0HM0N/HINC0H/X,Y,TIHE ,NP 

000.014 

000 


DIMENSION KOUTC24) 

000015 

000 


DATA PI/3.14159265/ 

000016 

000' 


DATA LP/6/ 

000017 

000 


RDPDG=PI/180.0 

000018 

000 


RADL1N=RDPDG*DEGLIN/T0TLIN 

000019 

000 


RADELE=RDPDG*DE6ELE/T0TIEL 

000020 - 

000 


NP=NLHK 

000021 

000 


DO 100 1=1,NLHK 

000022 

000 


ISCAN=1200-CIFIXCALInC I) )-1)/2 

000023 

000 


TIME C1)=PTIHEC1)-U SCAN*TMPSCL 

00002^ 

000 


IDIR=1 

000025 

000 

c 

CONVERT LAT , LON TO ROTATING EARTH CO-OR 

000026 

000 


CALL ERTOER CALA7 Cl ) ,ALOH Cl ),X1 ER , X2 ER , X3 ER , 1 D IR > 

000027 

000 

c 

CONVERT ROTATING EARTH TO INERTIAL EARTH COORDINATES 

000028 

. 000 


CALL ERT0STCX1ER,X2ER,X3ER,X1 ,X2,X3,I0IR,T1H£CI>> 

000029 

000 

c 

EARTH INERTIAL TO SATELLITE LOCAL VERTICAL 

000030 

000 


CALL STTOLV CXi ,X2,X3,I0IR,TIMECD) 

000031 

000 

c 

NOTE : ATTITUDE DATA NOT USED THUS LOCAL VERTICAL SAME AS BODY CENTE 

000032 

000 


X(l',l>=X1 

000033 

000 


XC2,I)=X2 

000034 

000 


XC3,I)=X3 

000035 

000 

100 

CONTINUE 

000036 

000 


DO 200 I=1,NLMK 

000037 

000 


ELEAN6=(P1CELE-AELECI) )*RA0ELE 

00C038 

' 000 


ALNANG=CPICLiN-ALINCI>)*RADLiH 

000039 

000 


Y Cl ,I)=-SINCEL£ANG)*COSCALNANG) 

000040 

000 


YC2,1)=-S1NCALNANG) 

000041 

000 


Y C3,I)=C0S CELEANG) *COS (ALNANG) 

000042 

000 

200 

CONTINUE 

000043 

000 


PRYC1)=0. 

000044 

000 


PRYC2)=0. 

000045 

000 


PRYC3)=0. 

000046 

000 


CALL HINMI2 CPRY ,PRY ,GN ORH , VA LU E, ITER) 

000047 

000 


PITCH=PRY Cl) 

000048 

000 


R0LL=PRY(2) 

000049 

000 


YAH=PRYC3) 

000050 

000 


ENCODE Cl 32, 702 ,MOUT)ITER ,GNORH .VALUE 

000051 

000 

702 

FORKATC1HO,5X, 'CONVERGENCE AT IT ER AT ION' , 16 , / /oX , 'GR AD I ENT NOBM= , 

000052 

000 


*E16.9,' S FINAL V ALU£=' , E16 .9 ) 

000053 

000 


CALL TQCHOUT) 

000054 

000 


RETURN 

000055 

000 


ENDS 


OF POOR QUALUJti 



A6INT/A0IPS 


X«IEGER*4 iriHE 

UlHShSiON KXMUO) ,HH1 nC10) >R1 C3),R2C3] 

CUNj«lQN/,uVCO-4/«A\rN,iNAV,lXR,XOAj:, rOfljlN,OSGIiIN,rOTXblj,Oe.CELE.Pi:Ll 
fN,eiCe‘JE.Xrtf>SCt,XOXR.IUDA5f ,l«,RlX,KlX,RlZ,RlDX.Rlui,RlDZ.PIl’CH.ROL 
»i.,lArt,Pll‘>lCj),lhwtj},Jl'4XC3J, Nl.C0e,FC2J.NRcUEFC2J .SCoLot X) .isCf-Ll C2) . 
4bbCoeFtU,2J ,SCGR0C2) ,SCL«R 1 C 2 ) ,eRC0fc,FCli,2U'^ASCtr,sCl,AS0,SCtiASl, 

♦ XELEM^I , 1 tufci'U . ASCUEF ( 16 J 
data Pl/3.14lii92D535/ 

DAIA OP/6/ 

C PROGftAfl TO XrtilxALlze ATS 6 NAVCOH 

RAD0EG=P1/180.0 
ME«U=1 

call XUCOMC NA(fCOM I’O BE PRINTEO? C1=YES, OEfAULTsKoj • , 42, 1 . 1 , HIM 

♦ ,HEMO,X,U 

XF(^lb^/0.DT.0) GU TO 999 

JPRl=MifiCXJ 

MbHU=l 

CALL XHCOrtC ENTER OAT NUMBER' , 17 , X . 1 , MXN .MENU . 1 . 1 J 

Xf (‘'ENU.GE.O}' GO TO 999 

1«AV=IX11 

TOrLIN=^400.0 

OfcGLXN-19.92 

XUrIEL=2400.0 

DbGEtic,=x0.07 

pxcbi'g=(Tori.XN+x.O)/2.o 

PlCEljE=UOl'X£L+1.0)/2.0 

TMp5Ct<=0.02 

XiR=74 

XUlR=74 

xday=mibcu 

XOOAYsiOAT 

ME«0=3 

CAiit* iNCOMl' ENTER 3 PICTURE START XXMES - HRMMSS ' , 37 . X . I . R^IN . MEN 
♦0,3. n 

XF( MENO.BE.O) go to 999 
xTIhEsRHXWCX) 

PTI'iX U=F1XMECXT1ME) 

XXIME=RpirNC2) 

PlXM(2J=n XMECXTiME) 
inHb=R,UNC3) 

PXIM(3)=f riMECITIHE) 

DO 100 1=1,3 

IMrtCn= 0.0 

1MX(I)=0.0 
lOO CONTINUE 
MbNU=4 

CALL IKCOHC EftlbR FIRST ORoXf POSXTXOn: rrttHMMSS), X. t, Z CKM)». 
♦52,l,l,RHXN,MtN0,3.l) 

IF (nbWU.Ob.O) GU TO 999 

irxME=RniNCu 

ii=b riMburxab) 

RKl J=Ri'iXiNC2) 

RX(2X=RMXN( 3) 



A6INT/A0IPS 


Hlt3J=ftMIWC4J 

ckIjL iwcohc E.fJreR stcowo oftuir position: rcHH«HSs), x. x, z ckhs' 

♦■,53.1,l,KrtlN, MENU , 3 , 1 J 
iF(MEftU.tiE.O) (.0 ro 999 
ITIME^RMINU) 

T2=f riMECIXlMED 
K2(U=KMIN{2) 

P^C2isKM£N(3) 
ftz{3J=KMlNC4i 
CAliL GAS0RB(RI,n,R2,T2J 
99 b CUNIINUE 

1ft jPRi.uE.O) GO TO 1000 
rtRIlECLP, 7f)l)NAVN .INAV ,llR,IDAr 

761 fORMATClHl, 'NAV«='fi5,T27,'lNAV='»Il0.r5 2,« xEAR=* . 12 . 177 • • DAX = ' . 13 
♦J 

«RlTE(bP,7b21 roTolN.OEGLlN.TOTlEO.DEGEOE 

762 FORMAi’t'fa. ‘rO‘Ibiri=' ,f 10.2.T27, ' OeGlIns ' ,F10. 2 . T52 . 'IOXlbb=* .F10.2, 
+T77. 'DEGeOfcs' .PlO.2) 

rfKllECL,P.7o3)PiCL.iN.PlCELE, IMPSCl 

763 fORHAlCf2, 'flCbiNs' ,‘’'10.2,127. ' PlCEO^s • , Flo.2 . T52 . 'laPSCOs' .F10.2) . 
HRilEtuP. /b4UOIR,IOOAI 

764 FURMAltri. 'OKbif X E AK= ' , 1 4, 127 , • ORBIT OAX=* .14) 
rtRllE(bP,76S) IM 

7b5' F0RMATtr2, 'ORBIT TIMt, TM=»,F15.4) 

WRIlEtbP,7b6JRXX,RlX.Rll 

766 F0kmAi(I 2, 'POSITION CK.a) ' . T27 , 'RIX=' , FI 5. 3 .152, ' R 1 ' ,F15 . 4, 1?7 , ' R 
41Z=> ,FlS.4) 

WRITE (LP, 767 )R1DX,R1DX.R1D2 

767 f 0RaATtl2, ' VbLOCiTX ' , 127, ' kl0X= ' , 6l5 . 9 . T52 . ' R10X=' . &15. 9 , T77 , ' RIOZ 
4=’.E15.9) 

PDsRlTCH/RAOOEG 

R0=80bb/RADD6G 

UisiA'^/i^ADOliG 

nRITE(bP,768)P0,RO,X0 

76a tORMAm,!. 'ATUfUOE (DEG) ' ,T27, 'PIICH=' ,£l5.9,r52. ' ROLb= ' , E15 .9 . T7 
*1 , • XA/{=' ,El5.9) 
rtRlT£(i,P,769)RTlM 

769 tORMATtra, 'PICTURE TIMES • , 1 27 . 3F2S . 4) 

770 FOrtMATax,5l25) 

771 F0RHAXC1X,5E25.9) 

WRITE CbP, / 70)NbCQe,F,Nkc0bF,«A6CEF 
WRirE{bP,77l)IMN.XHX 
WRITE (bP, 770lNbCJEf .NRCOEF 
bRII'ECI,P, 7713 SCbbO.sCLbl 
rtRITEtbP.77l)EbCOEf 
ttHilE(Up,771) SCORO.SCOrI 
ttRn£(uP,771J ERCOEF 
ftRlT£(i.P,770)NASCEf .lEOEMK.IEbEMX 
wRlTEtbP,771) SCbASO.aCbASI 
.milE(bP.77l) A5C0EF 
1000 CONTINUE 

CALL ReQU£S(RADS 0('ATSF2»)) 
eN|> 




ORIGINAL PAGE IS 
OF POOR QUALITY 



BCTOP F /HC lOA $ ****** 


aELT,L AF.BCTOPF/fICIDAS 

ELT007 RLIB62 1 2/ 2 2-1 7: 0 1 : 42 - ( 0, ) 


000001 

000 


SUBROUTINE BCTOPF ( X , Y , 2 , lOIR ) 

000002 

000 


DIMENSION At3,3> 

000003 

000 


COMMOK/NAVCOH/NAVN, INAV.lYRt 1I>AY jTOTLIN tOEGLIN.TOTlEL, DEGELE,PICLI 

000004 

000 


*N.PICELEf TKPSCLf lOYR, lODAYtTM, R1 X.RIY tR12 ,R1 OX, RlDVtRlDZ, PITCH, ROL 

000005 

000 


*L,YAW,PTIH(3>,TMN(3),TMX(3),NLCOEF(2) ,NRCO£F (2) rSCLLOCZ) ,SCLU (2) » 

000006 , 

000 


*£LC0EF(11 ,2> ,SCLR0(2) ,SCLR1<2) , ERCO E F<11 , 2) , N ASCEF , SCL ASO, S Ct-ASl , 

000007 

000 


*IEL£HN,1ELEHX,ASC0EF(16) 

000008 

000 

c 

PT(3> = YAW, PTC2) = ROLL, PT(1) = PITCH 

000009 

000 


CALL UHITCA) 

000010 

000 


CALL ROTATE CA, YAU, 3,1 ) 

000011 

000 


CALL ROTATECA, ROLL, 1,1) 

000012 

000 


CALL R0TATE(A,PITCH,2,1) 

000013 

000 


IFaDXR.E0,2)60 TO 10 

000014 

000 


XT = X*A(1 ,1 )+Y*A(1 ,2)+2*A(1 ,3> 

000015 

000 


YT = X*A(2,1)+Y*AC2,2)+Z*At2,3) 

000016 

000 


2T = X*A(3,1)+Y*A(3,2)+Z*AC3,3) 

000017 

000 


X=XT 

000018 

000 


Y=YT 

000019 

000 


Z=ZT 

000020 

000 


RETURN 

000021 

000 

10 

XT=X*A(1 ,1)+Y*AC2,1)+Z*A(3,1) 

000022 

000 


YT = X*A(1 ,2)+Y*A(2,2)+Z*AC3,2) 

000023 

000 


Z = X*A(1 ,3>+Y*AC2,3)+Z*A(3,3> 

000024 

000 


X=XT 

000025 

000 


Y=YT 

000026 

000 


RETURN 

000027 

000 


ENDS 


EMD ELT. 

3HDG,P ****** CNPS/IBHSSP ****** 



****** BCTOPF/MCIDA$ ****** 


3ELT,L 
ELT007 
000001 
000002 
000003 
000004 
000005 
000006 
000007 
OOOOOS 
000009 
000010 
000011 
000012 
000013 
000014 ’ 
000015 
000016 
000017 
000018 
000019 
000020 
000021 
000022 
000023 
000024 
000025 
000026 
000027 

END ELT 
aHD6,P 


AF.BCTOPE/rtCID AS 


RLIB62 1 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
‘ooo 
000 
000 
000 


2/22-17:01:42-(0t> 

SUBROUTINE 8CT OP F C X ♦ Y , 2 , IDIR > 

DIMENSION A'C3,3) 

’ COMMON/NAVCOK/NAVNtINAV ,IYRt IDAY ,TOTLIN ,0E6LIN,T0TIELr DEGELEjPICLI 
*NfPlCELE,TMPSCLtIOYRrI0D AY tTKtRlX,RlY,Rl2,R1DX,RlDY,RlDZ, PITCH, ROL' 
*L ,YAW ,PTIM C3) ,TMN (3) ,TMX C3> rNLCOEF (2 ) tNRCOEF C2) r SCLLO( 2) tSCLLi C2> r 
*CLCO£F C1 1f2>,SCLR0(2)«SCLR1(2),ERC0EFf11,2) ♦ NAS C £F , SCL AS 0» S CLASi » 
*IELEMN,IEL£MX, ASCOEFCI 6) 

C PT(3) = YAW, PTC2) = ROLL, PT(1> = PITCH 

CALL UNIT(A) 

CALL ROTATECA, YAW,3,1) 

CALL ROTATE (A , ROLL ,1,1 > 

CALL ROTATECA, PITCH, 2,1) 

IFC1DIR.EQ.2)60 TO 10 

XT - X*A(1 ,1)+Y*A<1,2>+2*A(1 ,3) 

YT - X*A<2,1 )+Y*A(2,2)+Z*A(2,3) 

2T = X*AC3,1)+Y*A(3,2)+Z*AC3,3) 

X=XT 

Y=YT ^ 

Z^ZT 

return 

10 XT=X*A(1,1)+Y*AC2,1)+2*A(3,1) 

YT X*AC1 ,2) + Y*A(2,2)+Z*A(3,2) 

2 = X*A(1 ,3)+Y*A(2,3)+Z*A<3,3> 

X=XT 

Y=YT 

return 

ENDS 


****** CNPS/IBHSSP ****** 





***^*ik UNIT/MCIDAS ****«<r 


3ELT,L AFeUNIT/MCIDAS 


ELT007 

RL1B62 12/22- 

000001 

000 

OOOOD2 

000 C 

000003 

000 C 

000004 

000 

000005 

000 

000006 

000 

000007 

000 

000008 

000 

000009 

000 

000010 

000 

000011 

QDO 


1 7:02:09-C0t) 

SUBROUTINE UNIT(A) 
THIS ROUTINE RETURNS IN 

DIMENSION A(9) ,BC9) 
DATA 0/1 .OrO.OtO.Q, 

^ 0»0y1«0f0»09 

* 0«0f0«0f1o0/ 

DO 10 1=1 ,9 
10 A(1)=BCI> 

RETURN 

ENDS 


A 3 X 3 IDENTITY MATRIX^ 


END ELT, 

BKD&fN 

SPIN 



****** STTOLV/HCrDAS ****** 


3ELT»L AF.STTOLV/MCID 
ELT007 RLI862 12/22-1 
000001 000 

000002 000 C 

000003 000 C 

000004 000 C 

000005 000 C 

000006 000 

000007 000 

000008 000 

000009 000 

000010 000 

000011 000 

000012 000 

000013, 000 

000014 000 

000015 000 

000016 000- 

000017 000 10 

000018 000 

000019 000 

000020 000 

000021 000 


AS 

7:02:08-(0»> 

SUBROUTINE STTOLV OC,V 
IF IDIR=1, POINTING VE 
FROM SrtT INERTIAL TO U 
IF I0IR=2, POINTING VE 
LOCAL 'VERTICAL TO SAT 
CALL ORBITCXS, YS,Z StTI 
CALL NRMLIZCXStVS.ZSfX 
X1=X 
V1=Y 
Z1=Z 

D=SQRT(XS**2+TS**2 J 

IF (ID1R.EQ.2) GO TO 1 

X*C-YS*X1+XS*Y1)/D 

Y=CXS*2S*X1+YS*ZS*Y1-Z 

Z=-CXS*X1+YS*Y1+ZS*Z1) 

RETURN 

X*-YS*X1/D+XS*ZS*Y1/P- 

y*XS*X1/D+YS*ZS*Yl/D-Y 

Z*-D*Y1-ZS*Z1 

RETURN 

ENtl$ 


>Z ,161 
CTOR ( 
OCAL V 
CTOR ( 
INERTI 
fie) 
NORM) 


0 

1 * D **2 


XS*Z1 

S*Z1 


R,TIHE> 

X,Y,Z) IS TRANSFORHEB 
ERTICAL FRAME. 

X,Y,Z) IS TRANSFORMED 
AL FRAME. 


)/D 


fROPI 


END ELT. 

3HD6,P ****** UNIT/HCIOAS ****** 



****** SE/MCIDAS 

c)ELT,L AF^SE/HCIDAS 

ELT007 RLIB62 1 2/ 2 2-1 7: 02 : 06-C0, ) 


000001 

000 


SUBROUTINE SE < PTIHE ,XL IN ,XEL E , ALAT, ALOW ) “ 

000002 

OOO 


COKHON/NAVCOM/NAVN fINAVt ITR, IDAY,T0TLIN,DE6L1N,T0TIEL, DEGELEtPICUl 

000003 

000 


*N ,PICELE,THPSC1.,I0 YR,10t)AY,TM,R1X ,R1 Y,R12 ,R1 DX,R1DY, RT DZ f PI TCH , ROL 

OOOOOA 

000 


*L,YAW,PTIM(3),TMNC3>,THX(3>,NLC0EF(2),NRC0EF<2) t SCLL0<2) fSCLtI (2) , 

000005 

000 


*ELC0EF(11 ,2) ,SCLR0(2J ,SCLR1 (2) ,ER COEF Cl 1 ,2) , NASC EF , SCLASO.S CLAS1 , 

000006 

000 


♦lELEHN jIELEHX, ASC0EFC1 6) 

000007 

OOO 


DATA lMT/0/ 

ooooos 

000 


CALL EDGCOR CPTIME fXLINtD ELLIN, DELELE) 

000009 

OOO 


ALIN=XLIN-DELLIN 

000010 

OOO ’ 


AELE=XtL£-DELELE 

000011 

OOO 


ISCAN=1200-<IFIX CXLlN-1) )/2 

000012 

OOO 


TIHE=PTIME+ISCAN*TKPSCL 

000013 

OOO 


I0IR=2 

000014 

OOO 


CALL PFT0TCCALIN,AELE,X1 ,X2, X3 ,IDIR ,INIT) 

00001 5 

OOO 


CALL BCrOPFCXI ,X2,X3,IDIR) 

000016 

OOO 

c 

BC SAME AS LV HERE. 

000017 

OOO 


CALL STTOLVCXi ,X2,X3,1DIR,TIME) 

000018 

OOO 


CALL ERTOST(X1ER,X2ERtX3ER,X1,X2,X3,1DIR,T1ME) 

000019 

OOO 


CALL ERTOER (ALAT , AL0N,X1 ER,X2ER, X3ER,I0IR> - 

000020 

OOO 


RETURN 

000021 

OOO 


END 


END ELI* 

aHDGrP ****** STTOL V/NCIDAS ****** 


-OF POOR QUALtl^.. 



****** SATEAR/MCIDAS ****** 


SELT,L AF.SATEAR/HCID AS 

ELT007 RLIB62 1 2/H2-1 7:02 :05-(0, ) 


000001 

000 


SUBROUTINE SATEA R CPICTIH ,XLI N , XELE , XLAX, XLON , ITY PE , INA V, BETAIH ,BET 

000002 

000 


*OOT,ATfRAC) 

000003 

000 


PT1ME=PICTIK 

000004 

000 


60 TO (1,2,3,4,S),1TYPE 

000005 

000 

1 

CALL SE(PTIHE,XLIN,XELE,XLAT,XLON> . 

000006 

000 


RETURN 

000007 

000 

2 

CALL ESCPTIHE rXLAT,XLON,XLIN fXELE) ■ 

000008 

000 


return 

000009 

000 

5 

CONTINUE 

000010. 

000 


RETURN 

000011 

000 

4 

CONTINUE 

■000012 

ooo 


■ RETURN 

000013 

000 

5 

CONTINUE 

000014 • 

000 


RETURN 

0000 1'5 • 

000 


END 


END ELT, 

3HDG,P ****** SE/HCIDAS ****** 


tolGlNAL PAGE TQ 
OP POOR QVALm 



**i*it* S/HCIDAS ****** 


aELT,L AF. S/HCIDAS 


ELT007 

RL1B62 12/22-1 

7s02 :0A-(0,) 

000001 

000 


FUNCTION SCPT) 

000002 

ODD 

c 


000003 

000 

c 

THIS FUNCTION RETURNS AS ITS VALUE THE VALUE OF THE OBJECTIVE 

000004 

000 

c 

FUNCTION S AT THE POINT "PT“, 

000005 

000 

c 


000006 

000 


DIMENSION XC3,150>tY(3,150>fTIHE<150> 

000007 

000 


COMMON/MINCOM/X, Y.TIMEtHP 

000008 

000 


DIMENSION PT(3),H(3»3) 

000009 

000 


sres»o 

OO0010 

000 


CALL UNiTCH) 

000011 

000 


CALL R0TATE(HtPT(3)f3,1) 

000012 

000 


CALL R0TATE(H,PT(2),1,1) 

000013 

ooo 


CALL R0TATECh,PTC1) »2,1) 

000014 

000 


DO 10 J=1,NP 

000015 

000. 


DO 10 1=1,3 

OOOQ16 

000 


10 SRES=SRES+(Ya,J>-H(I,1)*X(1 , J J -H Cl , 2) *X C2 , J ) -H Cl, 3) *X C3 , J ) )**2 

000017 

000 


S=SRES 

000018 

000 


RETURN 

000019 

000 


ENDS 


END ELT. 

3HDG»P ****** SATEAR/MCIDAS ****** 



****** ROTATE /«C IDA S ****** 


SELT,L 

ELT007 

000001 

000002 

000003 

000004 

000005 

000006 

000007 

000008 

000009 

000010 

000011 

000012 

000013 

000014 

000015 

000016 

000017 

000018 

000019 

000020 

000021 

000022 

000023 

000024 

000025 

000026 

000027 

000028 

000029 

000030 

000031 

000032 

000033 

000034 

000035 

000036 

G00037 

000035 

000039 

000040 

000041 

END £LT 
aHDGfP 


AF*ROTATE/MCID 
RLIB62 12/22-1 
000 

000 C 

000 C T 

000 C ” 

000 C R 

000 C I 

000 C I 

000 C 

000 
000 
ODD 
000 
000 
000 
000 
GOO 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
ooc 
000 
000 
000 
000 
000 
000 
000 
000 


AS ' 

7:02:03-(0t) 

SUBROUTINE ROT AT E C A ,R , IR , I DE RI V) 

HIS ROUTINE RETURNS IN THE PRODUCT OF THE INPUT MATRIX 

A” AND A MATRIX RM^ WHERE r IF '*ID£RIV” = 1, RM REPRESENTS A 
OTATION THROUGH AN ANGLE CIN RADIANS) ABOUT THE AXIS 

DERIV=2t the first derivative OF RM IS OPERATED ON Ay AND IF 
D£RIV^3, THE SECOND DERIVATIVE OF RM IS USED. 

DIMENSION A(3,3>tlNDX1 C3)rINDX2(3) 

DATA INDX1,INDX2/2,1t1 j3r3,2/ 

IR1=1NDX1C1R> 

IR2-INDX2CIR) 

CR=DCOS(R> 

SR=DSINCR) 

IFClDERIV .NE.1 )G0 TO 2 
DO 1 J=1,3 
T1=AC1R1,J) 

T2=A(IR2, J) 

AClkl ,J)=CR*T1+SR*T2 
ACIR2, J) = -SR*T1’^CR*T2 

1 CONTINUE 
RETURN 

2 IFCIDERIV.NE.2)G0 TO 4 
DO 3 J=1,3 

A(IR, J)^0.0 
T1=ACIR1, J) 

T2=ACIR2,J> 

A (IRI tJ) = -SR*T1 + CR*T2 
A CIR2 ,J)=-CR*T1-SR*T2 

3 CONTINUE 
RETURN 

4 CONTINUE 
DO 5 J , 3 
A(IR, J)=0.0 
Tl=A(IRlrJ) 

T2 = A (IR2tJ ) 

A (IRI ,J)=-CR*T1-SR*T2 
A(IR2, J)=SR*T1-CR*T2 

5 CONTINUE 
RETURN 
ENDS 


If 


****** s/mcidas ****** 





■ieis ifki'it 

PRTIAL/HCIOA S 


000056 

000 


CALL ROTATE(HBC,B,1 ,2) 

000057 

000 


CALL ROTATE(HBC,A,2,1) 

000058 

000 


Call unit(hcc) 

000059 

000 


CALL ROTATECHCCfCf 3t3> 

000060 

000 


CALL ROTATECHCC, 8,1,1) 

000061 

000 


CALL ROTATE(HCC,A,2,1) 

000062 

000 


D010J=1 ,NP 

000063 

000 


00101=1,3 

000064 

000 


T=Y<l,J>-fUI,1)*X<1,J)-Ha,2)*X(2,J)-H(I,3)*X(3,J) 

000065 

000 


TA=HA(I,1)*X(1,J)4-KA(I,2)*X(2,J) + HA{I,3>*X{3,J) 

000066 

000 


TB=HB(1,1)*XC1,J)+KB(I ,2>*X(2,J)+HBU,3)*X<3,J) 

000067 

000 


TC = KC(I,1)*X(1,J)+HC(I,2)*XC2,J)+HCU,3)*X{3,J) 

000068 

000 


TAA = KAAa,1)*XC1,J)+HAA(I,2>*X (2 , J ) +HAA < 1 ,3 ) *X (3 ,4 ) 

000069 

000 


TAB = HAb(I ,1)*X (1 ,4 >+HA8(I,2)*X (2 , J ) +HA8 ( I ,3 ) *X t3 ,4 > 

000070 

000 


TAC=HAC(I,1)*X(1,4)+HAC(I,2>*X(2,4)+HAC(I,3)*X(3,4) 

000071 

000 


T0B = HBea,1)*X(1,J) + H3B(I,2>*X (2 , 4 ) +HBB ( 1 ,3) *XC3,4 > 

000072 

000 


TBC=KBC(1,1)*X(1,4)+HBC(I,2)*X(2,4)+HBC{I,3)*X(3,4) 

000073 

000 


TCC = HCCa,1)*X(1,4)+HCC(l,2>*X<2,4) + HCC(l,3)*X(3,4> 

000074 

000 


PSPA=PSPA-T*TA 

000075 

000 


PSPB=PSPB-T*TB 

000076 

000 


PSPC=PSPC“T*TC 

000077 

000 


PSPASQ=PSPASQ-T*TAA+TA**2 

000078 

000 


PSPAPB=PSPAPB-7*TAg+TA*T8 

000079 

000 


PSPAPC=PSPAPC-T*TAC+TA*TC 

000080 

000 


PSP 8 SQ=PSP 0 SQ-T*TBe+TS **2 

000081 

000 


PSPbPC=PSPBPC-T*T8C+TB*rC 

000082 

000 

10 

PSPCSQ=PSPCSG-T*TCC+TC**2 

000083 

000 


GRAD < 1 )=2*PSPA 

000084 

000 


GRAD C2)=2*PSPB 

000065 

000 


GRAD C3) = 2*PSPC 

000086 

000 


HESS(1,1)=2*PSPASO 

00C0S7 

000 


HESS (1 ,2) =2*PSPAP8 

000088 

000 


HESS C2,1 )=2*PSPAPB 

000039 

OQO 


HESS (1 ,3>=2*PSPAPC 

000090 

000 


HESS(3,1)=2*PSPAPC 

000091 

000 


H£SS(2,2>=2*PSPSSQ 

000092 

000 


HESS (2,3)=2*PSPBPC 

000093 

000 


HESS (3,2)=2*PSPBPC 

000094 

000 


HESS C3,3>=2*PSPCSQ 

000095 

000 


RETURH 

000096 

000 


ENDS 

END ELT 




SHDG,P 

****** 

ROTATE/KCIDAS ****** 



**Air4* PRTI Al / f1 C I D A S ****** 


SELTft AF.PRTIAL/MCID AS 

ELT007 RLIB62 1 2/ 2 2-1 7:02 :03-(0 » ) 


000001 

000 


SUBROUTINE PRTl AL (PT, GRAD* HESS ) 

000002 

000 

c 


000003 

000 

c 

THIS ROUTINE RETURNS IN "GRAO” THE GRADIENT OF S AT “PT" AND IK 

000004 

000 

c 

"HESS" THE HESSIAN OF S AT "PT". 

000005 

000 

c 


000006 

000 


DIMENSION XC3, ISO) ,Y(3, 150) , TIME (ISO) 

000007 

000 


COMHON/HINCOH/X, YtTIME ,HP 

000008 

000 


DIHENS10NPT<3),GRADC3),HESS(3,3),H{3,3),KA{3,3),HB<3,3),HCC3,3), 

000009 

000 


*HAAC3,3),HAB(3,3) , HAC ( 3 , 3) , HBB <3, 3) , HBC< 3 , 3 ) , HtC t3 ,3) 

000010 

000 


A=P7(1) 

000011 

000 


B=PT(2) 

000012 

000 


C=PT(3) 

000013 

000 


PSPA=0 

000014 

000 


PSP3=0 

000015 

000 


PSPC=0 

000016 

000. 


PSPASQ=0 

000017 

000 


PSPAP8=0 

000018 

000 


PSPAPC=0 

000019 

000 


PSPBS9=0 

000020 

000 


PSPBPC=0 

000021 

000 


PSPCSQ=0 

000022 

000 


CALL UNITCH) 

000023 

000 


CALL R0TATECH,C,3,1) 

000024 

000 


CALL R0TATE(H,B,1 ,1) 

000025 

000 


CALL ROTATECH, A,2t1> 

000026 

odo 


CALL UNIT(HA) 

000027 

000 


CALL R0TATE(HA,C,3,1) 

000028 

000 


CALL RCTATE(HA,B,1 ,1) 

000029 

000 


CALL R0TATE(HA,A,2 ,2) 

000030 

000 


CALL UNITCHB) 

000031 

000 


CALL ROTATECHB ,C,3 ,1) 

000032 

000 


CALL R0TATE(HB,8,1 ,2) 

000033 

000 


-CALL ROTATECHB ,A, 2,1) 

000034 

000 


CALL UNITCHC) 

000035 

000 


CALL R0TATE(HC,C,3,2) 

000036 

000 


CALL R0TA7ECHC,B,1 ,1) 

-000037 

000 


CALL ROTATEChC,A,2,1) 

000038 

000 


CALL UNIT(HAA) 

000059 

000 


CALL R0TATECHAA,C,3,1) 

000040 

000 


CALL ROTATE (HAA,b, 1 ,1 ) 

000041 

000 


CALL R0TATE(HAA,A,2,3) 

000042 

000 


CALL UMTCHAB) 

000043 

000 


CALL R0TATE(HAB,C,3,1) 

000044 

000 


CALL R0TATEChAB,B,1,2) 

000045 

000 


CALL RCTATECHAb,A,2,2) 

000046 

000 


.CALL UNIT(HAC) 

000047 

000 


CALL ROTATECHAC,C,3,2) 

000048 

000 


CALL R0TATE(HAC,B,1 ,1) 

000049 

000 


CALL ROTATE(HAC,A,2,2) 

000050 

000 


CALL UNIT(HBB) 

000051 

000 


CALL ROTATE,<HBB,C,3,1) 

0000S2 

000 


CALL ROTATE (HBb,b,1 ,3) 

000053 

000 


CALL ROTATE(HBB,A,2,1) 

000054 

ood 


■ CALL UNIT(HBC) 

000055 

000 


CALL ROTAT£(HBC,C,3,2) „ 







****** PFTOTC/MCIDAS ****** 


Af .PFTOTC/MCID AS 
RLIB62 12/22-17; 
000 
000 
000 
000 
000 
000 

000 C 
000 c 
000 
000 
000 
000 
000 
000 

000 1 
000 ’ 

000 

000 

000 

ooo 


02r02-t0,> 

SUBROUTINE PFTOT C CXLI N ,X EUE , X , Y, Z »IoiR , i«iT ) 

COKMON/NAVCOH/NAVN,lNAV,lYR,IOAY,TOTLIN,DEGLIN,TOTlEL,OEGELE,PlCt.I 

*L,YAV,,PTIW(3) ,T«N{3),TKX C3 1 , NL COE F (2> tKRC0EF(2) , S CLLO ( 2) , SCLL 1 C2> , 
*ELCOEF(11,2),SCLRO(2>»SCLRU2),ERCOEF(11,2) , NASC E F , SCLASOt S CLAS1 , 
*1ELE«N,I£LEMX,ASC0EF(16> . 

IF IDIR = 1, X,Y,Z 70 EIN, ELE 
IF IDIft = 2, LIH, ELE TO X,Y,2 
DATA PI, RE, GR A CON /3 ♦14159265 ,6378 . 1 5 07436574/ 

IF(1NIT.EQ.2)G0 TO 1 

INIT = 2 

RDPOG s: PI/180. 

RADLIN = ROPDG*OEGLIN/TOTLIN 
RADELE = RDPDG*0£6ELE/T0TIEL 
if(ipir.eq.2>go to 10 

ELEANG = ATAH(X/Z) 

XLNANG = ASIN(V) 

XELE=PICELE+ELEANG/RADELE 
XLlh = PICLIN+XLHANG/RAOLIN 
RETURN 

£LEANG=(PICELE-XELEJ*RADELE 
XLNANG = (PICL1N-XLIN)*RADL1N 
X=-C0$(XLNAN6)*S1N (ELEANG) 

Y=-S1NCXLNAN6) 

Z=C0S (XLNANG)* COS (ELEANG) 

RETURN 

ENOS 


END eLt. ■ 

3Ht>G,P ****** PRTIAL/HCIDAS ****** 



ORSXT/MCIOAS 


aELT,t 

AF .ORBIT/MCIDA S 

ELT007 

RIIB62 12/22“ 

1 7:02;01-(0,> 

000001 

000 

SUBROUTINE ORB IT C X , Y, 2 »T ) 

000002 

000 

REAL KE,HU,H«INMO 

000003 

000 

COHHON/NAVCOK/NAVNf INAV,IYR, JO AY,T0TLIN,DE6L1N,T0TIEL»0EGELE,PICLJ 

000004 

OOG 

*N,PlCELE,TKPSCL,IOYR,100AY,THrRlX»RlY,ft1Z,Rl DX,R10Y,R1 OZ , PITCH , ROt 

000005 

000 

*L,YAWtPTIHt3),THN<3),TMX ( 3 ) , NLCOE F( 2 J , NRC OE F (2) , SCLLO( 2 > ,S C LL 1 C2)» 

000006 

000 

*ELC0EF(11 ,2),SCLR0(2) ,SCLR1 C2) ,£R COE F < 1 1 , 2) t NAS CEF , SCL ASO t S CLA Si f 

000007 

000 

*iELEMNtIELEnx,ASC0EF(16> 

OOOOOS 

000 

DATA Rc,KE,KU, EPSlLN/6 378. 15 fO .07436574 , 1 .0 , 1 .OE-7/ 

000009 

000 

SQRTKU=SQRT<«U) 

000010 

000 

R0=SQRT(R1X**2+R1 Y**2+R1Z**2 J 

000011 

000 

O0=(RlX*RlDX+RlYARlDY+RlZ*R1 DZJ/SQRTMU 

000012 

000 

V0SU«U = (R1 DX<r*2 + RlOY**2 + Rl DZ **2) / SQ RTHU 

000013 

000 

A=RO/ C2,-R0*V0SQHU) 

000014 

000 

CE=( A-ROJ/A 

000015 

000 

SE=DO/SaRTCA> 

000016 

000 

HKIN«0=KE*(T-TH)»SQRTHU/A**1 .5 

000017 

000 

GL=.5*HHI«M0 

000018 

000 

5 S1NGL=SIN(6L) 

000019 

000 

SNCSGL=SINGL*COS (GL) 

000020 

000 

XNUM=6l+SE*SINGL**2-CE*SNCSGL-.5*HHINR0 

000021 

000 

D£N0R = 1 + 2*SE*SKCSGL-CE*( 1-2* SI NGL*-* 2) 

000022 

000 

G=CGL*DENOR-XNUMJ /DENOn 

000023 

000 

IF<ABS(GL-G) .LT.EPSILNJGOTOTO 

000024 

000 

GL=G 

000025 

- 000 

GOTO 5 

000026 

000 

10 ECAN0M=2*6 

000027 

000 

C=A* (1-COS (ECANOM) > 

000028 

GOO 

S = SQRT(A J *SIN(ECANOH) 

000029 

000 

F=(RO-C) /RO 

000030 

000 

6=1/SQRTHU*(R0*S*D0*C) 

000031 

OOG 

X=(F*R1X+6*R10X>*RE 

000032 

OOG 

Y=( F*R1Y+G*R1DY)*RE 

000033 

OOG 

Z=(f *R1Z+G*fi1DZ>*RE 

000034 

000 

RETURN 

000035 

000 

ENDS 

END ElT 

* 


3HD6,P 

****** PACK 

****** 




****** OFSTG/MCIDAS ****** 


000113 

OOO 

c 

eUMP TABLE VALUES. 

0001U 

000 


MC0L=(NC0R-1)/5 + 1 

000115 

000 


J0UT=2 

000116 

000 


DO 650 IWH=1,NIC0L 

000117 

000 


ENC0DEC132,710,M0UTXIC00R(I«) .OFFSET ( IW ) ,W ATE CI U) , I b = IWWtNCOR 

000118 

000 


♦NICOL) 

000119 

000 


CALL TP(J0UT,t10UT5 

000120 

QOO 

650 

CONTINUE 

000121 

000 

710 

FORMAT ((5C5X,14,1X,F8.3,FS.4)>> 

000122 

000 


NCORKX=NCOR 

000123 

000 


CALL OYNASGC'OFFSTO'.LUNF) 

000124 

000 


CALL OPN(LUNF) 

000125 

000 


CALL WRITW(LUNF,1151,NC0RHX> 

000126 

000 


RETURN 

000127 

000 

990 

CONTINUE 

000128 

000 


CALL TQC72H OFFSETGEN ERROR RETURN. 

000129 

000 


* 

D00130 

000 


RETURN 

000131 

000 


ENDS 

000132 

000 

$F1 LEHA 

000133 

ooc 

DELETE GCCGCC,G0RP 

000134 

000 

SINCLUDE IFLD 

000135 

000 

SCATAL06 

000136 

000 

NAH £= 

GCCGCC r5,Rf WtD 

000137 

000 

TYPE= 


000138 

000 

LIB =ATSFLB ,LL 

000139 

000 

BEGIN 


OOOUO 

000 

SEO J 



END ELT. 

QHDGtP ****** OFSTVL ****** 



0FST6/MC1DAS ****** 


000056 

000 


CALL BSRCLUN) 

000057 

000 


CALL BSRCLUN) 

000058 

ODD 

c 

INPUT PAIR OF SCANS. 

000059' 

000 


CALL REAOWCLUNtMwStieOATA) 

OOOOoO 

000 


CALL RLADW(LUN,NW2,I00ATA) 

000061 

000 


CALL CRKTHR(195|1E0ATA,IEBRT) 

000062 

000 


JESCAN=1£BRTC195) 

000063 

000 


CALL CRKATS(NBP,IEOATA,IEBRT) 

000064 

000 


CALL CRKTHR (195,10DATA ,I08RT) 

000065 

000 


J0SCAN=I0BRT(1 95) 

000066 

000 


CALL CRKATSCNBP,10DATA,IOBRT) 

000067 

000 


IFCJESCAKI.NE.IESCAN.OR.JOSCAN.NE.IOSCAN) 60 TO 990 

000068 

000 

c 


000069 

000 

c 

60 correlations. 

000070 

000 

c 


000071 

000 


DO 600 KCOR^ItNCOR 

000072 

000 


IELE=I£LEST+ (KC0R-1)*IELINT 

000075 

000’ 


1COORCKCOR>=IELE 

000074 

000 


HAXLA6=(NLAG-1)/2 

000075 

000 


mNLA6=-HAXLAG 

000076 

000 


DO 500 KLA6=PHNLA6 ,MAXLAG 

000077 

000 


KAXBRT=0 

000078 

000 


«INBRT=5T2 

000079 

000 


NTSIZ=CNTSIZ-l)/2*2+1 

000080 

000 


FLP=0 

000081 

000 


DO 400 KVAL=1,NTSIZ 

000082 

000 


ICT=IELE+KVAL-1-NTSIZ/2 

000083 

000 


1 CS = I ELE + KVAL-I+KLAGi-NOMOFF-NTSIZ/2 

000084 

000 


ICTX = ICT + IDOF.F 

000085 

000 


ICSX=ICS+100FF 

000086 

000 


1PIXT=I08RT(ICTX) 

000087 

000 


iPlXS=lEBRT(ICSX)' 

000088 

000 


HAXt>RT=MAXO(MAXBRT,IPIXT) 

OOOOS9 

000 


HINBRT=MI(vOCKINBRT,IPIXT) 

000090 

000 


FLP=FLP+CIPIXT-IPIXS)**2 

000091 

000 

400 

CONTINUE 

000092 

000 

c 

STORE LPC2) MEASURE VALUE 

000093 

000 


FLPARY(KLA6-MINLAG+1)=FLP 

000094 

000 

500 

CONTINUE 

000095 

000 

c 

SEARCH TABLE OF VALUES OF LP(2> MEASURE FOR HIN. 

000096 

000 


10FF=N0H0FF+HINLAG 

000097 

000 


FLPMIN=FLPARY(1) 

000098 

000 


DO 510 K=2,NLAG 

000099 

000 


I F CFLPHIH.6T.FLPARYCK) ) lOF F =N0«0 FF ♦ (K-1 ) +H INLAG 

000100 

000 


IFCFLPMIN .GT .FLPARY (X) ) FLPH IN =FLPAR Y CK) 

000101 

000 

510 

CONTINUE 

000102 

000 

c 

WEIGHT OFFSETS BY BRIGHTNESS RANGE. 

000103 

000 


WT = FLOAT (MAX6RT-iniNERT)/512. 

000104 

000 


KATE (K COR) =W ATE (K COR) +WT 

000105 

000 


0FFSETCKC0R)=0FFSET(KC0R)+10FF*«T 

000106 

000 

600 

CONTINUE 

000107 

000 


IF(.NOT.ISLAST) GO TO 100 

000108 

000 

c 

COMPUTE AVARA6E OFFSETS. 

000109 

000 


DO 610 KCOR=1,NCOR 

000110 

000 


. lF(wATt(KCOR).LT.1 .OE-3) GO TO 610 

000111 

000 


OFFSET(KCOR)=OFFSET (KCOR)/HATE CKCOR) 

000112 

000 

610 

CONTINUE 



ilr * 0 T 6 / I D A S ****** 


SELT,t 

ELT007 

000001 

000002 

000003 

000004 

000005 

000006 

000007 

000008 

000009 

000010 

000011 

000012 

000013 

000014 

000015 

000016 

000017 

000018 

000019 

000020 

000021 

000022 

000023 

000024 

000025 

000026 

000027 

000028 

000029 

000030 

000031 

000032 

0D0033 

000034 

000035 

000036 

000037 

000038 

000039 

000040 

000041 

000042 

000043 

000044 

000045 

000046 

000047 

000048 

000049 

000050 

000051 

000052 

000053 

000054 

000055 


AF*0FST6/«CIDA 
RLIB62 12/22-1 
000 5J0 

000 SOP 
"000 SFO 
000 

000 C 
000 C 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
ooo' 

000 

000 

000 

000 

000 

000 

000 c 
000 
000 
000 
000 
000 
000 

000 c 
000 
000 

000 100 
000 
000 
000 

000 200 

000 

000 

000 

000 

000 

000 701 

000 

000 

000 c 
000 
000 
000 
000 
000 
000 
000 

000 c 


7:01 :59-C0,) 

8 ATS6 U3200 

TION.8t9,20 

RTRAN 

SUBROUTINE f^AlU 

OF ALTERNATE SCAN OFFSETS* CORRELATIONS DONE OK IR • 
ASSIGN TAPE TO UNIT 10. OUTPUT IS TO PRINTER AND UNIT 
LOGICAL hISSNG tISLAST 
DIMENSION I0DATA(1299>,1£DATA(1299> 

DIMENSION MINdO) ,MOUT (24) 

DIMENSION I0&RTC2598) ,IEBRTC2598) 

DIMENSION WATc (250) rOFFS£T(230), ICO0RC230) tPLPARY(9> 
COMMON/OF FDAT/NCORMX, I CDORrOFFSETfWATE 
DATA ISTAT6/0/ ,KLAG/9/,NCOR/220/fLPRT76/ 

DATA M1N/6HGCCGCC»8*0/ 

DATA NW1/98/ ,NW2/ 1299/ ,NW3/3 732/ 

DATA LUN/10/ ,HISSN6/. FALSE./ 

DATA LUNF/20/ 

DATA lFIRST/oOO/,INTSCN/ 8 / , LS IS CN/400/ 

DATA IELEST/100/ »NTSIZ/15/ tNOMOFF/ + 8/ , IELINT/1 0/ 

DATA NBP/2593/,ID0FF/1 98/ 

DATA IDSTRT/67/ 

CALL la(MIN) 

IREEL=0 

INITIALIZE SCAN SEARCH 
rriRST=(CiriRST+i)/2>^2 
INTSCN=ClNTSCN/2)*2 
LSTS CN = C (LSTSCN+1 )/2)*r2 
IESCAN=IFIRST+INTSCN 
CALL GTAP(14,IRE£L,LUN) 

CALL REWCLUN) 

SKIP HEADER RECORDS 
CALL READW(LUN,KW1,IEDATA> 

CALL READW{LUN,Krf1 ,IEDATA) 

IFCMISSNG) I£SCAN=IESCAK-2 
IFC.NOT.KISSNG) IE S CAN=I ES CA N-INTSCN 
10SCAN=IESCAN-1 
I SLAST-IESCAN.LE .LSTSCN 
CALL READW(LUN ^Nv.1 tIEOATA) 

CALL CRKTHR(195,1EDATA,IEBRT) 

JESCAN=1E8RT(195) 

MISSN6=J ESCAN.lt. I ESC AN 

IF(r»ISSNG) ENCODE (132, 701 tMOUT)IES CAN 

IFCMISSNG) CALL TQ(MOUT) 

F0RMAT(1X,^SCAN NUMBER^rlS,' NOT FOUND. 

IFCMISSNG) GO TO 100 

IF(J£SCAN.GT*IESCAN) GO TO 200 

EVEN SCAN FOUND. NOW CHECK FOR ODD SCAN. 

CALL R£ADW(LUN,N*1 rIODATA) 

CALL CRKTHR(195,10DATA tIOBRT) 

JOSCAN=IOBRT(195) 

MISSNG=JOSCAN.NE.IOSCAN 

IFCMISSNG) ENCODE (132 ,701 ,MOUT) 1 OS CAN 

IFCMISSNG) CALL TQCMOUT) 

IFCMISSNG) GO TO 100 
BACK UP TO READ IN DATA. 


20 « 


■pOO®' 


ss« 



OFSTG/AOIPS 


JO0f=2 

DO 650 Iwtt=l,«lCOO 

nRITE.C6P,710HICt)ORU»i) ,OFPSEICIWJ ,wArE.flWl.Xw=IWri,NCOR,NlCO£.) 
650 CUNTIHUE 

710 FURMftr {C5(5X,14,lX,F8.3,F8.4))) 

NCOKMXsNCOR 
GO TO 999 
990 COKTlNOe 

WRirE(S,7l5J 

715 F0RMAri2X, 'OFFSblGEN bKftOR RblURH.*) 

999 CONTIWOB 

CADO lOTPRWCGUNJ 

CAhL REOUbSCRADSOC '0FSTF2' )3 

END 



OFSTG/AOIPS 


C ■ tVEN SCAN puu«0» «0W CHtCK FOK UOD aCAN. 

CAI<b lOTPlN tLUN, iOOAlA.HrVl.LE, JSI AI J 
CAUli CKtvlKkU95.XUDArA.IOBRI3 
JUSCAN=10aRT(19s3 
MlSSRG=JOS>CAM.«E.lUSCAt.l 
If CrtISS«G3riRir£tIT. 701JIOSCAN 
iFCMlSSiiG) GO fO 100 

c Back, up io ri£ao in daia. 

call IO-lPSR(bUN,“2,ISJ’Ai'} 

C INPUT- P.Alft UF SCANS, 

CAUL IOIPlNCGUN,Ifc,DATA#Nrt2.l.e,lSrAT3 
CAOL CKKrHR(l9&,lEDAfA,ieBRl) 

J£SCAN = it26RTC19S3 

CaMj CKKAlStNBP,iE0AXA,ieBRX3 

CALL I0TPXN(L0N.IU0AXA,hrt2,L£.ISXAX} 

CALL CKKTHRtl95,l‘ODAXA.10BRX) 

JOSCAH = iObRU1953 

CALL CRKAiS{NaP,lO0AlA,i0BRX3 

If'£ JeSCAh,Nb.I£SCAh.OR.JDSCAN.Nli.lOSCAII3 GO IQ 990 
C 

C - DO correlaiions. 
c • 

DO bOO KC0R=1,NC0R 

l£Le=l£LESr+f KC0R-l3?i£LINX 

ICOOR(K.CORj=ieLE 

«AXLAG=(NLaG-13/2 

rtlNLAG=-HAXLAG 

DU 500 KLAG=.'1INLAG,MAXLAG 

MAXbRI-0 

MlNtjRl = 512 

NTSIZ=CNrSI2-U/2£'2 + l 
FDp=0 

DO 400 KVAL=l.HrSIZ 

lCT=l£L6-fKVAL-l-NXSiZ/2 

lCS=^eLE+XVAL-l■^KLAG+NOMOFf-NXSlZ/2 

ICrx=lCX+lDOFF 

1CSX=ICS+100FF 

lPlXX=103RrtlCtX3 

ZPlXS=l£t5Rr(ICSX) 

MAX8Rl-=»AX0('lAXBRX,lPlxr3 
MINBRI-HINOC AlNfaRX, lPlXX3 
f Lp=FLP+CiPixI-XPiXS3*^2 
400 CONflNOE 

C STORE LP(2) measure: VALUE 

f LPARiCXLAG-MlHLAGtl3=FLP 
500 COnI'INUE 

c Search iable of values of lpc23 measure for min, 
XOfF=KUmOFF+»InLaG 
FLP.|41N=FLPARU1) 

DO 510 Ksi.NLAG 

IFtFLPMXK.GT.FLPARX (K.) 3 IOFF=NOHOf F+ 1 K-i 3 +HINLAG ' 
lF(FLPMlN.Gr.FLPARX(K)3 FLPMiW=FLP ART tK3 
510 CONCXSUE 

c height offsets bx brightness Range. 

»I = FG0AU«AXBR1-MXNBRT3/S12. 
wATecKCOftJ=HAlB(KCORJtWT 
OFFbEr£KCUK3=OE'lfSErCKCOR) + XOFF*HX 
600 CONTINUE 

IF( .NOr.lSLAST3 GO TO iOO 
C COMPUTE AVaRAGE OFfStTS. 

00 610 KC0K=1,NC0R 

if(rtAl'ECKCOKJ.LT.l.OE-J) GO TO 610 
OFFSET CK.CORJ=OFFSEICKCOK3/WAT£(KCOR3 
610 CUNTXNOb 

c oomp table values. 

nicol=iNcOr-u/5+i 



n n o 


OFSTG/AOIPS 


program 10 generate tasge 

OF AOTBRnATE scan offsets. CORRfcLAnONS done on ir. 

A5&1GM TaPt. 10 Ortil TEN, OUTPUT IS TO PRXNieK AND GliOSAO COMMON 
LOGICAL! Ml&S.NG.ISLAST 
Oi»t,NSluN IDOaTACSoOU) .leuATALSbOO) 

EOuIVALEnCBUODAI'A, leOATA) 

OImENSIOn (4iNll0J,'40UI(37l 
dimension lueKT(2bSR) .IESRT(2S383 

DiHbNSlOW NAlEt^iOl ,0» FSET 1 23 0 ) , ICOOR 1 230 ) , FLP ART C9) 
COMHUN/UUFFER/ftAlE.Of t SET 
COM'lON/BUFf 1/NC0 KmX,1C00R 

DATA 15rftrG/0/,NLAG/9/.CJCOR22 20/,LPRX/62 
DATA NNl/l<l7/,NN2/19 4y/,Nri3/5S98/ 

DATA LuN/iO/ .HISSNG/.FAL sE./ 

Data lp/s/,11/5/ 

DATA IFiRijT/SOO/.IhlSCN/ 8/ , LST5CN/400/ 

DATA l£l,EST/10O/.KlSlB/15/*NOKOFf / + 8 / , lELlNT/lO/ 

DATA HtlP/^598/,10UFF/1982 
DATA lOSTKT/67/ 

C initialize scan SEARCH 

lEIftSt=C lIFiRST + i J/2H=2 

itirSCN=(lrtrsCN/23^2 
LSTSCNst ILSTSCN+1)/2)*2 
ltSCAN = lflRJ>TtINT5CN 
«CORWJC=0 
DO 10 1=1,230 
« ATE U 3=0. 

0FFS£rCI3=0.0 
lC0ORtl3=0 
10 CUNTINUE 

't*ilCOObC7A,702,riOUf3 

702 fORMATtlX, 'MOUNT TAPE. ENTER DRIVE NUMBER'! 

MtNU=l 

call iNCOMtHQOT.SS.ia.MlN.MEwO.l.l) 
lHUNlT=MlH(i3 
IFtMSNU.LT.O) GO TO 9^9 
If (M£NU.EQ.0 3I.4UniT = 0 
CALL ASNLUNCLUn, 'MM' ,IM0N1T3 
Call iotphw(lun3 
C skip HEADER RECORDS 

Call lDTPSR(LUN,t2,IST<lT3 
100 IfihlSSNGJ 1E5 CAn=1ESCAN-2 

IF ( .NOl .M1SSNG3 1ESCAN=XESCAN-I«TSCN 
lOSCANsIbSCA.N-l 
ISLASr=iE5CAN.Lb.LSTSCN 
200 CALL ICjrPiN(LUN,ILDATA,NWl,LE,lSTAT) 

CALL CHAiHKU95,iEUArA,I£BRX3 

UtSCAN=lE3KT(19S3 

rtlSSNG=J£sCAN.H.iESCA« 

IFCMlSSiNG) .'iRlIbClI,70l3lESCAN 
701 fORMAUlX, 'SCAN nuHbbK' ,15. ' NOT FOUND.') 

IFfMiSSUG) GO TO 100 

IF ( JE5CAN.GT.1E5CAN3 GO 10 200 



****** OFSTF/MCIOAS ****** 


000056 

000057 

000058 

000059 

000060 

000061 

000062 

000063 

000064 

000065 

000066 

000067 

000066 

000069 

000070 

'000071 

000072 

000073 

000074 

000075 

000076 

000077 

000078 

000079 

000080 

000081 

OODOS2 

000083 

000084 

000035 

000086 

000087 

000088 

000089 

000090 

000091 

000092 


000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

OOQ 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 


SCLASO=XO 

SCLAS1=X0 

lELEHH=rLE 

IELEMX=IRE 

IX=H*(R-1}/2 

00 29 1=0, « 

COEFCl+1 )=HORK (I+1+IX) 

29 ASC0EFU + 1> = C0£F(I+1) 

C STORE DATA 

CALL WCOM 

C PRINT COEFFICIENT, DEGREE SCALING 

ENCOO£(132,720,nOUT>H 
CALL TPCJOUT.KOUT) 

720 .FORKATCIOX, "DEGREE='', IZ> 

ENCOOEC132,721 ,HOUT)X0,XD 
CALL TP(JOUT,MOUT) 

721 F0RKATC2X, 'SCALING PARAMETERS. X0=',£l5.9,' XD= 

DO 30 1=0, H . 

ENCODE (132, 722, HOUT)I,COEF Cl +1> 

30 CALL TP(J OUT, MOOT) 

722 F0RHAT(10X,'d(',12,'>=',e15.9) 
return 

990 CONTINUE 

ENCOdeCI 52,720,H0UT)IER 

CALL TQ(MOUT)- 

RETURN 

730 F0RKATC1X, 'ERROR COOE=',I2) 

ENOS 

SFILEHA 

DELETE GCCGCC.GORP 
SINCLUDE ATSSSP 
SCATALOG 

NAHE=GCCGCC,5,R,W,0 

TYPE=FG 

LIB=ATSFLB,LL 

BEGIN 

3E0 J 


E15.9) 


END ELT. 

aHD6,P ****** OFSTG/HClOAS ****** 




OFSTF/WCIDAS *^***<t 


QlLTfl AF^OFSTF/hCIOAS 


£LT007 

RLIB62 12/22-1 7: 

01 :58-<0,> 

000001 

000 

SJOB 

OFSTFT 03200 

000002 

000 

SOPTION .8,9,13,20 

000005 

000 

SFORTRAN 

000004, 

000 


SUBROUTINE MAIN 

000005 

000 

c 

PROGRAM TO FIT POLYNOMIAL TO ALTERNATE SCAN OFFSET DATA STORED IN 

000006 

000 

c 

UNIT 1l) BY program OFFSETGEK. INPUT DATAr 

000007 

000 

c 

CARDI: FIRST-POINT-ELEHENT-COORDINATE LAST-POINT 

oooooa 

000 


LOGICAL OPTION 

000009 

OOQ 


DIMENSION lCORLK(2>,XSCAL£n0) ,OFSTLM ( 2 > , YSCALE ( 1 0) 

000010 

000 


DIMENSION DATl <690>,rfORK.(231 ) 

000011 

000 


DIMENSION COEF (21 ) ,UABEL (12) 

000012 

000 


DIMENSION KINdO) ,M0UT(24) 

^00013 

000 


COHK.ON/0 FFEAT/NCORHX, I COOR (2 30 ), OFFSET (230) ,UATE(230) 

000014 

000 


COMKON/NAVCOH/NAVN,INAV, lYR, ID AY , TOTLI N , D EG L I N , TOT I EL , D E GEL E , PI CLI 

000015 

000 


*N,PICELE,TKPSCLt10YR,10DAY,TK,R1X,R1Y,R1Z,R10X,R1DY,R10Z,PITCH,ROL 

000016 

OOG 


*L,YAW,PTIK(3),TWh(3),TMX (3 ) , NL COE F (2) , NRC OE F C 2) , SCLLO( 2) ,SCLL1 (2) , 

000017 , 

000 


*ELC0EF(11 ,2) ,SCLR0(2) ,SCLR1(2) .ERCOEF (1 1 ,2) , NAS C£F , SCL ASO, S CLAS1 , 

000018 

000 


*IELEHK,1ELEMX, ASC0EF(16) 

000019 

000 


DATA ICORLM/0,2400/,OFSTLK/5 .,15./ 

000020 

ooo 


DATA LP/6/ 

000021 

000 


data MlNOEG/15/ 

0D0022 

000 


DATA KXN/6HGCC6CC,8*0/ ,LUNF/10/ 

000025 

000 


CALL IQ(MIN) 

000024 

ooo 

c 

READ LIMITS OF USEFUL DATA (ICOOR LIMITS) FROM CARDS. 

000025 

000 


ILE=HIH(1) 

000026 

000 


IRE=HIN(2) 

000027 

ooo 


J0UT=1 

000026 

ooo 


1 F (OPTION (niN( 3) ,3H P>) J0UT=2 

000029 

000 

c 

READ IN OFFSETS AND WEIGHTS AND COORDINATESS (ELEMENT NUMBER) 

000050 

000 


CALL OYNASGC'OFFSTD'.LUNF) 

000051 

ooo 


CALL OPN(LUNF) 

000032 

ooo 


CALL R£ADW(LUNF,1151,NC0RMX) 

000033 

ooo 


DO 10 I=1,NC0RMX 

000034 

ooo 


IF(1COOR(I).£Q.ILE) ILEAE=I 

000035 

ooc 


1FCIC00R(I).EQ.IRE)IREAE=I 

000036 

000 

10 

CONTINUE 

000037 

ooo 


NPTS=1REAE-ILEAE+1 

00003S 

ooo 

c 

FIT LEAST SQUARES POLYNOMIAL. 

000039 

ooo 


IN=0 

000040 

ooo 


DO 25 1=ILEAE,IREAE 

000041 

ooo 


1N=IN+1 

000042 

ooo 


DATI (IH)=ICOOR (I) 

000043 

ooo 


DATl (IN+NPTS)=OFFSET(I) 

000044 

ooo 


DATI(IN+2*NPTS)=UATE(I) 

000045 

ooo 

25 

CONTINUE 

000046 

ooo 


CALL ApCK(OATI ,NPTS ,MI ND EG ,X D , XO , WORK, lER ) 

000047 

ooo 


IFUER.NE.D) GO TO 990 

000048 

000 


ETA=1 .UE-3 

000049 

ooo 


EPS=1 .OE-4 

000050 

ooo 


10P=-1 

000051 

ooo 


CALL APFS(HORK,N1NOEG,M,10P,EPS,ETA,IER) 

000052 

ooo 


1F(IER)990,28,28 

O0OO53 

ooo 

28 

CONTINUE 

000054 

ooo 


call rcoh 

000055 

ooo 


NASCEF=M 



OFSTF/AOIPS 


iELfMXslKli 

l>0 29 1=0»M 

CD6.rCI + U=rtORKCl + H-lXJ 

29 AScO£K(ItU=COc.FcI+n 

C PftlNl’ COEFFlClEN’i', OEGRee SCALING 

wRirt;(ir,720jM 

720 fuRMftrClOX, 'Ofc,GREfc=' .12J 

• «Riiecn»72nxo.xo 

721 fORMAl C2X, 'SCALING PAKAHEfERS. X0='.E15.9,' X0=',E15.9) 
1)0 30 1=0. M 

30 rtRlTliCir, 722)1, COEP(l + l) 

722 f ORHATCIOX, 'DC ' ,12, ')=' .E15.9) 

GO TO 999 

.990 COHriNOE 

hKixeUi’,730)IER 
999 continue 

CALL RbQUeS(KA050C'AISF2')) 

730 FORHAXdX. 'ERROR CODES' rl2) 

END 



o o o 


OFSTF/AOIPS 


C 


10 

C 


2 £> 


28 


PKOGRAM ru Fxr ROCjiUOrtrftL TO ALiTE-KNATE scan offset data stored in 

O.NIT lo BJ£ PROGRAM OFfoblGErt. INPUT DATA: 

CAROi; flR&X-pOiNT»eiiE.'l£.Ni-C0ORDiNATE LAST-POINT 
logical UPIXON 

DIMENSION ICORLMtal rXSCALEUOI ,OFSrLM(25 ,ISCALE(10> 

Ol'JENSlUN OAI1C690J ,w0BA(231) 

DIMENSION COEc (2U .LAbELC12) 
dimension HlriaO),MOUl(37J 
OIHEnSIo.I 1COUK1230) , OF FSEXC 230 ) . W ATE C230 J 
COMMOS/OUFfER/wAJE, OFFSET 
CUMHON/SOFFI/nCOkMX, ICOOR 

COMMDN/WAVCJH/NAVN, INA\^,12R,IDAI .TOTLIN.DSGLIN.TOTIEL.DEGELE.PICLI 
♦ N,PIC6Le, THPsCL, 10YK.IU0AI, TM,RlX,RlI,Ki'/;,HlDX.KlLY .hlOZ*PIfCri.ROL 
♦L» lAN»PriMC3J .TMNC3) . IriA C 3 ) , N LCOEl- ( 2 J f NKCOEF C 2 1 . SCLLO C . SCbLl { 2) . 
‘fELCOEFai,2J .SCLR0C2) .SCLKXC2) .tRCOEF(I1.2) .HASCEFfSCLASO.SCLASI. 
♦lELeMN.IELEMX. ASCOEFC16) 

DA I A ICORLMFO, 2AOO/,OFSILH/5. , 15./ 

DATA LP/6/ 

Data IT/S/ 

DATA MlNDEG/15/ 

HEN0=2 

CALL IrtCOMC’ ENTER LEFT. RIGHT ELEMENTS FROM PRIWIEK LISTING. 49 . 
>t3,l.MIN,l,l) 

IFfHENO.LE.O) GO TO 999 

READ LlKllS OF USEFUL OATa (iCOoR OimiTSI FROM CARDS. 

lLE=MlNtl J 

IK£=MI«C2I 

keAO in OFFSETS AND iiEIGHXS AND COOR01NATESS(ELEM£NT NUMBER) 

DO lU I = I,NCUF(MX 
IFCICOORCD.EQ.ILE) 1LEAE = 1 
If (ICOORCl J .Ea.lRE) IKE aE=I 
CONI XNOE 

NplS=lRBAE-lLEA£+i 

FIT LEAS! SQUARES POLYNOMIAL. 

lix^O 

DO 25 I=iLeAE, IREAE 
INslrttl 

oAritiNjsicooRtn 
U Am i N+ h PXS ) sOf FSeT { 1 ) 

DATl(lN+2#NPlSJ=rtATG(I) 

CONllNUE 

CALL A?Crt(OAlI.«Pi5FMlNDEG,XD,XO,WORK.IER) 

XFCI’ER.NE.O) GO 10 99u 
EtA=l.UE-3 


EpS=l.0E-4 

I0P=-1- 

Call APFS(f(OKK,nlNOEG,M,iaP,EPS,ErA,lBB) 

If {ieRJ990,28,28 

CONTINUE 

NASCc.F=H 

SCLrtSQ=xO 

SCLAfal=XD 

1ELEMN=1L£ 


S- 



ik-tcicirit * 

NRHLI2/HCI&AS 

* *** iir * 

QEtTiL 

AF«.NRMLIZ/HCID AS 


ELT007 

RLIB62 12/22-17: 

01 :57-(0O 

000001 

000 

SUBROUTINE NRMLIZ (VX,VY ,VZrVNORH) 

000002 

000 

V NORM -SORT (VX**2 + VY**2+VZ**2> 

D00003 

000 

VX=^VX/VNORM 

000004 

000 

VY^VY/VNORH 

000005 

000 

V2^V2/VNORM 

000006 

000 

RETURN 

000007 

000 

END^ 


END ELT. 

3HDG,P ****** OFFSETFIT ****** 



****** ftINKIZ/MCIDAS ****** 


SELTrL A F •MINH 1 Z /H CIt> AS 

ELT007 RLIS62 1 2 / 2 2-1 7:01 ; 56- (0 t ) 


000001 

000 


SUBROUTINE HINHIZ (PTIN ,PTOUT ,G NORR ,VAL , ITN> 

000002 

000 

c 


000003 

' 000 

c 

THIS ROUTINE PERFORMS A MODIFIED NEUTON KETHOD MINIMIZATION. IT 

000004 

000 

c 

BEGINS AT THE POINT "PT1N“ AND RETURNS IN "PTOUT" THE POINT 

000005 

000 

c 

SELECTED AS THE OPTIMAL POINT, IN "GNORH" THE NQRM OF THE GRADIENT 

000006 

000 

c 

OF S AT- “PTOUT”, ANO IN "VAL“ THE VALUE OF THE 09JECCTIVE FUNCTION 

000007 

000 

c 

S AT “PTOUT". A POINT XCK+1) IS DEEMED OPTIMAL WHEN 

000008 

000 

c 

ABS(S(XtK+1)>-SCX<K)))<=<10**-10)*ABSCStX(K))>. 

000009 

000 

c 


0,00010 

000 


LOGICAL ALTRET 

000011 

000 


DIMENSION PTIN (3) ,PT0UT(3) 

000012 

000 


DIMENSION HESS (3,3) ,GRAD(3) , PT (3> ,DIRC3) 

000013 

000 


DATA CONVRG ,ITERAT,EQUAL0/1 . OE-18,25,1 .OE-30/ 

000014 

000 

c 


000015 

000 


ITN = 0 

000016 

000' 


DO 5 1=1,3 

000017 

ODD 


5 PT(I)=PTIH (I) 

000018 

000 


OLDVAL=S (PT) 

000019 

000 


DO 50 1=1 rITERAT 

000020 

000 


1TN=ITN+1 

000021 

000 


CALL PRTIAL(PT, GRAD, HESS) 

000022 

000 


CALL INVERKHESS, HESS, 3, ALTRET) 

000023 

000 


IF(ALTRET) GO TO 16 

000024, 

000 


DO 10 J=1,3 

000025 

000 


DIR(J >=0 

000026 

000 


DO 10 K=1 ,3 

000027 

000 


10 DIR(J)=DIR(J)-HESS<J,K)*6RAD(K) 

000028 

000 


DDD=0 

000029 

000 


DO 15 J = 1 ,3 

000030 

000 


15 DDD = DDD + D1R (J )*GRADU) 

000031 

000 


IF (DDD.LT.-EOUALO) GOTO 25 

000032 

000 


IF(DDD.GT.+EQUALO)GOTO 19 

000033 

' 000 


16 ODD=0 

000034 

000 


DO 1 7 J=1 ,3 

000035 

000 


DDD = DDD + GRAD CJ )**2 

OOOD36 

000 


17 D1R(J)=-6RA0CJ> 

000037 

000 


GOTO 25 

000058 

000 


19 DO 20 J = 1 ,3 

000039 

000 


20 OIRtJ )=-DIfi (J) 

000040 

000 


25 CALL LSCPT, PTOUT, VALjDDD, DIR) 

■000041 

000 


1F(ABS(VAL-0LDVAL) .LE.C0NVR6*A9S(0LDVAL)) GO TO 60 

000042 

000 


OLDVAL=VAL 

000043 

000 


DO 30 J=1 ,3 

000044 

000 


30 PT(J)=PTOUTCJ) 

000045 

000 


50 CONTINUE 

000046 

000 


60 GNOR«=0 

000047 

000 


DO 65 1=1,3 

000048 . 

000 


65 6N0RH=GN0RH+6RA0(I)**2 

000049 

000 


GNORM=SQRT(GNORH) 

000050 

000 


RETURN 

000051 
END ELT* 

000 


ENDS 

aHDGfP 

****** 

NRMLl Z/MCIOAS ****** 



****** LS/HCIDAS ****** 


SELT,L 
ELT007 
000001 
000002 
OOOOG3 
000004 
000005 
000006 
000007 
000008 
000009 
000010 
000011 
- 000012 
000013 
000014 
000015 
000016 
000017 
000018 
000019 
000020 
000021 
000022 
000023 
000024 
000025 

END ELT 
8 HD 6 »P 


AF*LS/MCIDAS 
RLIB62 12/22-1 
000 

000 C 
000 C 
000 C 
000 C 
000 C 
000 C 
000 c 
000 c 
000 c 
000 c 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 


7:01 :56-CO, ) 

subroutine LSCXfY ,VAL,DD»D1R> 

THIS ROUTINE PERFORMS AN ARHIJO LINE SEARCH FROM THE POINT “X” 

IN THE DIRECTION ‘*DiR** AND RETURNS THE SELECTED POINT IN 

AND THE OBJECTIVE FUNCTION VALUE S<Y) IN ”VAL**. ON CALLr ”DD” 

IS THE UNNORKALIZED DIRECTIONAL DERIVATIVE <GR A D CS ( X ) 1 , DI R> • 

THIS LINE SEARCH ROUTINE RETURNS IN Y THE POINT X+ 2** (-N) * DIR 9 
WHERE N IS THE LEAST NONNEGATIVE INTEGER SUCH THAT 
-S(2**(-N)*D1R) represents AT LEAST 40% OF THE FUNCTIONAL DROP 
IN THE LINEARIZATION OF S AT X IN MOVING FROM X TO X+ 2** ( -N > *D IR o 

DIMENSION X (3> #YC3) r0lR<3) 

DATA FAC, WALK/5. uE-1,1 oQE-5/ 

RLAH = 1 .0 

TSTVAL = .4*RLAH*ABS (DD> 

OLDVAL-S (X) 

1 IF(RLAM.LT»WALK)RETURN 
DO 2 1-1 t3 

2 Y CI)-XCIMRLAH*DIR (I) 

VAL=S(Y) 

1FC0LDVAL-VAL.GE.TSTVAL)RETURN 

tstval-fac*tstval 

RLAM=FAC*RLAH 

GOTO 1 

ENDS 


****** 


MINMIZ/MCIDAS 


****** 


ORIGINAt. PAGElB 
OF POOR QUALira! 



LATSF/MClDAS ****** 


000170 

000 

ENDS 

000171 

000 

$FI UEHA 

000172 

OO'O 

DELETE LOATSf,GORP 

000173 

000 

SINCLUOE IFLD 

00017A 

000 

SCATAL06 

000175 

000 

NA«E=LDATSF,5,R,U,D 

000176 

000 

TYPE=FS 

000177 

000 

Llb=ATSFLB,LL 

000178 

000 

BEGIN 

000179 

000 

SED J 


END ELT* 

SHOGfP ****** LDATSf ****** 





****** 

LATSF/HCIOAS 

000113 

000 


000114 

000 


000115 

000 

921 

000116 

000 

701 

000117 

OOQ 


000118 

000 


000119 

000 


000120 

000 ' 


000121 

000 

922 

000122 

000 


000123 

000 


0001 24 

000 


000125 

000 


.000126 

000 


000127 ‘ 

000 


0001 2S 

000 


000129 

000 


000150 

000 


000131 

000 


000132 

000 


000133 

000 


000134 

000 


000135 

000 


000136 

000 

c 

000137 

000 


000138 

000 


000139 

000 


000140 

000 

50 

000141 

000 


000142 

000 


000143 

000 

100 

000144 

000 

c 

000145 

000 


000146 

000 


000147 

DOD 

150 

000148 

000 


000149 

000 


000150 

000 


000151 

000 


000152 

000 


000153 

000 


0001 54 

000 


000155 

000 


000156 

000 


000157 

000 


000158 

000 


000159 

000 


000160 

• 000 


000161 

000 


000162 

000 


000163 

000 


000164 

000 


000165 

000 


000166 

OOD 


000167 

000 


000168 

000 

100 

000169 

000 



****** 


CALL EHESS(3HREQ,0) 

RETURN 

WRITE (LP|701 ) 

FORMATCIX, 'E-O-F TERMINATES LOAD") 

CALL MARKOKCIAREA) 

CALL MARKOK (IRAREA) 

CALL REWCLUN) 

RETURN 

CONTINUE 

CALL EHESS(3HRE«J,0) 

RETURN 

ENO 

SUBROUTINE LINGR8 { INO AT A ,SOI R ,ELE , NELES ,BO F t OUTDAT) 

IMPLICIT INTE6ER(A-1) 

DIMENSION INDATAd ) tOUTDATd ) 

DIMENSION TOAT(2406),DLELE C672) 

COMMON/OFFSET/DLLLE 
DATA WOFF/6/ ,NE/2A06/ 

IFCBDF.NE.1) CALL EKES S ( 3HREQ , BD F ) 

IF(BDF.NE.I) CALL EXIT 
I02=ELE-1 

CALL CRKATS (NE , INDATA, TOAT) 

I FISDIR.EQ .1 > GO TO 100 

SHIFT EVEN SCANS 

DO 50 1=1, NELES 

DE=1 + I02+DLELE (D+WOFF 

DUTDAT{I)=TDAT(DE) 

CONTINUE 

CALL PACK (NELES, OUTDAT , OUTDAT) 

RETURN 

CONTINUE 

load odd scans 

DO 150 1=1, NELES 
0UTDAT(I)=TDAT(I+I02+U0FF) 

continue 

CALL PACKCNELES, OUTDAT, OUTDAT) 

RETURN 

END 

SUBROUTINE 6EN0FF (lUELE, NELES) 

COMHON/NAVCOM/NAVN , IN A V , lY R , 1 0 AV , TOTLIN , DEG LIN , TOTIEL , DE 6ELE , PICLI 
*N,PlCELE,TMPSCL,IOYR,IODAY,TH,RlX,Rly,RlZ,RlOX,RlOY,RlDZ,PITCH,ROL 

*L,YAW,PTIM(3) ,THN(3),TMX (2 ) , NLCO E F C2 ) ,NR COE F (2) , S CLLC ( 2) ,S C LLU2 ) , 
*£LC0EFd1 ,2) ,SCLRO (2) ,SCLR1 (2) ,ER COE F ( 1 1 , 2) , N AS CEF , S CL AS 0, S CLASI , 
*I£LEMN,1EL£KX,ASC0EF(1 6) 

COMHON/OFFSET/IDLELE(672) 


DATA NOHOFF/8/ 

CALL RCOH 

NEN0=IU£LE+NELES-1 
00 100 1=1, NELES 
IDLELE(1)=N0M0FF 
IELE=IUELE+I-1 

I F (IELE*LT . lELEMh .OR.I ELEHX .lt .1 ELE)GO TO 100 

T=I£LE*SCLAS1+SCLAS0 

CALL CNPS (Y,T, ASCOEF,NASCEF> 

IDLELE (I>=Y + 0.5 

CONTINUE 

RETURN 


ORIGINAL PAGE IS 
OF POOR QUALITY 



****:** LATSF/HCIDAS 


•k it it -k 


000056 

000 


:FlLE=IFIt.E-l 

' 000057 

000 


IFUFILE.LT.DGO TO S02 

000058 

000 


GO TO 922 

000059 

000 

502 

CONTINUE 

000060 

000 

c 

ADVANCE TO FIRST SCAN TO INPUT 

000061 

000 


CALL REAOUCLUN,NWA,IARRAY> 

000062 

000 


CALL READW<LUN,NWA,1ARRAY) 

000063 

000 


DO 510 1=1 ,1200 

000064 

000 


CALL REAOW(LUN,N-A,rARRAY) 

000065 

000 


CALL CRKTHR(19S,IARRAY,1AR0U> 

000066 

000 


INHS=IAR0W(19S) 

000067 

000 


IFCINMS.LE.IUMS) GO TO 511 

000065- 

000 

510 

CONTINUE 

,000069 

000‘ 


■ 60 TO 922 

* 000070 

000 

511 

CALL BSR(LUN) 

000071 

000 

c 

START MAIN LOOP TO READ DATA IN AND STORE 

000072 

000 

c 

READ A SCAN 

000073 

000 

601 

CALL READU (LUN,NWA rIARRAY} 

000074 

000 


CALL CRKTHR(195,1ARRAY,IAR0W) 

000075 

000 


1NMS=IAR0W(195) 

000076 

000 


If (INHS.LT.TLHS) 60 TO 690 

000077 

000 

c 

HAP SCAN NUMBER TO AREA ROU 

000078 

000 


L1=C1 200-INHS>*2+C3-2) 

000079 

000 


L2=(120D-1NHS)*2+(3~1) 

000080 

000 


LA1=L1-IULINE 

000081 

000 


LA2=L2-IUL1NE 

000082 

000 


LR1=H0D(LA1 tIBDF) 

000083 

000 


LR2=MOO (LA2 ,IBDFJ 

000084 

000 


IFCLRI.KE.O .AND. LR2.NE.0) GO TO 601 

000085 

000 


IROWI =LA1 /IBDF 

000086 

000 


IR0H2=LA2/IBDF 

000087 

000 


NS=NSECL(NELES) 

000088 

000 


ISD=H0D(INHS,2J 

000089 

000 


1SEC1=NS*IR0«1 

000090 

000 


ISEC2=NS*IR0W2 

000091 

000 

c 

PICK OUT LINE SEGMENT 

000092 

000 , 

c 

FIRST VISIBLE LINE 

000093 

000 


IF(IFVIS.AND.LRI.EQ.O) CALL LI NGRB (I ARRAY ( 2697) , IS D ,I.UELE,N ELES , IB 

000094 

000 


*DF,IAR0H) 

000095 

000 


IFUFVIS.AND.LRI.EQ.Q) CALL WR IT A Cl ARE A , I SEC 1 ,NS*1 1 2 , 1 AROU I 

000096 

000 

c 

second visible line 

000097 

000 


1FC1FVIS.AND.LR2.EQ.0)CALL LIN6RBCIARRAY(1297>,ISD,IUELE ,NELES ,IB0 

000098 

000 


*f ,IARO«) 

000099 

000 


- If (IFVIS.AND.LR2.EQ.0) CALL UR IT A Cl A RE A , I SEC 2 ,NS *1 1 2 , 1 AR OH) 

000100 

000 

c 

INFRARED 

000101 

000 


IFCIFIR) CALL LING RB C I Afi RA V C 97 ) , ISO , lOE LE ,NELES ,1B DF , I AR OH) 

000102 

000 


IFCIFIR.AND.LRI.EQ.O) CALL W RI TA C I R AR E A , I SEC 1 ,KS *1 1 2 , 1 AROH) 

000103 

000 


IFCIFIR. AND. LR2.EQ.0) CALL H RI TA C IR AREA , 1 SEC 2 ,NS*11 2 , 1 AROH) 

000104 

000 


GO TO 601 

000105 

000 

690 

CONTINUE 

000106 

000 


CALL MARKOK ClAREA) 

000107 

000 


CALL MARKOKCiRAREA) 

000108 

000 


CALL REUCLUN) 

000109 

000 


RETURN 

000110 

000 

910 

CONTINUE 

000111 

000 

912 

CONTINUE 

G00112 

000 

913 

CONTINUE 



****** LATSF/HCIOAS ****** 


SELTjL AF.LATSF/HCIDAS 

ELT007 RLIB62 12/22-1 7:01 £S3-(0, ) 


000001 

000 

SJOB 

LDATSF U3200 

000002 

000 

SOPTION .8r9t13»20 

000003 

. 000 

SFORTRAN 

000004 

000 


SUBROUTINE MAIN 

000005 

000 

C 

PROGRAM TO READ IMAGES FROM ATS-6 TAPES# 

000006 

000 


LOGICAL IFVIS , IFIR .OPTION 

000007 

000 


DIMENSION HIN( 12) , lARR AYC3732) ,IAR0U(672) 

000008 

000 


DATA LUN/10/,LP/6/, ICR/5/ 

000009 

000 


DATA MIN/6HLDATS Ft 10*0/,NWA/3732/ 

000010 

000 

c 

PICK OP INPUT PARAMETERS 

000011 

000 


CALL la(HlN) 

-000012 

000 

c 

INPUT: SSyyODO hhpmss area line element 

000013 

000 


10Ay=HIN (1 ) 

000014 

000 


1TIKE=MIN(2) 

000015 

000 


IAR£A=«IN (3> 

000016 

000 


LINE=H1N C4) 

000017 

000 


IELE=«IN (5) 

000018 

000 


1B0F=1 

000019 

000 


IREEL=0 

000020 

000 


KEY=3H 

000021 

000 

c 

COMPUTE DERIVED PARAMETERS 

000022 

000 


ISS=IDAY/100000 

000023 

000 


IF1LE=1 

000024 

000 


lVSS=(lSS/2>*2 

000025 

000 


1RSS=1VSS+1 

000026 

000 


1FV1S=ISS*E0.IVSS 

000027 

000 


1FIR-1SS.£Q*IRSS.0R.IAREA*LT#9 

G00028 

000 


IRAREA=1AREA 

000029 

000 


IF(IFIR.AND.IFVIS) IRAREA=1ARE A+8 

000030 

000 

c 

PICK UP AREA SIZE 

000031 

000 


CALL HOWBIGCIAREA.NLINS.NELES) 

000032 

00 0 


1F(NELES.GT,672) 60 TO 922 

000033 

000 

c 

COMPUTE UPPER LEFT COORDINATES 

000034 

000 


1ULINE=LINE 

000035 

000 


IUELE=IELE 

000036 

000 


1FC0PTI0N(KEY,3H O) lULlNE =L INE-CNLIN S / 2 ) *IBD F 

000037 

000 


IF(0PT10NCKEY,3H O) IUELE=IELE-(NELES/2)*IBDF 

000038 

000 


IUMS=1200-(1ULINL-1 )/2 

000039 

ooo 


1LLINE=1UL1NE+{NLINS-1 )*IBDF 

OOOOAO 

DOO 


ILMS=1i:00-(rLL INE-1I/2 

000041 

000 

c 

CREATE IR AREA FOR COMBINED LOAD 

000042 

DOO 


IF Cl FIR. ANO.IFVIS) CALL AR AS 12 Cl A R EA + 8 , NLINS , NELES ) 

000043 

ooo 

c 

ESTABLISH AREA DIRECTORY 

000044 

DOO 


I6E = 0 

000045 

ooo 


CALL ENAREACIAREA.IREEL, IDAY .ITIHE . lUL INE , lUELE , IBD F , IBD F . I GE) 

000046 

DOO 


I RREEL=IRSS*100000 + NOD CIREEL ,100000) 

000047 

ooo 


IFCIREEL.EQ.O) IRREEL=0 

000048 

DOO 


1RDAY=IR SS *100000 + HOD ( IDAY ,1 00000) 

000049 

ooo 


IFC1FV1S.AND.IFIR) CALL £N AR EA Cl AREA+8 , IRRE E L, 1 RO A Y , ITIH E, lUH NE , I 

000050 

ooo 


*UELL,IBDF,I8DF,IbE) 

000051 

ooo 

c 

GENERATE OFFSET TABLE 

000052 

ooo 


CALL GENOFFCIUELE, NELES) 

000053 

DOO 


IRR=HAX0C1 .MOD CIREEL, 10000)) 

000054 

ooo 


CALL GTAPC14,IRR,LUN) 

000055 

ooo 

c 

advance, TO file 


original 
OR POOR 



LATSP/AOIPS 


WKlTtCIT, 711i 

7U FURriAXC’ TAPE SEARCH STARTS «J 
CAfit, 101PiiKCIiil«,4-^, ISTAtJ 
uCi 510 lsl.1200 

call lOlPirt[Dil«,lARRAT,Hn'A.I.e.ISIAr) 
call CH*''TrtRU95.iAR«AT.lAROriJ 
I«rtS=lARUH£l95) 

IFUNMS.LE.XuMSJ'GO to 511 
510 CONTINUE 

GO TO 922 

511 call lOlPsRCLUN.-l.ISlAXl 

rtRlTEar,712J 

712 fORMAtC irtAGE LOAD SIARXS*) 

C start main loop 10 READ DATA IN AND STORE 

c read a scan 

601 CALL 101P1N(LUn,1ARRAT,NNA,LE.1S7AT1 

call CRRTHRCl95,iARRAT.IARON7 
lN*lS=IAR0ft(195) 

It (INHS.LT.IL mS) go To 690 
c MAP SCAN Number to area rOw 

L1=(120o-1NM5J^2^(2-2J 

L2=U200-lNMSJ^2+C3-n 

LAlsLl-’lOLINE 

LA2=L2'’I0LinE 

LB1=HODCLA1,1BOF7 

LR2 = (400CLA2.1BDn 

IflLRl.NE.O .AND. LR2.NE.07 GO TO 601 

lROftl=LAl/IBDF 

lROrt2=LA2/IbOF 

ISO = rfODUN«3,2) 

iSECl=hS<=IRQKl 

1SEC2=nS<.xrOk 2 
c PICK our LINE segment 

C URST visible LINE 

IF (IFVIS.AND.LRI.EQ.O) CALL LINGR8 t IARRAT C 37 457 , ISO , lUELE, NELES, iB 
^DF .IAR0«7 

IFflFvlS.ANO.LRl.El.OJ CALL /.RlXtC LRN vlS. I AROh , 0 7 
C SECOND VISIBLE LINE 

lFaFVl5.AND.LR2.Ea.0 7CALL L1NGR8 C I ARRAl a 945} . ISD. lUbLE. NL lES . IBD 
»F, JAR0N7 

IFdFVXS. AND. LR2.su. 0} CALL aRI 1 EC LRN VIS . 1 AROw. 07 
C INFRARED 

IFClflRI CALL LINGR3ClAKRAni457 .ISO.lUELE.NELES.iaOF.lAROW) 

IFCIf IR.ANO.LRl.EU.O) CALL NRlTb CLKNl R , I ARQW ,0} 

if(ifir.aad.lr2.eu.o) call wriieclRnir.iaroh.oj 

GO 10 601 
690 CONTINUE 

922 CUNllNUE 

CALl CLOSef CLrNvIS} 

CALL CLOStFCLRtURl 
CALL lOlPRWCLUNJ 
GO TO 999 
999 continue 

CALL KEU0ESCRAD50(*ArSF2')7 
END 



LATSF/AOIPS 


C 

10 


12 

13 


C 


C 


llihS3l-i;00-tiiJL,iNk,-n/2 ■ 
EbTABl^iSH area DiRECrORX 
DO iO 1=1.256 
lDlRtD=0 
lDIRaj = UI‘ 
lDIR(2)s»Sb ' 

I0IR(3J=' V 
1DIR(41=»1S' 

I0lK(3J='e.* 

10 Ik(oJ='H, • 

1D1R(71='1. » 
iDIRC»J=' • 

ID1R{9J=' At' 

10XR(IOJ='S6‘ 

10IRC11)=' ' 

XDlRtl2)='E.> 
lDll<a3)s'rf. ' 

1DIRC14J='T.' 

XD1R(15}=' » 

iOlRUbls' « 

DO 12 1=12. 1,-1 
lf(lDAl[.GX.MOOAl!(Ij)GO XO 13 
COHXIilUC 


rtUN=I- 

IOJR(l9)=llR4'100trtON 

101K(20 J=C10ai-H0uAX C«ON l)*100+lFlxa>Xl«tIPIC)/60.) 

I(i=PTl«(lPIC) 

lS = Prirt(li’ICj't-oO.-lMt-60, 

X&=MOUClS,bO) 

1M=V,00(XM.60J 
1D1RC21 J=lH*100<-IS 
10lft(24J=l 
1D1R(25)=1 
1D1R(26J=1PIC 
lDlftt2s/J = 'Ef ' 

101RC30)=' IC» 

1DIRCj1)=1 

l!>IRC32)=IjRNns 

loiRCiijssia 

IU1R(34)=513 

IDIRC37 1=512 

1UXR(38)=512 

I0lB(39J=lUt.be 

I01B(aOJ=1UHN£ 

10IB(41)=1 
Il>m(42) = l 
10IRC43)=1 
1DiB(44)=1 
IDIRC4bl=l 
lDIRt4bJ=l 
lDIRf47 1=1 
lt>lRC4bJ=l 

IDl R 151 1=10004(1(4 <11 NTti) 

1U1R(531=1 

NW0=(«EIitb4ll/2 

CALL, OPthCbRi'JVlS.NLlNS + l.NivOl 

CAbb OPtfUbRMlK.NblwS+l ,NWO) 

CAbb bcsb«tKl(ut<NVlb.XulR,2561 

XDlR(32J=bRfUR 

iDIR(3)=' ' 

XDIR(41='IR> , 

CAbb bij(.rtRrtbR(viR,10iR.256) 

oeneRAie offobr XAUbt. 

CAbb GbBUPf (lUtbE.HtbbS) 

CAbb ASnbONCboN, ''in' .IKUNin 
AOVAt^Ci!. 10 FlKbC bCfcpJ XU IRPUX 


ORIGINAL PAGE i: 
OF POOR QUALITY 



lATSF/AOIPS 


C PKOGRAri rO READ IMAGES FROM ATS-5 TAPES. 

iNTEGea^A XTIMX 
LOGICAL XFViS.lflR.OP'IION 
OlMENSloM .RMihUO) ,1ARBAX{5600J 

EOUIVALEiMCE {IDiRflARuw} 

DIMEWSIOU 10IRt256J ,1AR0WC'C56) 

OXMtNSlOW MOOAYC12) 

CUMHUW/«AVCOH/nAVM, INAV.lIR.lOAi ,TOTLIN.OEGHN,TOI1EL,DEGELE,PICLI 
♦ n,PlCELE, I’MPSCL.lUiR.lOOAX, TM.RIX, Rii.Ril, RIDA, HlDT. RIDE, PITCH. ROL 
Uw.PriMO) , iKivOJ ,irtAC3J ,NbCOEFt2) ,hRCOEF(2) .SCLL012J .SCLLl t2) . 
fELCOEf Cll.iJ .&CLR0C2) .SCLR1C2J .ERCOtt Cl 1 . 2 1 . H At>CLt , 5 >CLASq , bCLASX. 
♦lELEhN.XELErtX.ASCOEFCio) 

DAJ A KUDAlf/0, 31 ,39,90, 120,151.181,212.243,273.304,334/ 

Data LUN/t/.LP/e/.xcR/s/ 

DATA ll/5/,«ELES/5l2/ 

DATA EjLXhS/512/ 

C PICK OP IMPur PARAMETERS 

NrfA=5600 
Kt«U=l 

CALL IfiCOHC EtiTER PICTORE TIME (HrfMMSS) ' , ?8 • 1 . 1 , RHXN , MEMO . 3 . 11 

It (MENU.LE.O) GO TO 99» 

lllriX=KMllUl) 

Pi‘=FTi.MEUrXMXl 
DO 5 1=1.3 

5 IFCPr.EQ.PrXMCll) GO TO 6 
WRirECIT.7101irXMX 

710 FORMAItlX, 'NAVCOM NOT INIIXALXZEO FOR ',X7l 
GO 10 999 

6 CONlXiloe 
IP1C=I 
MbN0=2 

CALL INCOKC* EWIER START LINE AND ELEMENT 29 , 3 , 1 , MIN , MEHO. 1 . 1 ) 

If lMc.Na.LE.0) GO TO 999 
L1NE=M1NCI) 

Xfc.L£=rtX«(2) 

MEN(J = 1 

call XnCOM{I enter EOOtt AREAC1-7J ' .21, 5,1, HiN, menu, 1,11 
IF(hENU.LE.O) GO TO 999 
XAKEASMINC 11 

bRSVlSsIARtA + S+UPIC-D^T 

LRNIR=LRnyiS+21 

MENU=1 

call INCOMC Enter tape unit co.ll ' ,22,7.1.MIN,MENU,1,11 

IP (MENU. LE. 01 GO TO 999 

X«ONXl=MlN(ll 

moF=i 

. XfVXS=.TRUE. 

IFIRs.TROE. 

C COMPUiE UPPER LEFT COORDINATES 

1ULI«E=L1N£ 
lUELE=lELe 

lUajs = 12u0-(X0LIH£-l)/2 
ILHNE=iULXNE+(HLiN5“l)tIBDF 



****** INVERT/MCIDA S ****** 


selt,l af.invert/mcioas 

ELT007 RLIB62 1 2 /2 2-1 7; 01 : 5 2- ( 0 ♦ ) 


000001 

000 


SUBROUTINE INVERT (AAtB »N ,ALTRET) 

D00002 

000 

c * 


000003 

000 

c 

THIS ROUTINE RETURNS IN B THE INVERSE OF THE N-!> IHENSI ON AL MATRIX AA 

000004 

000 

c 

IF AA IS SINGULAR, A RETURN IS TAKEN THROUGH LABEL S. 

000005 

000 

c 

' 

000006 

000 


LOGICAL ALTRET,LOGTHP 

000007 

000 


DIMENSION AA(N,N),B{N,N),A(3,3> 

000008 

000 


DATA TV/1.0E35/ 

000009 

000 


ALTRET-. FALSE. 

000010 

000 


DO 1 1=1-, N 

000011 

000 


DO 1 J=1 ,N 

000012 

000 

1 

A(I,J)=AACI,J) 

000013 

000 


DO 2 1=1, N 

,000014 

000 


DO 2 4=1 ,N 

000015 

000 


B ( 1 , 4 ) =0 

000016 

000- 

2 

IFCI.EQ,J)B(I,J)=1.0 

000017 

000 


LIH=N-1 

000018 

000 


DO « 1=1, LIH 

000019 

000 


IF (ABSCA (I ,1) ) .LT. 1 .OE-30) CALL FLIP (A ,B , I , H , ALTRET) 

000020 

000 


L0GTMP=ALTRET 

000021 

000 


ALTRET=. FALSE. 

000022 

000 


IF(LOGTMP) GO TO 101 

000023 

000 


LIM2=I+1 

000024 

000 


00 4 4=LIH2,N 

000025 

000 


IF<TV*ABSCACI,I>).LE.ABS(ACJ ,1))) ALTRET=.TRUE. 

000026 

000 


1F(TV*ABS(A{I,I)).LE.ABS(A(4,I))) RETURN 

000027 

000 


FACTOR=A(J ,1)/A<1,I) 

000028 

000 


DO 3 K=LIH2,N 

000029 

000 

3 

AC4 ,K)=A( 4 ,K)-FACTOR*A (I ,K> 

000030 - 

000 


DO 4 K=1,N 

000031 

000 

4 

6(4 ,K)=B ( J,K)-FACTOR*B (I ,K) 

000032 

000 


DO 5 I=N,2,-1 

000033 

000 


L1M5=1-1 

000034 

000 


DO 5 4=L1M3,1 ,-1 

000035 

000 


1FCTV*ABS<A(I,1)).LE.ABS(A(J,I>)) ALTRET= .TRUE. 

000036 

000 


1F(TV*ABSCA(I,D) .LE.ABS(A(4,I))) RETURN 

000037 

000 


FACTOfi=A(4,I)/Aa,l) 

000038 

000 


DO 5 K=1,N 

000039 

000 

5 

B(4 ,K)=B(4,K>-FACT0R*B<I,K) 

00004Q 

000 


DO 6 1=1, N 

000041 

000 


FACTOR=A(I,l) 

000042 

000 


DO 6 4=1, N 

000043 

000 


I F(TV*ABS( FACTOR) .LE . A B S (S (I ,4 ) > ) A LTR ET= .TRUE . 

000044 

000 


IF (TV*ABS (FACTOR) .LE. ABS (8 (I ,4 ))) RETURN 

000045 

000 

6 

B(I,4)=B(I,4)/FACT0R 

000046 

000 


RETURN 

000047 

000 

101 

ALTRET=.TRUE . 

000048 

000 


RETURN 

000049 

. 000 


ENDS 


END ELT* 

QHDGrP ****** LATSE/MCIDAS ****** 
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****** GASORD/HCIDA S ****** 


AF^GASORBynCID AS 
RLIB62 12/22-1 7:C 
000 
000 
000 
000 

000 ^ 

000 ^ 

000 ’ 

000 ^ 

000 
000 
000 
000 
000 

000 5 

000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 

000 10 
000 

» 000 

I 000 

000 

: 000 

i 000 

^ 000 

; 000 

^ 000 

^ 000 

; 000 

> 000 

) 000 

i 000 

i 000 

5 000 


01 :51-(0») 

SUBROUTINE GASORB CRltTI |R2fT2) 

real K£,LfMfHU 

DIMENSION Rl(3)tR2(3> ^ 

CONMON/NAVCOM/NAVN tINAV , lY R, ID AT tTOTLIN fOEGLINtTOTIELf DEGELEtPlCLI 

*t4,PICELE»TKPSCL»IOYR,IODAY,T«TRlX,RlYtRlirfi1I>X»R1DYtRlDZtPlTCH»ROt 
*L,YAW,PTlM(3),THNC3),TMX(3),NLC0£F(2)tNRC0EFC2),SCLLO(2),SCLU C2), 

*£LC0EF(11,2) tSCLR0(2>f SCLRU2) ,ERC0EF(11 , 2 ) ♦ NASCEF , SCL AS 0,S CLAS1 , 
*IELEMN,IEL£HX, ASC0EF(16) 

DATA REtKEfF'U/6378.15,0.074365 74i1.0/ 

T«=T1 

DO 5 I=1t3 
R1(I)=RUI) /RE 
R2{i>=R2(I)/RE 
I CONTINUE 
R1X = R1 (1) 

R1Y=R1 (2) 

R1Z=R1{3) 

R1N0RK=SQRT (R1 (1 ) **2+R1 { 2)**2+Rl (3)**2) 
R2N0R«=SQRT(R2C1)**2+R2C2)»*2+R2(3)*»2) 

TAU=KE*CT2-T1) 

R1 R2NK=R1 N0RM*R2N0RH 

C0SAUM = (R1 (n*R2(1)+Rl(2)*R2C2) + Rl f3)*R2 <3) ) /R1 R2NB 
CSHF AN=SQRT( (1+COSANH) /2) 

t=(RlNORM+R2NORH>/(4*SeST(RlR2NM>*CSHFAN)-.5 
«=CHU*TAU**2)/ (2*SQRT(R1 R2NM)*CSHFAN)**3 
YU=1 
D YL=YU 

XL=H/YL**2-L 

CHFEAN=1-2*XL 

SHFEAN=SQlRT(4*XL*<1-XL)> 

ECAN0K=2*ATAN2 CSHFEAN,CH FEAN) 

XU = l£CANOr.-SIN CECANOH) >/SHFEAN**3 
Y0=1 +XU* CL+XL) 

IFCABS (YU-YU) .6E.1 ,0D-6)GOTO10 
YL=YU. 

A = C{TAU*SORT tHU) )/ (2*YL*SQRT (R 1R 2NH) *CS HF AN*SHFE AN ) ) **2 

F=1-A/RlN0RH*tl-C0S(ECAN0M)) 

6=TAU-A**1 .5 /SORT (HU)* ( E C ANOH-5 I N (ECANOH) ) 

R1DX=(R2(1)-F*R1(1))/G 
R10Y=(R2(2)-F*R1 C2))/G 
R1D2 = CR2 C3)-F*R1 (3))/G 
RETURN 
END 


END ELT. 

aHDGfP ****** INVERT/HClDAS ****** 



*«***« ftime/mcidas ****** 


SELT,L 

Af .FTIHE/MCIDAS 


ELT007 

RLIB62 12/22-17: 

01 :51-(0») 

000001 

000 

FUWCTXON FTIHE(ITIME) 

000002 

000 

IH=ITIME/10000 

000003 

000 

HRS=1H*60. 

000004 

000 

AHNS = FLOAT( (lTIM£-aT11EMOODO)*10000)/100) 

000005 

000 

SECS=FLOAT(«OO(ITIME»100)>/60. 

000006 

000 

FTIME=HRS-*-AMN$+SECS 

000007 

000 

RETURN 

000008 

ODD 

ENDS 

£ND ELT 

• 


SKD6,P 

****** GASORB/MCIDAS ****** 


omgwal^agbb 

OF POOE QtlAim 



«*«*** FLIP/MCIDAS ****** 


SELTfL AF .FLIP/MCI&AS 

EL7007 RL1B62 12/22-17:01:50-70,) 


000001 

000 


SUBROUTINE FLIP (A ,B ,I , N , ALTR ET ) 

000002 

000 

c 


000003 

000 

c 

THIS ROUTINE IS CALLED ONLY BY SUBROUTINE INVERT Af4D IS USED 

000004 

000 

c 

PERFORM An ELEMENTART ROU OPERATION ON MATRICES A AND 0, 

00G005 

000 

c 


000006 

000 


LOGICAL ALTReT 

000007 

000 


DIMENSION A<N, N> ,B (N,N) 

000008 

000 


altret=. false. 

000009 

000 


LIH=1+1 

000010 

000 


DO 2 L=LIW,N 

000011 

000 


IFC AB$(A{L,1)>.LT.1.0E-32)GOTO 2 

000012 

ooo 


DO 1 H=1,N 

000013 

000 


Aa,H)=A(I,M)+A(L,M) 

000014 

000 

1 

Ba,K)=B(I,M) + BtL,H) 

000015 

000 


return 

000016 

000 

2 

CONTINUE 

000017 

000 


altret=.true. 

000018 

000 


RETURN 

000019 

000 


ENOS 

END 6LT. 

3KD6,P 

****** 

fTIME 

/MCI DAS ****** 



****** FLALO/MCIOAS ****** 

aELT.L AF.FLALO/HCIOAS 

ELT007 RLIB62 1 2/22-1 7:01 :49-(0, > 

000001 000 FUNCTION FLALOCN) 

000002 000 IF (PI ,LT, 0)60 TO 1 

000003 000 N=H 

000004 000 X=1.0 

000005 000 GO TO 2 

000006 000 1 N=-H 

000007 000 -X=-1.0 

000008 000 2 FLALO=FL0AT<N/1000O)+FL0AT(tt0DCN/100,.100))/6O. + FL0ATCHOD(N,1O0))/ 

000009 000 * 3600. 

000010. 000 FLALO=FLALO*X 

000011 000 • RETURN 

000012 ' 000 ENO$ 

END, ELT. 

aHDGiP ****** FLIP ****** 





****** ES/HClOAS ****** 


9ELT,L Af.ES/HCIOAS 

ELT007 RHS62 12/22-1 7:01 :48-(0,) 

000001 000 SUBROUTINE ES (PTIME ,XL AT ,XLON, XLIN , XELE> 

000002 000 COMHON/MAV C0K/NA\«N,1NAV,IYR, 1 0 AY , TOT LIN, !> EG LI N , TOTI EL , DEGEL E j PICLI 

000003 000 *N,PICELE,T«PSCL,10YP,I0DAY,T«,R1X,R1Y,R1Z,R1DX,R1DY,R1DZ,P1TCH,R0L 

OOOOOA 000 *L,TAW,PTiM{3) ,THNC3),TMX (3) , NL COE F ( 2 > , NR COE F ( 2) , S C LLO < 2) ,SCLL1 <2) , 

000005 000 *£LCOEF(11,2),SCLftO(2),SCLR1{2>,ERCOEF(11,2) ,N ASC E F , SCL ASO , S CLAS1 , 

000006 000 *IELEMN,IELEKX, ASC0EFC1 6) 

000007 000 DATA lNIT/0/ 

000008 000 ALON=XLON 

000009 000 ALAT=XLAT 

000010 000 TIHE=PT1HE 

000011 000 I0IR=1 

000012 ' 000 ' ALINSV=0 

000013 000 AELESV=0 

000014 000 DO 100 11=1,5 

000015 000 CALL ERTOE R (ALAT, ALON , XI £R ,X2E R, X3ER, ID IR) 

000016 000 CALL ERTOST Cxi ER, X2ER , X3ER ,X1 ,X2 ,X3 ,IDIR, TIKE) 

000017 000 ' CALL STT0LVCX1 ,X2,X3,ID1R,TIHE) 

000018 000 C BC SAHL AS LV HERE, I 

000019 000 CALL BCTOPFCXI ,X2,X3,IDIR) 

000020 000 CALL P FTOTC CALI N , AELE , XI ,X 2 , X3 , 1 D IR , INIT) 

000021 000 CALL ED6C0R CPTIHE , ALl N ,0£LLIN,DELELE) 

000022 000 AELEC=AELE+DELEL£ 

000023 000 ALINC=ALIN+DELL1N 

000024 000 1FCABSCALINC-ALINSV).LE.0,5) 60 TO 150 

C00025 000 ALIWSV=ALINC 

000026 000 AELESV=AELEC 

000027 000 1SCAN=1200-C1FIXCALIN-1 ))/2 

000028 000 TIME=PT1KE+ISCAN*THPSCL 

000029 000 100 CONTINUE 

000030 000 150 CONTINUE 

000031 000 XLIN=ALINC 

000032 000 XELE=A£LEC 

000033 000 RETURN 

G00034 000 END 

END ELT, 

3HDG,P FLALO *****it 



****** ERTOST/MCIOAS ****** 


aELT,L AF.ERTOST/PICID AS 

ELT007 RL1B62 1 2/2 2-1 7:0 1 :A7-CO, ) 


000001 

000 


subroutine ERTOST (XE,rE,ZE,XtY,I,I01R»TlME) 

000002 

000 

c 

PARAMETER LIST 

000003 

000 

c 

XE,YfetZE=EARTH COORDS OF LANDMARK 

000004 

000 

c 

X,Y,Z=P0INT1N6 V£C IN INERTIAL COORDS TO LANDMARK 

000005 

000 

c 

IF 1D1R=1, EARTH TO POINTING 

000006 

000 

c 

IF IDlR=2t POINTING TO EARTH 

000007 

000 

c 

itime-hhmmss of current pointing vector 

000008 

ODD 


COMMON /NA VCON/NAVN ,INAVt IYRtID AY tTOTLIN,[>EGLlN|TOTlEL,DEGELE,PICLI 

000009 

000 


*N tPiCELEiTMPSCLtXOYR, I0DAY,TfirR1XtRlYtRUtRlDX»RlDYtRlDZ»PiTCH,R0L 

000010 

ODO 


*LtYAW,PTIM(3),TMn(3),TMX (3) ,NLC0EFC2),NRC0EF(2) tSCLL0(2 ) ,SCLl1 €2) ♦ 

000011 

000 


*ELC0EF(11,2),SCLR0(2),SCLRl(2)tERC0EF(11 , 2 ) , NAS C E F , SCL AS 0 r S CLASI , 

000012 

‘ 000 


*lELEMN,I£LEnX,ASC0EF(16) 

000013 

000 


DDD = FLOAT CHOD(1DAY,1000>> 

ODOOU 

000 


A==6378*15 

000015 

000 


B=635o.77 

000016 

000 


R1^99. 59477026 

000017 

000 


' R2=*985647336 

000018 

000 


R5=. 2506644773 

000019 

000 


PI=3. 14159265 

000020 

000 


RDPDG^PI/ISO.O 

000021 

000 


R-CR1+R2*DDD+R3*TIME) 

000022 

000 


R=AM0DCR,360* )*RDPD6 

000023 

000 


CALL ORBIT (XS ,YS, 2S jTIME) 

000024 

000 


CR=COS(R) 

000025 

ODO 


SR=S1N(R) 

000026 

000 


IFCIDIR .EQ« 2)&0 TO 10 

000027 

000 


X1=CR*XE-SR*YE 

000028 

000 


Y1=SR*Xe+CR*YE 

000029 

ODO 


21=ZE 

000030 

000 


X=X1-XS 

000031 

000 


Y=Yt-YS 

000032 

000 


2=^Z1-ZS 

000033 

000 


CALL NRHLIZ(X,Y,Z tRNORM) 

000034 

000 


return 

000035 

OOO 

10 

CONTINUE 

000036 

000 


F==B**2/A**2 

000037 

000 


AQ = F+C1.-'F)*Z**2 

000038 

000 ^ 


B0 = 2.* CCX*XS+Y*YS)*F+Z*Z$> 

000039 

OOO 


CQ=CXS**2 + YS**2)*F + ZS*^2-*B**2 

OOOOAO 

000 


RAD = BQ**2-4 #*AQ*CQ 

000041 

OOO 


IF (RAD *LT .D* ) GO TO 15 

000042 

OOO 


S=--(3Q + $QRT (RAD) )/ (2.* AQ> 

000043 

OOO 


X2 = XS-^X*S 

000044 

OOO 


Y2=Y5+y*S 

000045 

OOO 


Z2=ZS+Z*S 

000046 

OOO 


X£=CR*X2+SR*Y2 

G00047 

OOO 


YE=-SR*X2+CR*Y2 

000048 

ODO 


ZE = Z2 

000049 

OOO 


RETURN 

000050 

OOO 

15 

WRITE(6»1000) 

000051 

OOO 

1000 

F0RHATC1X,^THE SUBROUTINE ERTOST CIDIR=2) HAS^RECEiVED BAD 

000052 

OOO 


* INPUT DATA,^> 

000053 

OOO 


RETURN 

000054 

OOO 


END$ 




Ef?TOER/MCIDA S 


2ELT,L AF.ERTOER/HCIDAS 

ELT007 RLIB62 1 2/ 22--1 7: 01 :46- CO, ) 


000001 

000 


SUBROUTINE £R TOER <XL AT» XLON tX E, Y E »ZEt 10 R > 

000002 

000 


PI=3.U159265 

000003 

000 


RDPD6=PI/180.0 

000004 

000 


A=6378.15 

000005 

GOO 


B=6356.77 

000006 

000 


ESQ=(A-B>*{A+B>/A**2 

000007 

000 


IFUDk .EO. 2)GO TO 10 

000008 

000 


XET=XLAT*RDPOG 

000009 

000 


XLN=XLON*ROPDG 

0O0010 

000 


CLT=COS(XLT> 

000011 

' 000 


SLT=SIN (XLT) 

000012 

000 


CLN=COS (XLNJ 

000013 

000 


SLN=SIN(XLNJ 

000014 

000 


RR=A/S0RTt1.-ESQ*SLT**2> 

000015 

000 


XE=CLT*CLN*RR 

000016 

000 


YE=CLT*SLN*RR 

000017 

000 


Z£=SLT*RR* (1.-ESQ) 

000018 

000 


RETURN 

000019 

000 

10 

CONTINUE 

000020 

000 


IFtC (XE**2+YE**2 + ZE**2> .LT. B**2) .OR. C(XE**2+YE**2+ZE**2> 

000021 

000 


* .GT. A**2))G0 TO 15 

000022 

000 


XL0N=ATAN2(YE,XE)/R0PD6 

000023 

oco 


XLAT=ATAN (Z E / ( C 1 . -ESQ ) *S QR T ( XE ** 2+ Y E** 2 ) ) ) /R D PDG 

000024 

000 


RETURN 

000025 

000 

15 

XLAT=100. 

000026 

000 


XLON-200. 

000027 

000 


RETURN 

000028 

000 


ENDS 


END ELT. 

aHDG,P ****** ERTOST ****** 



- *★***★ 

ED6FT/HCIDAS 

ifk^*** 


000113 

000 

706 

F0R«ATC2X,'DI«ENSI0N=',I3) 


000114 

000 


CALL TQ(MOUT) 


000115 

000 


IX=IRES*(IRES-1>/2 


000116 

000 


DO 470 1=1, IRES 


000117 

000 


COEF (1)=W0RK (1+lX) 


000116 

000 


ENCODEa32,705,HOUni ,COEFCI) 


000119 

000 


CALL TQ(HOUT) 


000120 

00 0 

470 

ERCOEF(I,IPAIR)=COEF<I) 


000121 

000 

490 

CONTINUE 


000122 

000 


TKNCIPAIR+1 )=AMAX1 (THNL,THNR) 


000123 

000 

500 

TMX CIPAIR + 1 )=AKIN1 (THXL,THXR> 


000124 

' 000 

CONTINUE 


000125* 

000 


CALL WCOH 


000126 

000 


CALL EMESS(3HFIN,0) 


000127 ■ 

000 


RETURN 


000128 

000 

990 

CONTINUE 


000129 

000 


CALL TQC72H ERROR RETURN FROM APCH, APFS. 


000130 

000 


* 


000131 

000 


RETURN 


000132 

000 


ENDS 


000133 

000 

SFILEMA 


000134 

000 

DELETE 6CCGCC,G0RP 


000135 

000 

SINCLUDE ATSSSP 


000156 

000 

SCATALOG 


C00137 

000 

NAME= 

GCCGCC ,5,R,U,D 


000136 

000' 

TYPE= 

F6 


000139 

000 

LlD^ATSFLBtLL 


000140 

000 

BEGIN 



000141 

000 

$E0J 



END ELT. 
OHDGtP 

ic-kiz-irk* 

ERTOER 






original 
OF POOR 

page is 

qUALTPY 





EDGFT/MCIDAS 




000056 

000 


ISCNL(NU=1200-(lY1-1)/2 

000057 

000 


IUL(NL)=JU1N(9,JVEC5*(T2-T1) /60.0 

000058 

000 


60 TO .100 

000059 

000 

200 

CONTINUE 

000060 

000 


NR=NR+1 

000061 

000 


IFCnR.GT.300> 60 TO 100 

000062 

000 


1SCNRCNR)=1 200-(IYl-1>/2 

000063 

000 


1UR(NR>=JWINC9,JVEC)*(T2-T1) /60.0 

000064 

* 000 


60 TO 100 

000065 

000 

400 

CONTINUE 

000066 

ooo' 


1FCNL.lt. IP+D'GO TO 459 

000067 

000 


DO 410 1=1, NL 

000068 

000 


TS=PTIKClPAIfi+1>+ISCNLCI>*THPSCL 

000069 

000 


TMNL=AHIN1 CTKNLtTS) 

000070 

000 


THXL = Al<tAXl CTHXLfTS) 

000071 

000 


DATI CI)=ISCNLCI> 

000072 

000 

410 

DATI CI+NL>= lUL Cl)/ 100.0 

000073 

000. 


DATI C2*NL+1 )=-1 .0 

000074 

000 


CALL apchcdati,nl,ip,xo,xo,work,ier> 

000075 

000 


IfClER.LT.O) GO to 990 

000076 

000 


SCLL0C1PAIR)=X0 

000077 

000 


S CLL1 C1PA1R>=XD 

000076 

000 


EPS=1 .OE-4 

000079 

000 


I0P = + 1 

000080 

000 


•ETA=1 .OE-3 

000081 

000 


CALL APFS(W0RK,IP,1RES,I0P,EPS,ETA,IER> 

000082 

000 


IFCIER.lt. 0) GO TO 990 

000083 

000 


NLC0EFCIPAIR)=IRES 

000084 

000 


ENCODE Cl 32, 706,M0UT)IRES 

000085 

boo 


CALL T8CN0UT) 

000086 

000 


IX=lRES*ClRES-1>/2 

000087 

000 


DO 420 1=1 ,IRES 

000088 

000 


CO£FCl)=kORK(I+IX) 

000089 

000 


ENCODE Cl 32,70S,KOU7)1,COEF Cl> 

000090 

000 

705 

F ORH AT C2X, 'COE FC', 12, '> = •'•, £20.9) 

000091 

000 


CALL TG(HOUT) 

000092 

000 

420 

ELC0EFC1,1PAIR)=C0‘EFCI) 

000093 

000 

459 

CONTINUE 

000094 

000 


1FCNR.lt. IP+1> GO TO 490 

000095 

000 


DO 460 1=1, NR 

000096 

000 


TS=PTIMClPAlfi + 1 )+lSCNR CI)*THPSCL 

000097 

000 


THNR=AKINl CTMNR,TS) 

000098 

000 


TKXR=AKAX1 CTMXR ,TS) 

000099 

000 


DATI CI)=ISCNRCI> 

000100 

000 

460 

DATI CI+NR)=IUR Cl)/ 100.0 

000101 

000 


DATI C2*NR+1>=-1.0 

000102 

000 


CALL APCH CDATI ,NR , IP,XD ,XO,WORK,IER) 

000103 

000 


IFCIER.LT.O) GO TO 990 

000104 

000 


SCLRO C1PAIR)=X0 

000105 

000 


SCLRI C1PA1R)=X0 

000106 

000 


EPS=1 .oe-4 

000107 

000 


10P=+1 

000108 

000 


ETA=1 .OE-3 

000109 

000 


CALL APFSCW0RK,IP,IRES,10P,EPS,ETA,1ER) 

000110 

000 


IFCIER.LT.O) GO TO 990 

000111 

000 


' NRCOEF CIPA1R) = IR£S 

000112 

000 


ENCODE Cl 32, 706, KOUDIRES 



ifir-tc-irit it 

ED6FT/MCIDAS 


aELTtl 

AF,e06FT/«ClDAS 


ELT007 

RLIB62 12/22-17: 

01:44-C0r) 

000001 

000 

SJOB 

EOGFIT U3200 

000002 

000 

SOPTION .Bf9,15»20 

000003 

000 

SFORTRAN 

000004 

000 


subroutine main 

000005 

000 


DIMENSION JWlNm,10) 

000006 

000 


DIMENSION OAT I (601 ) ,U0RK(231 ) ,C0EFC21) 

000007 

000 


DIMENSION ISCNL(300),ISCNR (300) , lUL (300 > , lUR (300 ) ,NOUT ( 1 7) 

000006 

000 


DIMENSION MINCIO) 

000009 

000 


DIMENSION nOUT(4A) 

000010 

000 


COHMON/NAVCOH/NAVN,IHAV, lYR, IDAY,TOTLIN,DEGLIN,TOTIEL,DEGELE,PICLI 

000011 

000 


*N,PIC£L£,T«PSCL,IOYR,IODAY,TH,R1X,RlY,R1ZtR10X,R1DY,RlDZ,PITCHtROL 

000012 

000 


*L,YAU,PT1K(3) ,TMS(3>,THX (3) , NL COE F ( 2 > , NRC OE F (2) , SCLLO ( 2 ) , SC LLl (2>, 

000013 

000 


*ELC0EF(11 ,2) ,SCLR0(2) ,SCLR1(2) ,E R C OE F ( 1 1 ,2) , NAS CEF ,SCL AS 0,S CLASI < 

000014 

000 


*IELEHN,IELEHX,aSC0EF(16) 

000015 

000 


DATA THNL/UA0.0/,THXL/0 .0/,THNfi/14 40.0/,THXR/0.0/ 

000016 

000- 


data LUK/1 0/ , JSECSV/-1 / .MUMPER /I 0/ .NUHENT/1 1 / ,J OUT/2/ 

000017 

000 


DATA M1N/6H6CC6CC,8*0/ 

000016 

000 


CALL IQ(HIH) 

000019 

000 


IP=11 

000020 

000 


CAUL RCOH 

000021 

000 


ENCODE (132, 701 ,HOUT)PTIH 

000022 

000 

701 

FORMAT(5X,3F20.6) 

000025 

000 


CALL TP(1,H0UT) 

000024 

000 


CALL VARAS6 {'WIN'.LUN) 

000025 

000 


CALL OPN(LUN) 

000026 

000 


DO 500 IPAlR=1f2 

000027 

000 


NL*=0 

000028 

000 


NR=0 " 

000029 

000 


THNL=UAO.O 

000030 

000 


TMXL=0.0 

000031 

000 


THNR=1440.0 

000032 

000 


THXR=0.0 

000033 

000 


CALL SCRA(LUNtO) 

000034 

000 


IK=0 

000035 

000 

90 

CONTINUE 

000036 

000 

100 

CONTINUE 

000037 

000 


JSEC=1K/NUHPER 

000038 

000 


JVEC=HOD (IK ,NUHPER)+1 

000039 

000 


IF ( JSEC.NE . JSECSV) CALL RE AD W C LUN , 1 1 0 , J UIN) 

000040 

000 


IF(JSEC.NE.JSBCSV) JSECSV=JSEC 

000041 

000 


IK=IK+1 

000042 

000 


IF(JHIN(1 , JVbC) . to .0) GO TO 400 

000043 

000 


IH»M0D(JVJlN(4,JVtC) ,10) 

000044 

000 


IF(lh,GT.4) bO TO 100 

000045 

000 


T1=FTIME(JWIN(2,JVeC)) 

000046 

000 


T2=FTIHE (JWIN(3 , JVEC) ) 

000047 

000 


IF(T1 .NE.PTlMd ).0R. T2.NE.PTIH(1PAIR+1))G0 TO 100 

000048 

000 


IF(MOD(JHIN(4,JVEC)/100, 1000). NE. 100) 60 TO 100 

000049 

000 


I Y1 =JWIN (5 rJVEC) 

000050 

000 


lY2=JWIN(7,JVeC) 

000051 

000 


IFdYl .NE. 1Y2) 60 TO 100 

000052 

000 


IX1=JWiN(6,JVEC) 

000053 

000 


1FUX1 .GE.PICELE) GOTO 200 

000054 

000 


NL=NL+1 

000055 

000 


IF(NL.6T.300) 60 TO 100 IS 



****** EDGCOR/HCIDA S ****** 


aELT,L AF.EDGCOR/HCrO AS 
ELT007 RLIB62 12/22-1 7:01 

000001 000 SUBROUTINE EDGCOR (PTXHE , ALIN , D ELtIN lOELELE) 

000002 00 0 COHMON/NAVCOM/NAVN, INAVt IVRt ID AY , TOTLIN , D EGLIN, TOTI EL r PEGELE ,PICLI 

000003 000 *N,P1CEL£,THP-SCL,10YR,I0D AY ,TK,R1X,R1Y,R12,R1DX,R1DY,R1 DZ , PITCH fROL 

000004 000 *LtYAH,PTIH C3) ,T«NC3),THX C3 ) , NLCO E F ( 2) ,NRC OE F (2) , SCLL0(2> ,SCLLt <2) » 

000005 000 *ELC0EF(11 ,2 ) t S CLRO ( 2) , SC LR1 ( 2) ,£RCOEF (1 1 , 2) ,NASCEF ,SCLAS0,SCLAS1 , 

000006 OOO *1EL£HN,1ELEHX,ASC0EF<16> 

000007 poo DATA A/522.07/,LC/1235/ 

OOOOOS boo 1SCAN=1200-<IFIX(ALIN-1)>/2 

000009 000 TIHE=PTI«E+ISCAN*THPSCL 

000010 000 DO 20 XC=2,3 

000011 000 20 IFCTHNdO.LE. TIME, AND. TIKE. LE.THXCIO) GO TO 100 

000012 000 DELL1N=0.0 

000013 000 DELELE=0.0 

000014 000 RETURN 

000015 000 100 CONTINUE 

000016 000 iX=lC-l 

000017 000 XL = ISCAN*SCLL1 (IX)+SCULOCIX) 

000018 000 CALL CAPS CEL, X L , EL COE F ( 1 , IX) , MLCOE F (I X) > 

000019 000 XR=ISCAN*SCLR1 (IX)+SCLRO(IX) 

000020 000 CALL CNPS (ER , XR , ERCOE F (1 , IX) ,NRCOEF UX) ) 

000021 000 DEL£LE=(EL+ER)/2.0 

000022 000 KSC=1200-(LC-1 )/2 

000023 OOO Y=ISCAN-hSC 

000024 OOO XX=SQRTCA**2-Y**2) 

000025 OOO DELTA=IER-EL)/2.0 

000026 OOO RAD=Y**2-2*XX*0ELTA-DELTA**2 

00002? OOO DELLIN=SQRT<RAD)-Y 

000026 OOO RETURN 

000029 OOO ENDS 

END ELT. 

SHOG.P ****** EDGFT/flCIDAS ****** 



****** EATOST/MCIDAS *****<^ 


SELTtL Af .EATOST/MCIDAS 

ELTOOr RLI862 1 2 / Z 2-1 7:01 ; A3- (0» ) 


000001 

000 


SUaROUTlNE EATOSTCPICTIM,LNS»IES,LA,L0,lUNK,INAV,BETAfBETDOTf AF) 

000002 

000 


PTlKE=PrCTIH 

000003 

000 


IF(IUNK.Ea.l) GO TO 1 

00000<i 

000 


IfaUNK.EO.2) GO TO 2 

000005 

000 


RETURN 

000006 . 

000 

1 

CONTINUE 

000007 

-000 


XLIN^LNS 

00000& 

000 


XELE=IES 

000009 

000 


CALL S£(PTIME,XLiN,XELE»XLAT»XLOM) 

000010 

000 


LA = 1LAL0 CXLAT) 

000011 

000 


LO=1LALO (XLOHJ 

000012 

000 


RETURN 

000013 

000 

2 

CONTINUE 

oooou 

000 


XLAT=FLALO(LA) 

00001S 

000 


XLON=FLALO(LO) 

000016 

000 


CALL ESCPTIHE, XLAT,XLON,XLIN,XELE) 

000017 

000 


LNS=XLIN 

000018 

000 


IES=XELE 

000019 

000 


RETURN 

000020 

000 


END 


END ELT^ 

aHDGtP ****** EDGCOR/HCIDAS ****** 





****** STTOLV/HCIDAS ****** 


3ELT,L AF .STTOLV/HCIB AS 

ELT007 RLIB62 12/22-1 7:02 :08-{0, ) 

000001 000 SUBROUTINE STTOLV (X , Y , Z il DI R , TIME > 

000002 000 C If IDIR=1, POINTING VECTOR (X,T,Z) IS TRANSFORMED 

000003 000 C FROM SAT INERTIAL TO LOCAL VERTICAL FRAME. 

000004 OOO C IF I0IR=2, POINTING VECTOR (X,Y,Z) IS TRANSfORHEO FROR 

000005 ODD C LOCAL VERTICAL TO SAT INERTIAL FRAOE, 

000006 000 CALL ORB IT C XS , YS , Z S .TIKE ) 

000007 000 CALL NRHLIZ (XS .YS,ZS,XNORH) 

OOOOOS 000 X1=X 

000009 000 Y1»Y 

000010 000 Z1=Z 

000011 000 D=SQRT(XS**2+YS**21 

000012 000 IF CID1R.EQ.2J 60 TO 10 

000013 000 X=(-YS*Xl*XS*Yl)/0 

000014 000 Y=(XS*ZS*X1+YS*ZS*Y1-Z1*»**2)7D 

000015 000 Z=-(XS«X1+YS*Y1'»-ZS»Z1) 

000016 000. RETURN 

000017 000 10 X=-YS*Xl/D+XS*ZS*Y1/0-XS*21 

000016 000 r=XS*X1/0+YS*ZS*Y1 /D-YS*Z1 

000019 000 Z=-D*Y1-2S*Z1 

000020 000 RETURN 

000021 000 ENDS 

END ELT. 

3H06.P ****** UNIT/ KCIDAS ****** 



UNIT/HCIOAS ****** 


aELT,L AF.UNIT/NCIDAS 


ELT007 

RL1B62 12/22- 

000001 

000 

000002 

000 C 

000003 

000 C 

000004 

000 

000005 

000 

000006 

000 

000007 

000 

000008 

000 

000009 

000 

000010 

000 

000011 

000 


1 7:02:09-(0,) 

SUBROUTINE UNIT(A> 
THIS ROUTINE RETURNS IN 

DIMENSION AC9) pBC9> 
DATA B/1 40,0. 0,0. 0» 

^ 0 . 0 , 1 . 0 , 0 . 0 , 

* 0. 0,0.0, 1.0A 

DO 10 1=1,9 
10 ACD^BCl) 

RETURN 

ENDS 


A 3 X 3 IDENTITY MATRIX 


* END ELI. 
cvHDG,N 
SPIN 


F 



